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HOW TO USE THIS MANUAL 



This manual is a programmer's guide to the 8087 Support Library. The Library is a 
package of software tools that enable the 8087 or the full 8087 emulator to provide a 
full range of floating point capabilities to ASM-86 and PL/M-86 programmers. 

This manual is primarily a reference guide, to be used when you are writing code for 
ASM-86 and PL/M-86 programs. We have made the documentation for each pro- 
cedure and function as self-contained as practical, to avoid excessive page-flipping. 

To obtain an initial overview of the library, read Chapters 1 and 2 and the first parts 
of Chapters 3, 4, and 5. Read Chapter 6 if you are concerned with compliance to the 
IEEE Standard. Familiarize yourself with the contents of the appendices. 

To find specific facts about individual procedures and functions, look in the 
reference pages, headed by large names, at the ends of Chapters 3, 4, and 5. 

Note that the reference pages contain much material discussing unusual values such 
as denormals, unnormals, pseudo-zeroes, etc. You should remember that these 
values are very uncommon. For most applications they never occur, and you may 
ignore the paragraphs that discuss them. 

For quick reference when you are an experienced user of the 8087 Support Library, 
consult Appendices D and E. 



Related Publications 

The following manuals will be useful to you when you are reading this manual. Par- 
ticularly crucial is the documentation on the 8087 itself, to be found in the Assembly 
Language Reference Manual and the 8086 User's Manual Numerics Supplement. 
These manuals all describe the Series-Ill versions (8086-based) of the software 
products. 

• iAPX 86, 88 Family Utilities User's Guide, Order Number 121616 

• PL/M-86 User's Guide, Order Number 121636 

• 8086/8087/8088 Macro Assembly Language Reference Manual, Order Number 
121627 

• 8086/8087/8088 Macro Assembler Operating Instructions, Order Number 
121628 

• The 8086 Family User's Manual Numerics Supplement, Order Number 1 21 586 



Nomenclature for the 8086 Family 

Three naming conventions are used in this manual to describe the 8086 family of 
microprocessors: 

1 . The iAPX Convention describes the system containing one or more components 
of the family. It must contain an 8086 or 8088 microprocessor. The table below 
shows the components used in the most common iAPX systems. 

2. The 8086/8087/8088/8089 Convention describes a specific component within 
an iAPX system. 
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3. The NDP/NPX Convention describes a numeric functional unit which can be 
realized by more than one hardware configuration. 

• NPX (Numeric Processor Extension) refers either to the 8087 component or 
the software 8087 full emulator. 

• NDP (Numeric Data Processor) refers to any iAPX system which contains 
an NPX. 
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8086 
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iAPX 88/20 
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iAPX 88/21 
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iAPX 88/22 




1 




2 



Notational Conventions 



The programming examples given in this manual are presented as chunks of code 
that can be inserted intact into a PL/M-86 or ASM-86 program. They contain no 
formulas of general syntax. 

The examples of LINK86 and LOC86 invocations follow these conventions common 
to Intel manuals: 

• EIH33 background is used to denote user input to the operating system. 

• <cr> indicates a carriage return. 
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CHAPTER 1 
INTRODUCTION 



The 8087 Support Library greatly facilitates the use of floating point calculations in 
ASM-86 and PL/M-86. It adds to these languages many of the functions that are 
built into applications programming languages, such as Pascal-86 and 
FORTRAN-86. 

The full 8087 emulator, E8087, and the interface libraries E8087.LIB, 8087. LIB, 
and 87NULL.LIB, allow a choice of run-time environments for modules containing 
floating point code. The emulator is the functional and programming equivalent to 
the 8087, implemented in 16K bytes of 8086 software. 

The decimal conversion module, DCON87.LIB, aids the translation between 
human-readable numbers of many decimal formats, and binary numbers in all the 
formats supported by the NPX. It also provides translation between binary formats. 

The common elementary function library, CEL87.LIB, provides a useful assortment 
of transcendental, rounding, and other common functions involving floating point 
numbers. 

The error handler module, EH87.LIB, makes it easy for you to write interrupt pro- 
cedures which recover from floating point error conditions. 



Why the Support Library Is Needed 

There are three areas in which many applications will require support software to 
overcome functional limitations of the NDP. 

First, the NDP is exclusively a binary arithmetic machine. Although the NDP 
manipulates a variety of data types, it cannot operate directly on floating-point 
decimal quantities. Software must be invoked to convert decimal inputs to binary 
formats, and binary answers to decimal output. 

Second, the 8087 instruction set is limited in its implementation of elementary 
functions. There is a tangent function, for example, but no sine, and the tangent is 
defined only for inputs between and tt/4. The implementation of a full set of 
elementary functions with complete domains is left to software. 

Third, the handling of unmasked exceptions can complicate application software. 
Although the masked mode of error detection and recovery is adequate for most 
applications, there is sometimes need for customized error recovery. Deciphering 
NPX status information and restoring the NPX to a recovery state can be quite 
difficult. The task can be made easier with software. 

The 8087 Support Library provides the software necessary to eliminate the above 
deficiencies. The Support Library makes the NDP a complete, accurate, easy-to-use 
floating point computer. 



System Software Needed to Use the Library 

Since the Support Library is called by ASM-86 and PL/M-86 programs, you will 
need either the ASM-86 assembler or the PL/M-86 compiler. 



1-1 



Introduction 8087 Support Library 



You will need to compile your PL/M-86 modules under the MEDIUM or LARGE 
models of computation, since all of the procedures of the Support Library are FAR 
procedures. 

You will also need LINK86 and LOC86 to convert the object modules generated by 
your programs into programs that can be loaded or programmed into memory (or 
programmed into ROM) and executed. 

All of the above programs come in versions that run on an 8080/8085-based 
development system under ISIS-II. They also come in versions that run on an 8086- 
based development system via the Series III RUN utility. 

For consistency, we give all examples of LINK86 invocations as they appear running 
the 8080/8085 version of LINK86. LINK86 and the Support Library are assumed to 
be on drive 0; all user-supplied files are assumed to be on drive 1 . 

No run-time system software is needed for the Support Library. Once you have 
obtained the final executable program, the program will run on any NDP system. 
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CHAPTER 2 

THE FULL 8087 EMULATOR 

AND INTERFACE LIBRARIES 



Overview 

The full 8087 emulator is a software module, executable on an iAPX86/10 or 
iAPX88/10, which completely and exactly duplicates all the functions of the 8087. 
You can use the emulator to develop prototypes for systems that will have an 8087; 
or you can use the emulator if you need the 8087's floating point capabilities but do 
not need its speed,, The emulator is contained in the file E8087, which you include in 
a LINK86 statement with your program modules. 

Included with the emulator are three interface libraries, E8087.LIB, 8087. LIB and 
87NULL.LIB. By selecting one of these three libraries to appear in your LINK86 
statement, you choose the run-time environment of your program: running with the 
emulator, with the 8087 component, or with no NPX capability. 

The interface libraries also contain the procedures INIT87 and INITFP, which 
initialize the NPX under the various environments. 

The full 8087 emulator is an expanded update of the partial 8087 emulator, PE8087, 
which we offer with PL/M-86. If you are adding E8087 to a LINK86 statement that 
includes PE8087, you must eliminate PE8087 from the statement. 

Whenever we mention "the emulator" in this manual, we are referring to the full 
8087 emulator. 



How to Use the Emulator in Your Programs 

To invoke the emulator in your source programs, you simply issue instructions as if 
you had an 8087. There is absolutely no difference between source programs that use 
the emulator and source programs that use the component. 

Furthermore, there is no difference between unlinked object modules (or libraries) 
that use the emulator, and those that use the component. CEL87.LIB, 
DCON87.LIB and EH87.LIB can all be made to run either with the emulator or 
with the 8087 component. 

The changes that you make in switching between the emulator and the component 
occur entirely in the invocation of LINK86. The files E8087 and E8087.LIB are 
linked to your object modules when you are using the emulator; the one file 
8087. LIB is linked when you are using the component. 

The source and object module compatibility between the emulator and the compo- 
nent is accomplished by a set of special external variables. These variables are issued 
automatically by all Intel translators that handle 8087 code, and are invisible to the 
programmer except in some LOC86 listings. The names of the variables are given in 
Appendix B. 



Emulation of Unmasked Exceptions 

When the 8087 detects an error, it checks its control word to see if the error is 
masked or unmasked. If unmasked, an interrupt signal is generated on one of the 
output pins of the 8087. 
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The emulator, in the same unmasked error situation, makes a call to interrupt 16 
decimal. Thus, an INTERRUPT procedure is executed; the address of the procedure 
is given in a four-byte pointer at memory location 00040 hex. 

Note that with the 8087 component you can configure your system and program 
your interrupt controller so that exactly the same interrupt is generated when the 
8087 sees an unmasked error. Thus, the source code even for your exception handler 
can remain exactly the same. 



Emulator Memory and Stack Requirements 

The emulator, of course, occupies 8086 memory that is not needed when the 8087 
component is used. The memory used comes in four segments: 

1 . The code of the emulator, which can be located in read-only memory, is in a 
segment called aqmCODE. It occupies about 16K bytes. 

2. The emulator uses a non-reentrant data block, which must be in RAM. The 
block is about 150 bytes long, and is in a segment called aqmDATA. 

3. The emulator uses 8086 software interrupt 16 decimal, and interrupts 20 
through 31 decimal. The addresses of these interrupts are stored in memory 
locations 00040 hex through 00043 hex, and 00050 hex through 0007F hex, by 
INIT87 or INITFP. While you are still using the emulator, these locations 
cannot be changed. 

4. The emulator uses about 60 bytes of the same 8086 stack as used by the calling 
routine. This stack must be in a segment called STACK. 



How to Initialize the NPX 

When power is first applied to an NDP system, the NPX must be brought to a 
known state. This is accomplished by calling one of the 8087 initialization pro- 
cedures. If the emulator is being used, the addresses of interrupts 20 through 31 
must be loaded before any floating point instructions are executed. 

We have provided two 8087 initialization procedures, INIT87 and INITFP, in each 
of the three interface libraries. The 8087. LIB versions initialize the 8087 component 
and set the 8087 control word. The E8087.LIB versions load interrupts 20 through 
31 before the 8087 initialization. The 87NULL.LIB versions satisfy 8087 externals 
with null procedures. 87NULL.LIB should only be used with software that contains 
no NPX instructions. 

The difference between INIT87 and INITFP is the 8087 control word default 
settings. INIT87 masks all exceptions. INITFP masks all exceptions except "I". 
Although you may later modify the control word, you should first call either INIT87 
or INITFP. 



PL/M-86 Usage of INIT87 or INITFP 

Following is the declaration of 1NIT87, which must appear at the top of your 
initialization module, and the call, which must be invoked before any floating point 
arithmetic is performed. The statements for INITFP are the same, with INITFP 
instead of INIT87. 

Since INIT87 and INITFP are FAR procedures, you must compile the PL/M-86 
module under the MEDIUM or LARGE controls. 
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IN IT87 : PROCEDURE EXTERNAL; 
END; 

CALL INIT87; 



ASM-86 Usage of INIT87 or INITFP 

Following is the declaration of INIT87, which must appear at the top of your 
initialization module, and the call instruction, which must appear within the 
initialization procedure. The statements for INITFP are the same, with INITFP 
instead of INIT87. 

; The following line must appear outside of all SEGMENT- 
; ENDS pairs 

EXTRN INIT87: FAR 



CALL INIT87 



; Set up the NPX 



Linkage of the Emulator and the Interface 
Libraries to User Programs 

If your target system has an 8087 component, link the library 8087. LIB to your 
object modules. Neither the emulator (E8087) nor the partial emulator (PE8087) 
should appear in the LINK86 command. 

If your target system does not have an 8087 component, link the emulator E8087 
and the interface library E8087.LIB to your object modules. The partial emulator 
PE8087 should not appear in the LINK86 command. 

If you issue a call to INIT87 or INITFP but do not have any other floating point 
instructions, you should link 87NULL.LIB to your object modules. 

The module E8087 may appear anywhere in the sequence of input modules to 
LINK86. The interface libraries must appear after all modules that contain 8087 
code. 

Following is the suggested order for the object modules in your LINK86 command. 

Your object modules 

DCON87.LIB if you are using it 

CEL87.LIB if you are using it 

EH87.LIB if you are using it 

8087. LIB if you are using the component, or 

E8087,E8087.LIB if you are using the emulator. 



Examples: Suppose your program consists of two modules, MYMODl.OBJ 
and MYMOD2.0BJ, created with either ASM-86 or PL/M-86. If the program is 
to be executed on a system using the emulator, issue the command 



LINK86 :F1 :MYM0D1 .OBJ, : F1 : MYM0D2 . OBJ , &<cr> 
3 :F0:E8087,; : FO : E8087 . LI B TO : F1 : MYPROG . LNK 



To modify the program so that it runs on a system using the 8087 component, you 
do not have to change your source files. Simply change the LINK86 command to 
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>> 



LINK86 :F1 :MYM001 .OBJ, : F1 : M Y MO D 2 . OBJ , &<cr> 
: FO : 8087 .LIB TO :F1 :MYPR0G . INK 



If your program calls INIT87 or INITFP, but contains no other floating point 
instructions, issue the command 



>> 



LINK86 :F1 :MYM0D1 .OBJ, : F1 : M Y M|0 D 2 . OBJ , &<cr> 
: FO : 87NULL .LIB TO : F1 : MYPROG . LNK 
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CHAPTER 3 
THE DECIMAL CONVERSION LIBRARY 



This chapter describes the use of DCON87.LIB, a library of procedures which con- 
vert floating point numbers from one format of storage to another. 



Overview of DCON87 Procedures 

Intel supports three different binary formats for the internal representation of 
floating point numbers (see Appendix C). DCON87.LIB provides translation 
between these formats and an ASCII-encoded string of decimal digits. 

The binary-to-decimal procedure mqcBIN DECLOW accepts a binary number in 

any of the three formats. Because there are so many output formats for floating 
point numbers, mqcBIN DECLOW does not attempt to provide a finished, for- 
matted text string. Instead, it provides the "building blocks" for you to use to con- 
struct the output string which meets your exact format specification. 

The decimal-to-binary procedure mqcDEC BIN accepts a text string which con- 
sists of a decimal number with optional sign, decimal point, and/or power-of-ten 
exponent. It translates the string into the caller's choice of binary formats. 

An alternate decimal-to-binary procedure mqcDECLOW BIN provides an input 

format similar to the output of mqcBIN DECLOW. The low-level input interface 

is provided for callers who have already broken the decimal number into its consti- 
tuent parts. 

The procedures mqcLONG__TEMP, mqcSHORT_TEMP, mqcTEMP_LONG, 

and mqcTEMP SHORT convert floating point numbers between the longest 

binary format, TEMP REAL,' and the shorter formats. The conversions are per- 
formed so that the outputs are NaN's if and only if the inputs are NaN's. This 
improves upon the conversion algorithm used by the 8087 when it loads shorter for- 
mats to and from its stack. 

All DCON87 procedure names begin with the three letters "mqc". We have added 
this prefix to reduce the chance that a DCON87 name will conflict with another 
name in your program. To make the names more readable, we have put the "mqc" 
prefix in lower case letters" throughout this chapter. 

DCON87 contains alternate public names for its procedures which are used by some 
Intel translators. These names are listed in Appendix F. 



Declaring DCON87 Procedures in Your 
ASM-86 Programs 

In each source program module which calls a DCON87 procedure, you must declare 
that procedure as external before it is used. 

Following are the declarations you should make at the top of your ASM-86 pro- 
gram. You should include only those procedures you use. 

All DCON87 procedures must be declared FAR. 
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This ASM-86 code also includes declarations which ease the construction of 
parameters. The ASM-86 examples given with the procedures later in this chapter 
assume that these declarations have been made. 

To make reference easier, we have duplicated the procedure declarations in the 
examples for the individual functions. You should not duplicate the declarations in 
your programs. 

; ASM-86 declarations for using the DCON decimal 
; conversion library 

; the following EXTRN statements must appear outside of 
; all SEGMENT-ENDS pai rs. 

EXTRN mqcBIN_DECLOW: FAR 

EXTRN mqcDEC_BIN: FAR 

EXTRN mqcDECLOw_BIN: FAR 

EXTRN mqcLONG_TEMP: FAR 

EXTRN mqcSHORT_TEMP: FAR 

EXTRN mqcTEMP_LONG: FAR 

EXTRN mqcTEMP_SHORT: FAR 

SHORT_REAL EQU 
L0NG_REAL EQU 2 
TEMP_REAL EQU 3 

BIN_DECLOW_BLOCK STRUC 

BIN_PTR DD ? 

BIN_TYPE DB ? 

DEC_LENGTH DB ? 

DEC_PTR DD ? 

DEC_EXPONENT DW ? 

DEC_S IGN DB ? 
BIN_DECLOW_BLOCK ENDS 

DEC_BIN_BLOCK STRUC 

DD ? ; The names of these fields are the same as in 

DB ? ; BIN_DECL0W_BL0CK. 

DB ? 

DD ? 
DEC_BIN_BLOCK ENDS 

DEC L0W_BIN_BL0CK STRUC 

DD ? ; The names of these fields are the 

DB ? ; same as in B IN_DEC L0W_B LOCK 

DB ? 

DD ? 

DW ? 

DB ? 
DECLOW BIN BLOCK ENDS 



Declaring DCON87 Procedures in Your 
PL/M-86 Programs 

Following are the declarations you should make at the top of your PL/M-86 pro- 
gram. You should include only those procedures you use. 
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This PL/M-86 code also includes declarations which ease the construction of 
parameters. The PL/M-86 examples given with the procedures later in this chapter 
assume that these declarations have been made. 



To make reference easier, we have duplicated the procedure declarations in the 
examples for the individual functions. You should not duplicate the declarations in 
your programs. 



To use DCON87.LIB, you must compile your PL/M-86 modules under the 
MEDIUM or LARGE models of computation, since DCON's procedures are FAR 
procedures. 



See the section entitled "An Alternate Method of DCON87 Error Detection" later 
in this chapter for a discussion of when you might want to declare the procedures as 
BYTE procedures (to receive an error status). 



/* PL/M-86 declarations for using the DCON decimal 
conversion library * / 

mqcBIN_DECLOW: PROCEDURE (BLOCK PTR) EXTERNAL; 

DECLARE BL0CK_PTR POINTER; 
END mqcBIN_DECLOW; 

mqcDEC_BIN: PROCEDURE (BL0CK_PTR) EXTERNAL; 

DECLARE BL0CK_PTR POINTER; 
END mqcDEC_BIN; 

mqcDECLOW_BIN: PROCEDURE (BL0CK_PTR) EXTERNAL; 

DECLARE BL0CK_PTR POINTER; 
END mqcDECLOW_BIN; 

mqcLONG_TEMP:: PROCEDURE ( LONG_REAL_PTR , TEMP_REAL_PTR) 

EXTERNAL; 

DECLARE (LONG_REAL_PTR,TEMP_REAL_PTR) POINTER; 
END mqcLONGJTEMP; 

mqcSHORT_TEMP: PROCEDURE (SHORT REAL_PTR , TEMP_REAL_PTR) 

EXTERNAL; 

DECLARE (SHORT_REAL_PTR,TEMP_REAL_PTR) POINTER; 
END mqcSHORT__TEMP; 

mqcTEMP_LONG:: PROCEDURE (TEMP_REAL_PTR , LONG_REAL_PTR) 

EXTERNAL; 

DECLARE (TEMP_REAL_PTR,LONG_REAL_PTR) POINTER; 
END mqcTEMP_LONG; 

mqcTEMP_SHORT: PROCEDURE (TEMP_REAL_PTR , SHORT_REAL_PTR) 

EXTERNAL; 

DECLARE (TEMP_REAL_PTR,SHORT_REAL_PTR) POINTER; 
END mqcTEMP_SHORT; 

DECLARE SHORT_REAL LITERALLY '0'; 
DECLARE LONG REAL LITERALLY '2';: 
DECLARE TEMP REAL LITERALLY , 3';: 
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DECLARE BIN_DECLOW_BLOCK STRUCTUREC 
BIN_PTR POINTER, 
BIN_TYPE BYTE, 
DEC_LENGTH BYTE, 
DEC_PTR POINTER, 
DEC_EXPONENT INTEGER, 
DEC_SIGN BYTE); 

DECLARE DEC_BIN_BLOCK STRUCTURE ( 
BIN_PTR POINTER, 
BIN_TYPE BYTE, 
DEC_LENGTH BYTE, 
DEC PTR POINTER); 



Stack Usage of DCON87 Procedures 

DCON87 requires 176 bytes of the 8086 stack for its internal storage. This amount is 
allocated in the public STACK segment by the DCON87 modules. 

All DCON87 procedures save the entire state of the 8087 upon entry and restore it 
upon exit. 



The Accuracy of the Decimal Conversion Library 

DCON87 guarantees that an 18-digit decimal number, when converted into 

TEMP REAL and then back to decimal, will remain the same. However, DCON87 

cannot guarantee that a TEMP REAL number, when converted to decimal and 

then back to TEMP REAL, will remain the same. Such accuracy would require 

arithmetic beyond the TEMP REAL precision, which is too slow. The loss of 

accuracy is at worst less than 3 bits. 

The IEEE Standard (described in Chapter 6) specifies the accuracy for decimal con- 
versions of SHORT_REAL and LONG_REAL numbers. Since DCON87 uses 
TEMP REAL arithmetic to perform its conversions, it exceeds the IEEE specifica- 
tions. 

Following is the range of SHORT REAL and LONG_REAL numbers for which 

the IEEE standard specifies best possible conversion accuracy, for conversion both 
to and from the decimal format. The maximum number of decimal digits for best 
accuracy is given in the second column and the maximum magnitude for the decimal 
exponent is given in the third column. The exponents are given assuming the decimal 
point is to the right of the significand. For example, the largest number for totally 
accurate SHORT_REAL to decimal conversion is 999999999el3. The smallest 
positive number is le-13. 

Conversion Significand Exponent 

Decimal to SHORT_REAL 9 digits 13 

Decimal to LONG__REAL 19 digits 27 

SHORT_REAL to Decimal 9 digits 13 

LONG__REAL to Decimal 17 digits 27 

Following is a wider range of numbers for which the best possible conversion is not 
required. For numbers in the wider range, the IEEE standard allows an error of up 
to 0.47 times the unit in the last place. 
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Conversion Significand Exponent 

Decimal to SHORT_REAL 9 digits 99 

Decimal to LONG_REAL 19 digits 999 

SHORT_REAL to Decimal 9 digits 54 

LONG_REAL to Decimal 17 digits 341 



Error Reporting in DCON87.LIB 



DCON87 incorporates the 8087 error-reporting mechanism. Whenever an error is 
detected, an appropriate 8087 exception bit is set. Using 8087 instructions, you can 
set your choice of "unmasked" bits, and provide procedures which will be called 
whenever an unmasked exception is detected. This way, you do not have to test for 
errors every time you call a DCON procedure. 

Whenever an unmasked exception is detected by any DCON87 procedure, it places 
certain quantities into 8087 storage before calling your specified trap handler. 

• The opcode 00101 101 1 10 (16E hex) is placed into the 1 1-bit OPCODE register 
of the 8087. This is the opcode for the "FLDCW" instruction, which does not 
itself change the OPCODE register. You may therefore interpret the presence of 
the "FLDCW" opcode, together with an instruction address which is not all '1' 
bits, as meaning that the exception was discovered by DCON. (An all '1' bits 
instruction address signifies a CEL87 function. See Chapter 4.) 

• The 20-bit absolute address of the DCON87 procedure which caused the 
exception is placed into the instruction address register of the 8087. 

• The 20-bit absolute address of the input number is placed into the operand 
address register of the 8087. For mqcDEC__BIN and mqcDECLOW_BIN, this 
is a pointer into a DCON87 internal buffer on the 8086 stack, containing a 

string of decimal digits stripped for input to DECLOW BIN. For all the other 

procedures, this is the pointer to the original input binary number. 

An Alternate Method of DCON87 Error Detection 

If you prefer to perform error testing yourself, an alternate method is possible. All 
DCON87 procedures return the 8087 exception byte in the AL register. 

You can access the AL register in PL/M-86 programs by declaring the DCON87 
procedures to be of type BYTE. For example, you could change the declaration of 
mqcBIN_DECLOW to: 

mqcBIN_DECLOW: PROCEDURE (BL0CK_PTR) BYTE EXTERNAL; 

DECLARE BL0CK_PTR POINTER; 
END mqcBIN_DECLOW; 

The call to mqcBIN_DECLOW would then take the form: 

STATUS = mqcBIN_DECLOW(aBIN_DECLOW_BLOCK) ; 

After executing the above statement, you can test exception conditions by examining 
STATUS. 

Format of Input Decimal Numbers to mqcDEC_BIN 

The decimal number which is input to DEC BIN takes the form of a string of con- 
secutive ASCII bytes of memory, whose pointer is passed in the parameter block. 
The string may contain blanks after, but not before or within, the number. 
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The number may begin with a plus sign or a minus sign. Lack of a sign is equivalent 
to a plus sign. 

After the optional sign is the significand. The significand is a sequence of one or 
more decimal digits, with an optional decimal point. The decimal point may appear 
before, within, or after the decimal digits. The significand must be present, with at 
least one decimal digit. The upper limit to the number of digits is determined by the 
limit of 255 bytes to the entire string; however, only the 20 most significant non-zero 
digits of a long significand will affect the output answer. 

Following the significand is an optional exponent. This consists of one of the letters 
"E", "D", or "T", followed by an optional plus or minus sign, followed by a 
string of one or more decimal digits. The signed number indicates the power of ten 
which must be multiplied by the significand to yield the desired input number. The 
letters "E", "D", and "T" all have the same meaning (to identify the exponent). 
The letter can be in upper case or lower case. 

Examples of Input Decimal Strings 

The following strings all represent the number 100: 

100 

+100 

100. 

100.00 

1 00 .0000000000000000000060992 1 6948 

1E2 

1T2 

ld2 

1.E2 

.1E3 

100000d-3 

Here are some valid strings for other numbers: 

3.14159265 

-17 

-34.25T-009 

5.6E7 

.001 

-.4 

-0 

The following strings are invalid: 

1 E2 ; embedded space not allowed 

E2 ; missing significand 

,E2 ; significand must contain at least one decimal digit 
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mqcBIN_DECLOW (Convert binary number to decimal 
string, low level interface) 



Parameter: BIN__DECLOW_BLOCK__PTR, a double-word which points to a 
13-byte block of memory. The block consists of six fields containing the following 
variables: 

BIN PTR is a double-word which points to the input binary number. The amount 

of memory space occupied by the number is determined by BIN TYPE. 

BIN TYPE is a one-byte input code which indicates the type (SHORT REAL, 

LONG_REAL, or TEMP_REAL) of the number pointed to by BIN_PTR. The 
codes for the types are 0, 2, and 3, respectively. The caller must insure that 

BIN TYPE is set to the correct value — any other value will give an undefined 

result. 

DEC LENGTH is a one-byte input which indicates the number of decimal digits to 

be output at DEC PTR. All one-byte values are legal — if zero is given, only the 

DEC_EXPONENT and SIGN are returned, with no string at DEC_PTR. A zero 

DEC LENGTH combined with unusual inputs will yield an Invalid Result error, as 

described below. 

DEC PTR is a double-word which points to the memory location which will hold 

the output string of ASCII decimal digits. Other possible output values are shown in 
the table below DEC_SIGN. 

DEC^EXPONENT is a word variable, output by mqcBIN_DECLOW. The 

decimal point is assumed to be to the right of the digits output at DEC PTR. 

DEC EXPONENT is then the power of ten which, when multiplied by the number 

at DEC PTR, will yield the correct number. For NaN's and infinities, 

DEC_EXPONENT is 32767. For and -0, DEC_EXPONENT is 0. 

DEC_SIGN is a byte variable, output by mqcBIN_DECLOW. It is set to ASCII 
" + " (2B hex) if the input number is positive, and ASCII "-" (2D hex) if the 
number is negative. Other possible output values are shown below. 

Decimal outputs for unusual input numbers: 

Input number DEC_SIGN Decimal digits buffer 

NaN "." "." followed by blanks 

+ infinity "+" "+" followed by blanks 

-infinity "-" "-" followed by blanks 

+0 "0" blanks 

-0 "-" "0" followed by blanks 



Function: mqcBIN DECLOW converts the floating point number of type 

BIN TYPE, located at BIN PTR, into a decimal representation stored in 

DEC_PTR, DEC_EXPONENT, and SIGN. The caller will use these quantities to 
construct the final output format of the number. 

mqcBIN_DECLOW provides only 18 decimal digits. If DEC_LENGTH is greater 
than 18, the DEC_PTR buffer will be filled with 18 decimal digits, followed by 
ASCII underscore (5F hex) characters. The decimal point will be assumed to be to 
the right of the rightmost underscore. Thus, the underscore can be interpreted as 
"unknown decimal digit". 
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Errors: There are three possible errors for mqcBIN DECLOW: 

1. Invalid input. This error occurs only when the input DEC LENGTH is given 

as 0, and the input number is one of: -0, any infinity, or NaN. In these cases, the 
empty decimal digits buffer which the caller requested would have contained 

information other than decimal digits. The outputs DEC SIGN and 

DEC_EXPONENT are set appropriately. If the 8087 "I" exception is 
unmasked, the trap handler is called with the input number on the 8087 stack. 

2. Denormalized input. This error occurs only when BIN TYPE is 

TEMP REAL, and the leading significand bit is zero. DEC SIGN and 

DEC EXPONENT are given the correct output values. The output decimal 

digit string is given leading zeroes to indicate the extent of denormalization. If 
the 8087 "D" exception is unmasked, the trap handler is called with the input 
number on the 8087 stack. 

3. Inexact result. The correct answer is output. If the 8087 "P" exception is 
unmasked, the trap handler is called with the input number on the 8087 stack. 



PL/M-86 usage example: 

mqcBIN_DECL0W: PROCEDURE (BL0CK_PTR) EXTERNAL; 

DECLARE BL0CK_PTR POINTER; 
END mqcBIN_DECL0W; 

DECLARE REALJ/AR3 REAL; 
DECLARE DIG ITS_BUFFER (6) BYTE; 
DECLARE FINAL_BUFFER(15) BYTE; 

/* Assume that REAL_VAR3 contains the hex value 
C0E9A36D. */ 

BIN_DECL0W_BL0CK.BIN_PTR = 3REAL_VAR3; 
BIN_DECL0W_BL0CK.BIN_TYPE = SH0RT_REAL; 
BIN_DECL0W_BL0CK.DEC_LENGTH = 6; 
BIN_DECL0W_BL0CK.DEC_PTR = QDI G ITS_BUF FER ; 
CALL mqcBIN_DECLOW(aBIN_DECLOW_BLOCK); 

/* The ''building blocks' 1 returned by mqcBIN_DECL0W are 
assembled into a six-digit scientific notation 
format */ 

FINAL_BUFFER(0) = BIN_DEC L0W_B LOCK . DEC_S I GN ; 
FINAL_BUFFER(1) = DIG I T S_B U F F E R ( ) ; 
FINAL_BUFFER(2) = ' . ' ; 

CALL M0VB( aDIGITS_BUFFER(1) f 3F INAL_BUF FER (3) , 5 ); 
FINAL_BUFFER(8) = 'E'; 

/* Assume DEC I MA L_W0R D_0UT translates the first INTEGER 
parameter into a string of decimal digits placed at the 
second POINTER parameter */ 

CALL DECIMAL_W0RD_0UT (B IN_DEC L0W_B LOCK . DEC_EXP0NENT + 5, 
3FINAL_BUFFER(9)); 

/* FINAL_BUFFER now contains the string "-7.30120E0" */ 



BIN_DECLOW 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• n 3 i r s * 

EXTRN mqcBINJJECLOW: FAR 



REALJ/AR3 

DIG ITS_BUFFER_3 

FINAL_BUFFER 
BLOCK 3 



DD 0C0E9A36DR 

DB 6 DUP (?) 

DB 15 DUP (?) 

BIN_DECLOW_BLOCK < REAL_VAR3, 

SH0RT_REAL,?,DIGITS_BUFFER_3,?,? > 



MOV BL0CK_3.DEC_LENGTH, 6 

MOV DX, SEG(BL0CK_3) 

PUSH DX 

MOV DX, OFFSET (BL0CK_3) 

PUSH DX 

CALL mqcBIN_DECLOW 



plug in the correct 

Length 

upper part of param 

address 

goes onto stack 

Lower part of param 

address 

goes onto stack 



Now DIG ITS_BUFFER contains the string "730120", DEC_SIGN 
equals "-", and DEC_EXPONENT equals -5. This 
information can be used to construct a number in any 
format desi red . 
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mqcDEC_BIN (Convert decimal string to binary number) 



Parameter: DEC_BIN_BLOCK_PTR, a double-word which points to a ten- 
byte block of memory. The block of memory contains four fields which hold the 
following variables: 

BIN PTR is a double-word which points to the output binary number. The amount 

of memory space occupied by the number is determined by BIN TYPE. 

BIN TYPE is a one-byte input code which indicates the type (SHORT REAL, 

LONG_REAL, or TEMP_REAL) of the number pointed to by BIN„_PTR. The 
codes for the types are 0, 2, and 3, respectively. The caller must insure that 

BIN TYPE is set to one of the indicated legal values — any other value will give an 

undefined result. 

DEC LENGTH is a one-byte input which indicates the length of the memory field 

pointed to by DEC PTR. The length must be between 1 and 255; will give an 

undefined result. 

DEC PTR is a double-word which points to the input string of ASCII decimal 

digits. The format of the input string is given earlier in this chapter. The caller 
should already have verified correct input syntax before calling mqcDEC BIN. 



Function: mqcDEC BIN converts the string pointed to by DEC PTR into a 

binary floating-point number of type BIN TYPE, and leaves the number at the 

memory location pointed to by BIN PTR. 



Errors: There are six possible errors: 

1. Illegal input. If the input string is not of the specified format, the output is 
undefined. No warning is issued. 

2. Overflow beyond the TEMP_REAL format. If the 8087 "O" exception is 
masked, a correctly signed infinity is stored in the output buffer. If it is not 
masked, the 8087 trap handler is called with INDEFINITE on the 8087 stack, 
and nothing is written to the output buffer. 

3. Overflow beyond the output format, but not the TEMP REAL format. If the 

8087 "O" exception is masked, the correctly signed infinity is stored in the out- 
put buffer. If the 8087 "O" exception is unmasked, the trap handler is called. In 
this case, the top of the 8087 stack contains the output number with the signifi- 
cand rounded to the output format's precision. 

4. Underflow of the TEMP_REAL format. If the 8087 "U" exception is masked, 
then the output buffer contains the result of gradual underflow. If "U" is 
unmasked, INDEFINITE is left on the 8087 stack, the output buffer is 
unchanged and the 8087 trap handler is called. 

5. Underflow of the output format, but not of the TEMP REAL format. The 

output buffer contains the result of gradual underflow. If the 8087 "U" excep- 
tion is unmasked, the trap handler is called. In this case, the top of the 8087 
stack contains the output number with the significand rounded to the output 
format's precision. 

6. An inexact result. The output buffer contains the result, rounded to the output 
format. If the 8087 "P" exception is unmasked, the trap handler is called with 

the TEMP REAL result on the 8087 stack. Since this condition is not usually 

considered an error, the "P" exception is usually masked and the "P" error bit 
is ignored. 
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PL/M-86 usage example: 

mqcDEC_BIN: PROCEDURE (BL0CK_PTR) EXTERNAL; 

DECLARE BL0CK_PTR POINTER; 
END mqcDEC_BIN; 



DECLARE REAL_VAR REAL; 
DECLARE INPUT_BUFFER (120) 
DECLARE ACTUAL ADDRESS; 



BYTE 



/* Assume that at this point a decimal string (e.g. 
"-3.4E2") has already been placed into INPUT_BUFFER, 
and a buffer length (6 for the above number) has been 
placed into ACTUAL */ 



DEC_BIN_BLOCK.BIN 
DEC_BIN_BLOCK.BIN 
DEC_BIN_BLOCK . DEC 
DEC BIN BLOCK. DEC 



PTR = 3REAL_VAR; 
TYPE = SHORT_REAL; 
LENGTH = ACTUAL; 
PTR = 3INPUT BUFFER; 



CALL mqcDEC_BIN ODEC_B IN_B LOCK) ; 

/* RE A L_VAR now equals the encoded SHORT_REAL value 
the sample string given, REAL_VAR = -340 */ 



For 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 
; pairs: EXTRN mqcDEC_BIN: FAR 

INPUT_BUFFER DB 100 DUP (?) 

ACTUAL DW 

REAL_VAR DQ 

BL0CK_1 DEC_BIN_BLOCK < RIEAL_VAR, L0NG_REAL, ?, 

& INPUT_BUFFER > 

ENTRY: 

Assume that at this point a decimal string (e.g. 
"-3.4E2") has already been placed into INPUT_BUFFER, 
and a buffer length (6 for the above number) has been 
placed into ACTUAL. 



MOV AX, ACTUAL 

MOV BLOCKJ .DEC_LENGTH, AL 

MOV DX, SEG(BLOCKJ) 

PUSH DX 

MOV DX, OFFSET (BLOCK_1 ) 

PUSH DX 

CALL mqcDEC_BIN 



get length of st ri ng 

place byte into param 

block 

top half of param 

block pt r 

goes onto stack 

bottom half of param 

block pt r 

goes onto stack 

convert string to 

bi nary 



REAL_VAR now equals the encoded LONG_REAL value. For 
the sample given, REAL_VAR = -340. 
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mqcDECLOW__BIN (Convert decimal string, low-level inter- 
face, to binary number) 



Parameter: DECLOW_BIN_BLOCK_PTR, a double-word which points to a 
13-byte block of memory. The block of memory contains six fields which hold the 
following variables: 

BIN PTR is a double-word which points to the output binary number. The amount 

of memory space occupied by the number is determined by BIN TYPE. 

BIN TYPE is a one-byte input code which indicates the type (SHORT REAL, 

LONG_REAL, or TEMP_REAL) of the number pointed to by BIN_PTR. The 
codes for the types are 0, 2, and 3, respectively. The caller must insure that 

BIN TYPE is set to one of the indicated legal values — any other value will give an 

undefined result. 

DEC LENGTH is a one-byte input which indicates the length of the memory field 

pointed to by DEC PTR. The length must be between 1 and 21, inclusive. Any 

other value gives an undefined result. 

DEC PTR is a double-word which points to the input string of ASCII decimal 

digits. The digits give the significand of the number to be converted. Leading zeroes, 
the sign, the decimal point, the exponent, and lagging zeroes have already been 
removed. 

If DEC LENGTH equals 21, then there were more than 20 digits in the original 

significand. In this case, the given significand consists of the first 20 digits only. The 
twenty-first digit is the least significant non-zero digit of the original significand. 

DEC EXPONENT is a word which gives the base 10 exponent of the number to be 

converted. It is assumed that the decimal point has been shifted to the immediate 

right of the last digit of the significand. DEC EXPONENT is then the power of ten 

which, when multiplied by the significand, yields the number to be converted. 
Negative powers are represented in two's-complement form. 

DEC SIGN is a byte which gives the sign of the number to be converted. The possi- 
ble values are ASCII "+" (2B hex) and ASCII "-" (2D hex). 



Function: mqcDECLOW BIN accepts a decimal number which has already been 

converted from a higher-level format to the lower-level format described under 
DEC_PTR, DEC_EXPONENT, and DEC_SIGN above. It converts the number 

into a binary floating-point number of type BIN TYPE, and leaves the number at 

the memory location pointed to by BIN PTR. 



Errors: All errors (invalid input, overflow, underflow, and inexact result) are the 
same as for mqcDEC BIN, and have the same consequences. 



PL/M-86 usage example: 

mqcDECL0W_BIN: PROCEDURE (BL0CK_PTR) EXTERNAL; 

DECLARE BL0CK_PTR POINTER; 
END mqcDECL0W_BIN; 

DECLARE REAL_VAR2 REAL; 

DECLARE N_DIGITS BYTE; 

DECLARE DI6ITS_BUFFER(21) BYTE; 
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/* Assume that at this point the string "371" has been 

placed into D I G I TS_BU F FER , and N_D I G I TS has been set to 
the string length 3 */ 



/* The following code assumes D I G I T_BU F F ER contains a 
decimal string representing an integer number of 
pennies. It places the number of dollars, in SHORT 
binary format, into RE A L_V AR2 . */ 

PTR = £DREAL_VAR2; 
TYPE = SHORT_REAL; 
LENGTH = N_D I G ITS ; 
PTR = 3DIG ITS_BUF FER; 

EXPONENT = -2; 
-SIGN = '+'; 
C A L L mq c D E C LO W_B I N ( 3 B I N_D E C LO W_B LO C K ) ; 

/* REAL_VAR2 now equals the encoded SHORT_REAL value 
3.71 */ 



REAL 



B] 


EN 


DECLOW 


BLOCK 


BIN 


B 


[N 


"declow 


BLOCK 


BIN 


B 


[N 


DECLOW 


BLOCK 


DEC 


B] 


EN 


DECLOW 


BLOCK 


DEC 


B 


[N 


DECLOW 


BLOCK 


DEC 


B 


[N 


DECLOW 


BLOCK 


DEC 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• n 3 i r s * 

EXTRN mqcDECLOW_BINI: FAR 



N_DIG ITS 

DIG ITS_BUF FER 

REAL_VAR2 
BLOCK 2 



DB 

DB 21 DUP (?) 

DQ ? 

DEC L0W_BIN_BL0CK < REAL_VAR2, 

L0NG_REAL,?,DIGITS_BUFFER,-2, '+•> 



assume that the string ■ '371 • ■ has been placed into 
D I G ITS_BUF FER r and N_D I G ITS has been set to the string 
length 3 . 



MOV AL,N_DIGITS 

MOV BL0CK_2.DEC_LENGTH, AL 

MOV DX, SEG(BL0CK_2) 

PUSH DX 

MOV DX, OF FSET (BLOCK_2) 

PUSH DX 

CALL mqcDECLOW_BIN 



load string length 

p lace i nto param 

block 

top half of param 

block pt r 

goes onto stack 

bottom half of 

param block pt r 

goes onto stack 

convert string to 

binary 



; The encoded L0NG_ REAL value 3.71 has now been placed 
; into REAL VAR2. 
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mqcLONG_TEMP (Convert LONG_REAL to TEMP_REAL) 



Parameters: LONG_REAL PTR points to the eight bytes of memory which 

hold the input LONG REAL binary floating point number. 



TEMP REAL PTR points to the ten bytes of memory into which the 

TEMP REAL answer is placed. 



Function: The input LONG_REAL number is extended to the TEMP_REAL 
number which represents the same value. Unusual inputs (NaNs, infinities, -0) 
translate to the same unusual outputs, with zeroes padding the end of the output 
significand. 



Errors: The only possible error is a denormalized input value. The corresponding 
denormalized output value is left in the output buffer. If the 8087 "D" exception is 
unmasked, the trap handler is called with the output on the 8087 stack. 



PL/M-86 usage example: 

mqcLONG_TEMP: PROCEDURE ( LONG_REAL_PTR , TEMP_REAL_PTR) 
EXTERNAL; 
DECLARE (LONG_REAL_PTR,TEMP_REAL_PTR) POINTER; 
END mqcL0NG_TEMP; 

DECLARE TEMP_REAL_ARRAY (8) STRUCTURE 

(S IGNI F ICAND (8) BYTE, EXP0_SI6N(2) INTEGER); 
DECLARE LONG_REAL_ARRAY (8) STRUCTURE (NUM(8) BYTE); 
DECLARE I BYTE; 

/* The following code expands an array of L0NG_REAL 
variables into an array of TEMP_REAL variables */ 

DO I = TO 7; 

CALL mqcL0NG_TEMP ( 3L0NG_REAL_ARRAY ( I ) , 

3TEMP_REAL_ARRAY(I) ); 
END; 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• n a i r s * 

EXTRN mqcLONG_TEMP: FAR 

TEMP_NUM DT ? 
LONG NUM DD ? 
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; The following code transfers the LON G_RE A L number at 
; LONG NUM to the TEMP REAL number at TEMP NUM 



MOV DX f SEG(LONG_NUM) 

PUSH DX 

MOV DX, OFFSET ( L0NG_NUM) 

PUSH DX 

MOV DX, SEG(TEMP_NUM) 

PUSH DX 

MOV DX, OFFSET(TEMP_NUM) 

PUSH DX 

CALL mqcLONG_TEMP 



top half of parameter 

pointer 1 

goes onto stack 

bottom half of 

parameter pointer 1 

goes onto stack 

top half of parameter 

pointer 2 

goes onto st ac k 

bottom half of 

parameter pointer 2 

goes onto stack 



3-15 



SHORT_TEMP 

mqcSHORT_TEMP (Convert SHORT_REAL to TEMP_REAL) 



Parameters: SHORT REAL PTR points to the four bytes of memory which 

hold the input SHORT REAL binary floating point number. 



TEMP REAL PTR points to the ten bytes of memory into which the 

TEMP REAL answer is placed. 



Function: The input SHORT_REAL number is extended to the TEMP_REAL 
number which represents the same value. Unusual inputs (NaNs, infinities, -0) 
translate to the same unusual outputs, with zeroes padding the end of the output 
significand. 



Errors: The only possible error is a denormalized input value. The corresponding 
denormalized output value is left in the output buffer. If the 8087 "D" exception is 
unmasked, the trap handler is called with the output on the 8087 stack. 



PL/M-86 usage example: 

mqcSHORT_TEMP: PROCEDURE (SHORT_REAL_PTR , TEMP_REAL_PTR) 
EXTERNAL; 
DECLARE (SHORT REAL_PTR , TEMP_REAL_PTR) POINTER; 
END mqcSHORT_TEMP; 

DECLARE TEMP_REAL_ARRAY (8) STRUCTURE (NUMdO) BYTE); 
DECLARE SHORT_REAL_ARRAY (8) REAL; 
DECLARE I BYTE; 

/* The following loop expands an array of SHORT_REAL 
variables into an array of TEMP_REAL variables. */ 

DO I = TO 7; 

CALL mqcSHORT_TEMP ( 3SH0RT_REAL_ARRAY ( I ) , 

aTEMP_REAL_ARRAY(I) ); 
END; 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• d a i r s * 

EXTRN mqcSHORT_TEMP: FAR 

SH0RT_X DD ? 
TEMP X DD ? 
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; The following code converts the SHORT_REAL number at 
; SH0RT_X into a TEMP_REAL value, and leaves it at TEMP_X 



; It is assumed that both numbers are in the DS segment 
PUSH DS 
MOV AX, OFFSET (SH0RT_X) 



PUSH AX 
PUSH DS 

MOV AX, OFFSET (TEMP_X) 

PUSH AX 

CALL mqcSHORTJTEMP 



top half of f i rst 

parameter 

bottom half of first 

parameter 

pushed onto stack 

top half of second 

parameter 

bottom half of second 

parameter 

pushed onto stack 

convers ion is now 

comp lete 
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mqcTEMP__LONG (Convert TEMIP_REAL to LONG^REAL) 

Parameters: TEMP REAL PTR points to the ten bytes of memory which hold 

the input TEMP REAL binary floating point number. 

LONG REAL PTR points to the eight bytes of memory into which the 

LONG REAL answer is placed. 



Function: The TEMP_REAL number is rounded to the nearest LONG REAL 

number. If there are two LONG REAL numbers equally near, the one with a zero 

least significant bit is chosen. 



Errors: There are five possible errors: 

1. Overflow. If the exponent of the TEMP REAL input exceeds the output 

capacity, the correctly signed infinity is placed in the output buffer. If the 8087 
"O" exception is unmasked, the trap handler is called. In this case, the number 
left on the 8087 stack is the input number with the exponent unchanged and the 
significand rounded to the output's precision. 

2. Underflow. The output buffer is filled with the result of gradual underflow, 
rounded to fit the output format. If the 8087 "U" exception is unmasked, the 
trap handler is called. In this case, the 8087 stack is left with the input number's 
significand rounded to the output precision; with the too-small exponent 
unchanged. 

3. Inexact. This occurs when the input significand has non-zero bits outside the 
field of the output format. The correctly rounded answer is left in the output 
buffer. If the 8087 '?" exception is unmasked, the trap handler is called with 
the input number on the 8087 stack. 

4. Unnormalized invalid input. This occurs when the input is unnormalized, there 
is no underflow and rounding to the output precision does not normalize the 
answer. INDEFINITE is left in the output buffer. If the 8087 "I" exception is 
unmasked, the trap handler is called with the input number on the 8087 stack. 

5. Invalid nonzero NaN truncation. The only NaN inputs faithfully represented in 
the smaller format are those for which the significand bits being truncated are 
all zero. For other NaNs, invalid operation is signalled. The output buffer is 
filled with the truncated result only if the truncated significand is non-zero. If 
the significand would be zero, it is made non-zero by substituting the highest 
nonzero byte of the original significand into bits 5 through 12 of the new 

LONG REAL significand. If the 8087 "I" exception is unmasked, the trap 

handler is called with the input number on the 8087 stack. 



PL/M-86 usage example: 

mqcTEMP_LONG: PROCEDURE (TEMP_REAL_PTR, LONG_REAL_PTR) 
EXTERNAL; 
DECLARE (TEMP_REAL_PTR, LON G_R E A L_PT R ) POINTER; 
END mqcTEMP_LONG; 

DECLARE TEMP_REAL_ARRAY (8) STRUCTURE (NUMCIO) BYTE); 
DECLARE L0NG_REAL_ARRAY (8) STRUCTURE (NUM(8) BYTE); 
DECLARE I BYTE; 
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/ * The following code transfers an array of 
numbers to an array of L0NG_REAL numbers 

DO I = TO 7; 

CALL mqcTEMP_LONG(aTEMP_REAL_ARRAY(I), 

3L0NG_REAL ARRAY(D) ; 
END; 



8 TEMP 
*/ 



REAL 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

■ n a i r s * 

EXTRN mqcTEMP_LONG: FAR 

TEMP_REAL_ARRAY DT 8 DUP (?) 
LONG_REAL_ARRAY DQ 8 DUP (?) 

The following code transfers the array of 8 TEMP_REAL 
numbers in TEMP_REAL_ARRAY to the 8 L0NG_REAL numbers in 
L0NG_REAL_ARRAY. 

It is assumed that both arrays are in the current DS 
segment . 



PUSH BP 

MOV BP, OFFSET ( LON G_R E A L_A R R A Y ) 

MOV DX, OFFSET(TEMP_REAL_ARRAY) 



MOV CX, 8 

L00P8: 

PUSH CX 

PUSH DX 

PUSH DS 

PUSH DX 

PUSH DS 

PUSH BP 



CALL mqcTEMP_LONG 

POP DX 

ADD DX, 10 

ADD BP, 8 

POP CX 

DEC CX 

JNZ L00P8 

POP BP 



BP shou Id a Iways 
be preserved 

DX will not be 
preserved 
initialize loop 
counter 



preserve 
preserve 
pointer 
top ha If 
paramete 
bottom h 
TEMP par 
top ha If 
paramete 
bottom h 
LONG par 
perform 
restore 
pointer 
i nc remen 
TEMP_REA 
i nc remen 
LONG_REA 
restore 
count do 
loop bac 
conversi 
restore 
original 



counter 
TEMP_REAL 

of TEMP 
r 

alf of 
ameter 

of LONG 
r 

alf of 
ameter 
conversion 
TEMP_REAL 

t to next 
L 

t to next 
L 

counter 
wn 

k if more 
ons 
the 
BP 
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mqcTEMP_SHORT (Convert TEMP_REAL to SHORT_REAL) 

Parameters: TEMP REAL PTR points to the ten bytes of memory which hold 

the input TEMP REAL binary floating point number. 

SHORT REAL PTR points to the four bytes of memory into which the 

SHORT REAL answer is placed. 



Function: The TEMP REAL number is rounded to the nearest SHORT REAL 

number. If there are two SHORT REAL numbers equally near, the one with a zero 

least significant bit is chosen. 



Errors: The errors (overflow, underflow, inexact, invalid unnormalized, and 

invalid nonzero NaN truncation) are the same as for mqcTEMP LONG, with the 

same actions taken, except for the subcase of error 5 in which the highest nonzero 
byte of the original significand must be copied to the result significand to preserve its 
character as a NaN. In this case, that byte occupies bits through 7 of the new 
SHORT__REAL significand. 



PL/M-86 usage example: 

mqcTEMP_SHORT: PROCEDURE (TEMP_REAL_PTR , SHORT_REAL_PTR) 
EXTERNAL; 
DECLARE (TEMP_REAL_PTR,SHORT_REAL_PTR) POINTER; 
END mqcTEMP_SHORT; 

DECLARE TEMP_REAL_ARRAY (8) STRUCTURE 

(SIGNI F ICAND (8) BYTE, EXP0_SIGN(2) INTEGER); 
DECLARE SHORT_REAL_ARRAY (8) REAL; 
DECLARE I BYTE; 

/* The following loop transfers an array of 8 TEMP_REAL 
numbers (which could, for example, represent an image 
of the 8087 stack) to a more compact array of 
SHORT_REAL numbers. */ 

DO I = TO 7; 

CALL mqcTEMP_SHORT ( 3TEMP_R E A L_AR R A Y ( I) , 

3SH0RT_REAL_ARRAY(I) ); 
END; 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• d a i r s * 

EXTRN mqcTEMP_SHORT: FAR 

TEMP_VAR DT ? 
SHORT VAR DD ? 
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; The following code transfers the TEMP_REAL number at 
; TEMP VAR to the SHORT REAL number at SHORT VAR 



MOV DX, SE6(TEMP_VAR) 

PUSH DX 

MOV DX, OFFSET(TEMP_VAR) 

PUSH DX 

MOV DX, SEG(SHORT_VAR) 

PUSH DX 

MOV DX, OFFSET(SHORT_VAR) 

PUSH DX 

CALL mqcTEMP_SHORT 



top half of parameter 

poi nter 

goes onto stack 

bottom half of 

parameter pointer 

goes onto stack 

top half of parameter 

pointer 

goes onto stack 

bottom half of 

parameter pointer 

goes onto stack 
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8087 Support Library 



Linkage of DCON87.LIB to Your Program Modules 

The final action you must take to use DCON87.LIB in your programs is to insert the 
file name DCON87.LIB into a LINK86 command. 

DCON87 requires either the 8087 component or the 8087 emulator. If the compo- 
nent is present, you must also link in 8087. LIB. If the emulator is used, you must 
link in both E8087 and E8087.LIB. 

Following is the suggested order for the object modules in your LINK86 statement. 

Your object modules 

DCON87.LIB 

CEL87.LIB if you are using it 

EH87.LIB if you are using it 

8087. LIB if you are using the component, or 

E8087, E8087.LIB if you are using the emulator. 

As an example, if you are linking your PL/M-86 modules MYMODl.OBJ and 
MYMOD2.0BJ into a program using the 8087 emulator, issue the command 



> 
>> 



LINK86 :F1 : MYM0D1 .OBJ , : F1 : MYM0D2 . OBJ , &<cr> 

3 : FO : DC0N87 .LIB, :F0:E8087, : f : E8087 . LI B &<cr> 
TO : F1 : MYPROG . LNK<c r> 



If you have a single ASM-86 generated object module :Fl:MYPROG.OBJ to be 
executed in a system with an 8087 component, issue the command 



>> 



LINK86 :F1 :MYPR0G.0BJ, : FO : D C |l 8 7 . LI B , : FO : 8087 .LIB &<cr> 
TO :F1 :MYPR0G.LNK<cr> 
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CHAPTER 4 

THE COMMON ELEMENTARY 

FUNCTION LIBRARY 



Overview 

This chapter describes the functions of CEL87.LIB, a library that contains 
commonly-used functions of floating point arithmetic. 

The 8087 chip gives direct support for the most time-consuming part of the com- 
putation of every CEL87 function; but the forms of the functions and the ranges of 
the inputs are limited. CEL87 provides the software support necessary to let the 8087 
provide a complete package of elementary functions, giving valid results for all 
appropriate inputs. 

Following is a summary of CEL87 functions, grouped by functionality. 

Rounding and Truncation Functions: 

mqerlEX rounds a real number to the nearest integer; to the even integer if there is a 

tie. The answer returned is real. 
mqerIE2 is as mqerlEX, except the answer is a 16-bit integer. 
mqerIE4 is as mqerlEX, except the answer is a 32-bit integer. 
mqerlAX rounds a real number to the nearest integer; to the integer away from zero 

if there is a tie. The answer returned is real. 
mqerIA2 is as mqerlAX, except the answer is a 16-bit integer. 
mqerIA4 is as mqerlAX, except the answer is a 32-bit integer. 
mqerlCX chops the fractional part of a real input. The answer is real. 
mqerIC2 is as mqerlCX, except the answer is a 16-bit integer. 
mqerIC4 is as mqerlCX, except the answer is a 32-bit integer. 

Logarithmic and Exponential Functions: 

mqerLGD computes decimal (base 10) logarithms. 

mqerLGE computes natural (base e) lograithms. 

mqerEXP computes exponentials of base e. 

mqerY2X computes exponentials of any base. 

mqerY12 computes an input real to the AX integer power. 

mqerY14 is as mqerY12, except the input integer is DXAX. 

mqerYIS is as mqerY12, except the input integer is on the 8086 stack. 

Trigonometric and Hyperbolic Functions: 

mqerSIN, mqerCOS, mqerTAN compute sine, cosine, and tangent. 
mqerASN, mqerACS, mqerATN compute the corresponding inverse functions. 
mqerSNH, mqerCSH, mqerTNH compute the corresponding hyperbolic functions. 
mqerAT2 is a special version of the arc tangent function. It handles all the 
quadrants for rectangular-to-polar conversion in two dimensions. 

Other Functions: 

mqerDIM is FORTRAN'S positive difference function. 

mqerMAX computes the maximum of two real inputs. 

mqerMIN computes the minimum of two real inputs. 

mqerSGH combines the sign of one input with the magnitude of the other input. 

mqerMOD computes a modulus, retaining the sign of the dividend. 

mqerRMD computes a modulus, giving the value closest to zero. 



4-1 



The Common Elementary Function Library 
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All CEL87 procedure names being with the four letters "mqer". We have added this 
prefix to reduce the possibility that a CEL87 name will conflict with another name in 
your program. To make the names more readable, we have put the "mqer" prefix in 
lower case letters throughout this chapter. 



Declaring CEL87 Procedures in Your 
ASM-86 Programs 



In each source program module that calls a CEL87 procedure, you must declare that 
procedure as external before it is used. 

Following are declarations you should make at the top of your ASM-86 program. 
You should include only those procedures you use. 

To make reference easier, we have duplicated the EXTRN statements in the 
examples for the individual functions. You should not duplicate the EXTRN 
statements in your program. 

All CEL87 procedures must be declared FAR. 

; All of the following EXTRN statements should appear 
; outside of all SEGMENT-ENDS pairs 



EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 
EXTRN 



mqer ACS 
mqer AS N 
mqer AT 2 
mqerATN 
mqer COS 
mqerCSH 
mqerDIM 
mqer EXP 
mqer IA2 
mqerIA4 
mqer I AX 
mqer IC2 
mqerIC4 
mqerlCX 
mqerIE2 
mqer IE4 
mqer I EX 
mqerLGD 
mqerLGE 
mqerMAX 
mqerMIN 
mqerMOD 
mqerRMD 
mqerSGN 
mqer SIN 
mqerSNH 
mqer TAN 
mq e r T N H 
mqerY2X 
mqerYI2 
mq e r Y 1 4 
mqerYIS 



FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
FAR 
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Declaring CEL87 Procedures in Your 
PL/M-86 Programs 

Following are the declarations you should make at the top of your PL/M-86 pro- 
grams. You should include only those procedures you use. 

To make reference easier, we have duplicated the procedure declarations in the 
examples for the individual functions. You should not duplicate the declarations in 
your programs. 

To use CEL87.LIB, you must compile your PL/M-86 modules under the MEDIUM 
or LARGE models of computation, since CEL87's procedures are FAR procedures. 

mqerACS: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerACS; 

mqerASN: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerASN; 

mqerAT2: PROCEDURE (Y f X) REAL EXTERNAL; 

DECLARE (Y,X) REAL; 
END mqerAT2; 

mqerATN: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerATN; 

mqerCOS: PROCEDURE (THETA) REAL EXTERNAL; 

DECLARE THETA REAL; 
END mqerCOS; 

mqerCSH: PROCEDURE (THETA) REAL EXTERNAL; 

DECLARE THETA REAL; 
END mqerCSH; 

mqerDIM: PROCEDURE (Y,X) REAL EXTERNAL; 

DECLARE (Y f X) REAL; 
END mqerDIM; 

mqerEXP: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerEXP; 

mqerIA2: PROCEDURE (X) INTEGER EXTERNAL; 

DECLARE X REAL; 
END mqerIA2; 

mqerlAX: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerlAX; 

mqerIC2: PROCEDURE (X) INTEGER EXTERNAL; 

DECLARE X REAL; 
END mqerIC2; 

mqerlCX: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerlCX; 
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mqerIE2: PROCEDURE (X) INTEGER EXTERNAL; 

DECLARE X REAL; 
END mqerIE2; 

mqerlEX: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerlEX; 

mqerLGD: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerLGD; 

mqerLGE: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerLGE; 

mqerMAX: PROCEDURE (Y f X) REAL EXTERNAL; 

DECLARE (Y,X) REAL; 
END mqerMAX; 

mqerMIN: PROCEDURE (Y,X) REAL EXTERNAL; 

DECLARE (Y,X) REAL; 
END mqerMIN; 

mqerMOD: PROCEDURE (Y f X) REAL EXTERNAL; 

DECLARE (Y,X) REAL; 
END mqerMOD; 

mqerRMD: PROCEDURE (Y,X) REAL EXTERNAL; 

DECLARE (Y,X) REAL; 
fnh mqerRMD; 

GN: PROCEDURE (Y,X) REAL EXTERNAL; 

LARE (Y,X) REAL; 

qerSGN; 

mqerSIN: PROCEDURE (THETA) REAL. EXTERNAL; 

DECLARE THETA REAL; 
END mqerSIN; 

mqerSNH: PROCEDURE (THETA) REAL EXTERNAL; 

DECLARE THETA REAL; 
END mqerSNH; 

mqerTAN: PROCEDURE (THETA) REAL EXTERNAL; 

DECLARE THETA REAL; 
END mqerTAN; 

mqerTNH: PROCEDURE (THETA) REAL EXTERNAL; 

DECLARE THETA REAL; 
END mqerTNH; 

mqerY2X: PROCEDURE (Y,X) REAL EXTERNAL; 

DECLARE (Y,X) REAL; 
END mqerY2X; 

mqerYIS: PROCEDURE (Y,I) REAL EXTERNAL; 

DECLARE Y REAL, I INTEGER; 
END mqerYIS; 



END mqerRMD; 

mqerSG 

DECLttrcc \i 
END mqerSGN; 
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How to Call CEL87 Functions in Your 
PL/M-86 Programs 

CEL87 functions operate with numbers on the 8087 stack; however, this fact is invis- 
ible to you when you are programming in PL/M-86. There are no PL/M statements 
which explicitly load and unload numbers from the 8087 stack. The code for this is 
implicitly generated whenever arithmetic involving REAL numbers is programmed. 

To invoke a CEL87 function after you have declared it, simply use the function in an 
assignment statement (or in an expression in any other context). All CEL87 func- 
tions, if declared as shown in the examples, will cause the PL/M compiler to load 
the input parameters when CEL87 is invoked, and to place the answer in the correct 
destination after the CEL87 procedure is complete. 

PL/M-86 does not support the formats LONG_REAL and TEMP_REAL. 

Numeric constants and variables used by PL/M are 32-bit SHORT REAL 

numbers, in 8086 memory, which are converted to TEMP REAL when placed on 

the 8087 stack. When PL/M stores the result of a CEL87 function that returns a 

value on the 8087 stack, the function is converted from TEMP REAL to 

SHORT REAL. The conversion takes place with an FSTP instruction generated 

by the compiler. 

It is possible for the FSTP instruction to generate an "O" overflow or a "U" 

underflow error. (Indeed, because TEMP REAL is such a large format, this is 

where "O" and "U" errors are most likely to occur.) These errors are beyond the 
control of the CEL87 function; hence the outlines given in this chapter of what 
values are placed into 8087 registers do not apply for these errors. 

Consult the 8086 Family User's Manual Numerics Supplement for the error 
specifications of the FSTP instruction. 

Stack Usage of CEL87 Procedures 

CEL87 requires 50 bytes of the 8086 stack for its internal storage. If, however, there 
are situations in which CEL87 is interrupted by a program which itself calls CEL87 
(for example, an error trap handler), then an additional 50 bytes of 8086 stack are 
needed for each recursion level possible. CEL87 declares a 50-byte STACK segment 
within its program modules, hence the first 50 bytes will be automatically allocated 
byLINK86andLOC86. 

CEL87 requires 4 numbers on the 8087 stack. Input parameters are counted in the 
four numbers. This means that if a CEL87 function has input parameters at the top 
two stack positions ST(0) and ST(1), the parameters will be destroyed by the func- 
tion. Also, the numbers ST(6) and ST(7) must be empty for the function to execute 
correctly. If the 8087 stack overflows during a CEL87 function, the results are 
undefined. 

All CEL87 functions are reentrant in that they use no fixed 8086 memory locations 
to store internal variables. However, since the 8087 stack contains only eight 
elements, it is advisable for an interrupt routine to explicitly preserve the 8087 stack 
upon entry and to restore it upon exit. Any procedure that interrupts a CEL87 func- 
tion must save the entire 8087 stack, since it is impossible to tell which 4 stack 
elements are being used by CEL87. 
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Register Usage of CEL87 Procedures 

CEL87 conforms to the register-saving conventions of all 8086 high level languages 
provided by Intel. 

According to these conventions, the 8086 registers that can be destroyed by CEL87 
functions are AX, BX, CX, DX, DI, SI, and ES. The 8086 registers that are 
unchanged by CEL87 functions are BP, SS, and DS. The SP register is changed only 
by mqerYIS, which accepts a parameter on the 8086 stack and returns with the 
parameter popped from the stack. 

All CEL87 functions save the 8087 Control Word upon entry, and restore it upon 
return. During computation, all CEL87 functions except mqerDIM use 
TEMP_REAL precision, and a rounding mode of CEL87's choosing. mqerDIM 
uses the precision and rounding in effect at the time mqerDIM is called. 

The 8087 exception flags (error bits) are treated by all CEL87 functions as "sticky". 
This means that the only change to error bits that CEL87 may make is from off to 
on. If an error bit is already on when CEL87 is called, it will stay on. This allows you 
to perform a sequence of 8087 calculations that involve CEL87 functions with some 
of the 8087 errors masked. At the end of the calculation, if the exception flag for a 
masked error is still zero, you can be sure that the error was never encountered. 

Error Reporting in CEL87.LIE5 

CEL87 incorporates the 8087 error-reporting mechanism. Whenever an error is 
detected, an appropriate 8087 exception bit is set. Using 8087 instructions, you can 
set your choice of "unmasked" bits, and provide procedures that will be called 
whenever an unmasked exception is detected. This way, you do not have to test for 
errors every time you invoke a CEL87 function. 

There are only 7 CEL87 functions which leave a meaningful value in the "P" preci- 
sion error bit. These are mqerIC2, mqerIC4, mqerlCX, mqerIE2, mqerIE4, 
mqerlEX, and mqerSGN. For all other CEL87 functions, the "P" exception should 
be masked, and the "P" exception bit returned should be ignored. 

Whenever an unmasked exception is detected by any CEL87 function, the 8087 calls 
the trap handler. The values on the 8087 stack are decribed under each individual 
CEL87 function. In addition, CEL87 sets the following values whenever the trap 
handler is called: 

• The return address from the CEL87 function, converted to 20 bits, is placed in 
the operand register of the 8087. 

• The 8087 control word is restored to the value it had when the CEL87 function 
was called. 

• The 8087 instruction address register is set to 0FFFFF hex. 

• The 8087 opcode register is set to the code described under each individual 
CEL87 function. The description gives the code in the form of three hex- 
adecimal digits. The first (most significant) digit gives the number of values 
remaining on the 8087 stack that relate to the CEL87 function that caused the 
trap. The bottom two digits identify the CEL87 function. 

The trap handler can identify an interrupt as generated by CEL87 by examining the 
8087 instruction address register to see that it is set to 0FFFFF hex. If it is, then the 
identity of the CEL87 function is determined by examining the 8087 opcode register. 
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mqerACS x = arc cosine(x) 



Input parameter: (x) is the top number on the 8087 stack. 



Function: mqerACS returns the number, in radians, whose cosine is equal to the 
input value. Results are in the range to n. Valid inputs are from -1 to 1, inclusive. 
All zeroes, pseudo-zeroes, and denormals return the value n/2. Also, unnormals less 
than 2" 63 return the value tt/2. 



Output: The answer replaces the input on the 8087 stack. 



Errors: Since cosines exist only in the range -1 to 1, inclusive, all inputs to 
mqerACS outside of this range are invalid. Thus an "I" error is given for infinities, 
NaNs, values less than -1, and values greater than 1. An "I" error is given also for 
unnormals not less than 2" 63 . 

If "I" is unmasked, the trap handler is called with the input still on the stack, and 
the 8087 opcode register set to 175 hex. If "I" is masked, the answer is the input for 
NaNs; the answer is the value INDEFINITE for other invalid inputs. 



PL/M-86 usage example: 

mqerACS: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerACS; 

DECLARE HYPOTENUSE REAL; /* Longest side of a right 

t ri ang le */ 
DECLARE ADJACENT_SIDE REAL; /* The other side, next to 

angle computed */ 
DECLARE (THETA_RADIANS, TH ET A_D E.G R E E S ) REAL; /* Two forms 

of the 
answer */ 
DECLARE PI LITERALLY '3.14159265358979'; 

HYPOTENUSE = 10.; ADJ ACENT_S IDE = 5.; /* Test values */ 

/* The following lines calculate the value of an angle of 
a right triangle, given the length of the hypotenuse 
and the adjacent: side. mqerACS returns the value in 
radians; the value is then converted to degrees. * / 

THETA_RADIANS = mqerACS( ADJ ACENT_SI DE / HYPOTENUSE); 
THETA_DEGREES = (180. /PI) * T H ET A_R A D I AN S ; 

/* Now THETA_,DEGREES = 60 ~ it is a 30-60-90 degree 
t r i ang le */ 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

* d a i r s * 

EXTRN mqerACS: FAR 
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HYPOTENUSE 
ADJACENT SIDE 



THETA_RADIANS 
THETA_DEGREES 
RAD TO DEG 



DQ 10.0 ; Longest side of a right 

; triangle 
DQ 5.0 ; the other side, next to 
; angle computed 
; The above initial values are test 
; values 
DQ ? 
DQ ? 

DT 4004E52EE0D31E0FBDC3R ; the constant 

; 180/PI 



The following lines compute the angle of a right 
triangle, just as in the above PL/M example, except 
with LONG REAL variables. 



FLD ADJACENT_SIDE 
FDIV HYPOTENUSE 

CALL mqerACS 

FST THETA_RADIANS 

FLD RAD_T0_DEG 

FMUL 

FSTP THETA DEGREES 



(x) = ADJ ACENT_S IDE 

(x) = ADJACENT_SIDE/ 

HYPOTENUSE 

angle is now in (x) 

radians result is saved 

(x) converted to degrees 
degrees are saved, 8087 
; stack popped 



; Now THETA_DEGREES = 60 ~ It is a 30-60-90 
; t ri ang le. 
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mqerASN x = arc sine(x) 



Input parameter: (x) is the top number on the 8087 stack. 



Function: mqerASN returns the number, in radians, whose sine is equal to the 
input value. Results are in the range -n/2 to tt/2. Valid inputs are from -1 to 1, 
inclusive. All zeroes, pseudo-zeroes, and denormals return the input value. Also, 
unnormals less than 2~ 63 return the input value. 



Output: The answer replaces the input on the 8087 stack. 



Errors: Since sines exist only in the range -1 to 1, inclusive, all inputs to mqerASN 
outside of this range are invalid. Thus an "I" error is given for infinities, NaNs, 
values less than -1, and values greater than 1. An "I" error is given also for 
unnormals not less than 2" 63 . 

If "I" is unmasked, the trap handler is called with the input still on the stack, and 
the 8087 opcode register set to 174 hex. If "I" is masked, the answer is the input for 
NaNs; the answer is the value INDEFINITE for other invalid inputs. 



PL/M-86 usage example: 

mqerASN: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerASN; 

DECLARE HYPOTENUSE REAL; /* Longest side of a right 

t r i ang le */ 
DECLARE 0PP0S I T E_S I DE REAL; /* The other side, opposite 

the angle computed */ 
DECLARE (THETA_RADIANS f THETA_DEGREES) REAL; /* Two forms 

of the 
answer */ 
DECLARE PI LITERALLY '3.14159265358979'; 

HYPOTENUSE = 10.; 0PP0S ITE_S IDE = 5.; /* Test values */ 

/ * The following lines calculate the value of an angle of 
a right triangle, given the length of the hypotenuse 
and the opposite side. mqerASN returns the value in 
radians; it is then converted to degrees. */ 

THETA_RADIANS = mqerACSC 0PP0S ITE_S IDE / HYPOTENUSE); 
THETA_DEGREES = (180. /PI) * T H ET A_R AD I AN S ; 

/* Now THETA_DEGREES = 30 -- it is a 30-60-90 degree 
t r i ang le */ 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

; pairs: 

EXTRN mqerASN: FAR 
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HYPOTENUSE 
OPPOSITE SIDE 



THETA_DE6REES 
RAD TO DE6 



DQ 10.0 ; Longest side of a right 

; triangle 
DQ 5.0 ; the other side, next to 
; angle computed 
; The above initial values are test 
; values 
DQ ? 

DT 4004E52EE0D31E0FBDC3R ; the constant 

; 180/PI 



The following lines compute an angle of a right 
triangle, just as in the above PL/M example; except 
with LONG REAL variables. 



FLD 0PP0SITE_SIDE 
FDIV HYPOTENUSE 

CALL mqerASN 

FLD RAD_T0 DEG 

FMUL 

FSTP T.HETA DEGREES 



; (x) := OPPOS ITE_SIDE 

; (x) = 0PP0SITE_S IDE/ 

; HYPOTENUSE 

; angle is now in (x) 

; (x) converted to degrees 

; degrees are saved, 8087 

; stack popped 



; With the test inputs, THETA_DEGREES = 30; it is a 
; 30-60-90 triangle 
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mqerAT2 



x = arc tangent(y/x) 



Input parameters: (x) is the top number on the 8087 stack; (y) is the next number 
on the 8087 stack. 



Function: mqerAT2 performs one half of a rectangular-to-polar coordinate con- 
version. If the inputs (x) and (y) are interpreted as the rectangular coordinates of a 
point in a plane, mqerAT2 returns the angle, in radians, which the point is displaced 
from the positive X-axis. See figure 4-1. The angle ranges from n (points just below 
the negative X-axis) to it (points on or just above the negative X-axis). 

The formulas used to calculate the angle are as follows: 

If (x) > 0, return (arc tangent(y/x)). 
If (x) = 0, return mqerSGN (tt/2, y). 
If (x) < 0, return (arc tangent(y/x)) + mqerSGN (n, y). 

It is legal for one (but not both) of the inputs to be an infinity. The point is then 
thought of as "infinitely far out" along one of the axes, and the angle returned is the 
angle of that axis. For points near the X-axis, the angle retains the sign of (y). Thus, 
x = +INFINITY returns with the sign of (y); x = -INFINITY returns n with the 
sign of (y); y = +INFINITY returns +tt/2; y = -INFINITY returns -tt/2. In all these 
cases, it is legal for the non-infinite input to be unnormal. Also note that the distinc- 
tion between +INFINITY and -INFINITY is made even when the 8087 is in projec- 
tive (unsigned infinity) mode. 

mqerAT2 accepts denormal inputs. Its action depends on whether the 8087 is in nor- 
malizing mode, as indicated by the setting of the "D" error masking bit. If the 8087 
is in normalizing mode ("D" is unmasked), then any denormals are assumed to be 
zero. If the 8087 is in warning mode ("D" is masked), the denormals are replaced 
with unnormals with the same numeric value. Note that even though mqerAT2 
checks the "D" masking bit, it does not set the "D" error bit; nor does it call the 
"D" trap handler. 



Y-AXIS 



(X,Y) 



mqerAT2(Y,X) 



X-AXIS 



Figure 4-1. Rectangular-to-Polar Coordinate Conversion 121725-1 



4-11 



AT2 



In some cases, mqerAT2 accepts an unnormal input with no error. This can only 
happen when (y) and (x) represent the coordinates of a point which is extremely close 
to one of the axes. A point is considered close enough to the X-axis if the absolute 
value of (y/x) is less than 2~ 63 . A point is considered close enough to the Y-axis if the 
absolute value of (x/y) is less than 2" 63 . If the point is near the positive X-axis, the 
value returned is the tiny ratio (y/x) with the sign of (y). If the point is near one of 
the other three axes, the value returned is the angle of that axis: n with the sign of (y) 
for the negative X-axis, n/2 for the positive Y-axis, and -n/2 for the negative 
Y-axis. Under no circumstances is it legal for both inputs to be unnormal. 



Output: The 8087 stack pops once, with the answer replacing the two inputs. 



Errors: The first thing mqerAT2 does is check for input NaN's. If either input is a 
NaN, there is an "I" error. If "I" is masked, the input NaN is returned (if both 
inputs are NaN's, the larger NaN is returned). 

An "I" error is given in a number of other cases. If both inputs are zero (i.e., the 
point on the coordinate plane is the origin), the angle cannot be defined; so an "I" 
error is given. If an unnormal input is given which does not fit any of the legal cases 
described above, an "1" error is given. For all these "I" errors, if "I" is masked, 
the value INDEFINITE is returned. 

If "I" is unmasked for any "I" error, the trap handler is called with the original 
inputs on the 8087 stack, and the 8087 opcode register set to 277 hex. 

If the magnitude of the result is too tiny to be represented in TEMP REAL format, 

a "U" underflow error results. If "U" is masked, mqerAT2 returns a gradual- 
underflow denormal if possible; if not possible. If "U" is unmasked, the trap 
handler is called with the 8087 opcode register set to 277 hex; but the inputs are not 
left on the 8087 stack. Instead, the correct answer is placed on the stack, with the 
too-small exponent presented in "wrapped" form. To obtain the true exponent, 
subtract decimal 24576 from the given exponent. 



PL/M-86 usage example: 

mqerAT2: PROCEDURE (Y,X) REAL EXTERNAL; 

DECLARE (Y,X) REAL; 
END mqerAT2; 

FSQRT87: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END FSQRT87; 

/* The above procedure is simply a PUBLIC ASM-86 procedure 
which consists of the FSQRT instruction, followed by a 
RET instruction. * / 

DECLARE REC_X, REC_Y REAL; 

DECLARE (P0LAR_THETA, P0LAR_R) REAL; 

REC_X = 3.; REC._Y = 3. * FSQRT87C3.) ; /* Test values */ 

/* The following statements convert the rectangular 

coordinates REC_X and RE C_Y into the polar coordinates 
POLAR R and POLAR THETA. */ 



4-12 



AT2 



POLAR_R = FSQRT87 (REC_X * REC_X + REC_Y * REC_Y); 
POLAR_THETA == mqe rAT2 ( REC_Y , REC__X) ; 

/* Now POLAR R = 6 and POLAR THETA = PI/3. */ 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

; pairs: 

EXTRN mqerAT2: FAR 



POLAR_THETA 

P0LAR_R 

REC_X 

REC Y 



DQ ? 

DQ ? 

DQ 3.0 

DQ 3.0 



FLD REC_Y 

FLD REC_X 

CALL mqe r AT 2 
FSTP POLAR_THETA 

FLD REC_Y 

FMUL ST, ST 

FLD REC_X 

FMUL ST, ST 

FADDP ST(1) f ST 

FSQRT 

FSTP POLAR R 



rectangular X coordinate 

initialized to test value 

rectangular Y coordinate 

initialized to test value 

Y-coordinate parameter goes on 

stack 

X-coordinate parameter goes on 

stack 

angle is computed 

and stored 

Y onto stack 

Y is squared 
X onto stack 
X is squared 
sum of squares 
radius is on stack 

radius is stored and stack is 
popped 



; Now POLAR_THETA = PI/4 and 
; P0LAR_R = 3 * Squareroot (2) 
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mqerATN x = arc tangent(x) 

Input parameter: (x) is the top number on the 8087 stack. 



Function: mqerATN returns the number, in radians, whose tangent is equal to the 
input value. Numbers are chosen between -: n/2 and n/2. All zeroes, pseodo-zeroes, 
and denormals return the input value. Also, unnormals less than 2" 63 return the in- 
put value. 

Infinite inputs are valid, whether the 8087 is in affine (signed infinity) of projective 
(unsigned infinity) mode. The input -INFINITY returns the value -n/2; the input 
+INFINITY returns the value +n/2. 



Output: The answer replaces the input on the 8087 stack. 



Errors: An input NaN gives an "I" error. Also, an unnormal not less than 2" 63 
gives an "I" error. If "I" is masked, the answer returned is the input for NaN's; the 
answer is the value INDEFINITE for illegal unnormals. If "I" is unmasked, the 
trap handler is called, with the input still on the 8087 stack and the value 176 hex in 
the 8087 opcode register. 



PL/M-86 usage example: 

mqerATN: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerATN; 

DECLARE 0PP0SITE_SIDE REAL; 

DECLARE ADJACENT_SIDE REAL; /* Shorter sides of a right 

t ri ang le */ 
DECLARE THETA_RADIANS REAL; 

0PP0SITE_SIDE = 1.; ADJ ACENT_S IDE = 1.; /* Test values */ 

/* The following line computes the angle of a right 
triangle, given the lengths of the two shorter 
si des . */ 

THETA_RADIANS = mqe rATN (0PP0S ITE_S I DE/ADJ ACENT_S I DE) ; 

/* THET A_RAD I ANS now equals the value of the angle: 
PI/4. */ 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 
EXTRN mqerATN: FAR 
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OPPOSITE_S IDE DQ 5.0 

ADJACENT SIDE! DQ 5.0 



short side, opposite angle 

computed 

short side, next to angle 

computed 

The above initial values are test 
va lues 



THETA RADIANS DQ ? 



The following lines c omp ute an angle of a right 
triangle, just as in the above PL/M example, except 
with LONG REAL variables. 



FLD OPPOS ITE_SIDE 

FDIV ADJACENT_SIDE 

CALL mqerATN 

FSTP THETA RADIANS 



(x) := 0PP0SITE_SIDE 

(x) = 0PP0SITE_SIDE/ 

ADJACENT_SIDE 

ang le is now in (x) 

radians are saved, 8087 

stack popped 



; With the test inputs, THETA_RADIANS = PI/4. 
; It is a 45-45-90 triangle. 
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mqerCOS x = cosine(x) 

Input parameter: (x) is the top number on the 8087 stack. 



Function: mqerCOS returns the trigonometric cosine of x, where x is an angle 
expressed in radians. All input zeroes, pseudo-zeroes, and denormals return the 
value 1 . Also, unnormals whose value is less than 2~ 63 return the value 1 . 



Output: The answer replaces the input on the 8087 stack. 



Errors: An "I" error is given for input infinities and NaN's. An "I" error is also 
given for unnormals which do not represent values less than 2" 63 . 

If "I" is unmasked, the trap handler is called with the input still on the stack, and 
the 8087 opcode register set to 172 hex. If "I" is masked, the answer is the input for 
NaNs; the answer is the value INDEFINITE for other invalid inputs. 



PL/M-86 usage example: 

mqerCOS: PROCEDURE (THETA) REAL EXTERNAL; 

DECLARE THETA REAL; 
END mqerCOS; 

DECLARE (P0LAR_R, P0LAR_THETA) REAL; 
DECLARE REC_X REAL; 

DECLARE PI LITERALLY '3.141 59265358979 ' ; 
DECLARE DEG_T0;RAD LITERALLY 'PI/180.'; 

P0LAR_R = 2.; P0LAR_THETA = 30.; /* Test values */ 

/* The following line computes the X-coordinate of a 

polar-to-rectangular conversion. The input angle is in 
degrees, so it must be converted to radians. */ 

REC_X = P0LAR_R * mqerCOS (P0LAR_THETA * DEG_T0_RAD); 

/* Now in the test case, REC_X = the square root of 3. */ 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• d a i r s * 

EXTRN mqerCOS: FAR 
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POLAR_THETA 
POLAR R 



REC_X 

DE6 TO RAD 



DQ 30.0 
DQ 2.0 

; the above initializations are test 

; values. 
DQ ? 

DT 3FF98EFA351294E9C8AER ; the constant 

; PI/180. 



The following lines compute theX-coordinate of a 
polar-to-rectangular conversion, as in the PL/M 
example above; except that the variables are 
LONG REAL. 



FLD POLAR_JHETA 

FLD DEG_T0_RAD 

FMUL 

CALL mqerCOS 

FMUL POLAR_R 

FSTP REC X 



degrees angle onto 8087 stack 

converted to radians 

cosi ne i s taken 

answer scaled to correct radius 

X-coordinate stored and stack is 

popped 



; With the test case, RE C_X is now the square root of 
: 3. 
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mqerCSH x = hyperbolic cosine(x) 

Input parameter: (x) is the top number on the 8087 stack. 

Function: mqerCSH returns the hyperbolic cosine of x, where x is an angle 
expressed in radians. All input zeroes, pseudo-zeroes, and denormals return the 
value 1 . Also, unnormals whose value is less than 2~ 63 return the value 1 . 

Either infinity returns the value +INFINITY with no error. 
Output: The answer replaces the input on the 8087 stack. 

Errors: An "I" error is given for input NaN's. An "I" error is also given for 
unnormals which do not represent values less than 2~ 63 . 

If "I" is unmasked, the trap handler is called with the input still on the stack. If "I" 
is masked, the answer is the input for NaNs; the answer is the value INDEFINITE 
for invalid unnormals. 

mqerCSH will give an "O" overflow error if the input is greater than about 11355, 
or less than about -11355. When "O" is masked, the value +INFINITY is returned. 
When "O" is unmasked, the trap handler is called with the input still on the 8087 
stack. 

When the trap handler is called, mqerSNH first sets the 8087 opcode register to 16F 
hex. 

PL/M-86 usage example: 

mqerCSH: PROCEDURE (THETA) REAL EXTERNAL; 

DECLARE THETA REAL; 
END mqerCSH; 

DECLARE (INPUT_VALUE, 0UTPUT_VALUE) REAL; 

INPUTJ/ALUE = 1.3; /* Test value */ 

OUTPUT_VALUE = mqerCSH ( INPUTJ/ALUE) ; 

/* For the test case, 0UTPUT_VALUE now is approximately 
1.97091 */ 

ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

* d a i r s * 

EXTRN mqerCSH: FAR 

INPUTJ/ALUE DQ 1.56 ; test value 
OUTPUT VALUE DQ ? 



FLD INPUTJ/ALUE 
CALL mqerCSH 
FSTP OUTPUT VALUE 



input goes onto 8087 stack 
hyperbolic cosine is taken 
answer is stored in 8086 memory, 
8087 is popped 



; For the test value, 0UTPUTJ/ALUE now is about 2.48448 
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mqerDIM x = max(y-x, + 0) 



Input parameters: (x) is the top number on the 8087 stack; (y) is the next number 
on the 8087 stack. 



Function: mqerDIM first compares (y) and (x); if (x) is greater than (y), the answer 
+0 is returned. If (y) is greater than (x), the positive value (y - x) is returned. 

When the 8087 chip is in affine (signed infinity) mode, and when both inputs are not 
the same, then infinite inputs are allowed. The results are as expected from the 
above comparison and subtraction. That is, if (x) is +INFINITY or (y) is 
-INFINITY the answer is +0; if (x) is -INFINITY or (y) is +INFINITY, the answer 
is+INFINITY. 

mqerDIM uses the 8087 precision and rounding modes which were in effect at the 
time mqerDIM is called. Thus the results can vary, depending on the settings of 
these bits in the 8087 control word. 



Output: The 8087 stack pops once, with the answer replacing the two inputs. 



Errors: The first error which is detected is the "D" error, for unnormal or denor- 
mal inputs. If the 8087 is in warning mode ("D" is masked), the calculation con- 
tinues; but the "P" error bit remains set. If the 8087 is in normalizing mode ("D" is 
unmasked), the trap handler is called. The trap handler is called directly from the 
interior of mqerDIM; the 8087 opcode register contains the code for the instruction 
which found the denormal: either FCOM or FSUB. The (y) and (x) inputs are left on 
the 8087 stack. Most trap handlers will replace the denormal input(s) with normal- 
ized numbers, reexecute the FCOM or FSUB instruction, and continue through 
mqerDIM. The trap handler provided by EH87.LIB (described in Chapter 5) will do 
this, and replace denormals with zero. 

mqerDIM next checks for NaN inputs. If either input is a NaN, an "I" error is 
given. If "I" is masked, the answer returned is the input NaN (the larger NaN if 
both inputs are NaN's). 

An "I" error is also given if the 8087 chip is in projective (unsigned infinity) mode, 
and either input is infinity; or if the 8087 chip is in affine (signed infinity) mode, and 
both inputs are the same infinity. In these cases a masked "I" yields the value 
INDEFINITE for an answer. 

For all "I" errors, when "I" is unmasked, the trap handler is called with the inputs 
still on the 8087 stack, and the 8087 opcode register set to 265 hex. 

In the case (y) > (x), when the answer is calculated by the subtraction (y - x), the 
subtraction can cause either an "O" overflow error or a "U" underflow error. If 
the error is masked, the answer returned is +INFINITY for overflow. For 
underflow, a gradual-underflow denormal is returned if possible; is returned if not 
possible. If the error is unmasked, the trap handler is called, with the value 165 hex 
in the 8087 opcode register. The inputs are not on the 8087 stack. Instead, the 
correct answer appears on the 8087 stack, with the out-of-range exponent given 
"wrapped around". For underflow, the correct exponent is the given exponent 
minus the decimal value 24576; for overflow, it is the exponent plus 24576. 
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PL/M-86 usage example: 

mqerDIM: PROCEDURE (Y,X) REAL EXTERNAL; 

DECLARE (Y f X) REAL; 
END mqerDIM; 

DECLARE (COSTS, RECEIPTS) REAL; 
DECLARE (PROFIT, LOSS) REAL; 

COSTS = 4700.00; RECEIPTS = 5300.00; /* Test values */ 

/* The following lines return the profit or loss, given 
the costs and receipts. The positive difference goes 
into PROFIT or LOSS as appropriate, and the other value 
i s set to zero . */ 

PROFIT = mqerDIM(RECEIPTS, COSTS); 
LOSS = mqerDIMUGSTS, RECEIPTS); 

/* For the test case, PROFIT is now 600.00, and LOSS 
is 0. */ 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• d a i r s " 

EXTRN mqerDIM: FAR 

COSTS DQ 800.00 
RECEIPTS DQ 650.00 

; the above initializations are test cases 
PROFIT DQ ? 
LOSS DQ ? 

; The following lines compute profit and loss, just as 
; in the PL/M example. 



FLD RECEIPTS 
FLD COSTS 
CALL mqerDIM 
FSTP PROFIT 
FLD COSTS 
FLD RECEIPTS 
CALL mqerDIM 
FSTP LOSS 



first parameter to mqerDIM 
second parameter to mqerDIM 
positive difference is computed 
answer is stored and stack is popped 
now perform function the other way 



other answer goes to LOSS, and stack 
is again popped 



; For this test case, LOSS = 150.00 and PROFIT is 
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mqerEXP x = e x 

Input parameters: x is the top number on the 8087 stack. 



Function: mqerEXP raises the constant e (2.718281828459+) to the power of the 
input number. This is a standard function because it is used to derive other func- 
tions, notably the hyperbolics. Also, on many computers it is slightly easier to com- 
pute than exponentials based on other constants. The exponential is valid for both 
positive and negative inputs. 

All input zeroes, pseudo-zeroes, denormals, and unnormals less than 2~ 63 yield an 
answer of 1 . 

If the 8087 is in affine (signed infinity) mode, infinite inputs are valid. +INFINITY 
returns itself as the answer; -INFINITY returns 0. 



Output: The answer replaces the input on the 8087 stack. 



Errors: All input NaN's, and input unnormals greater than 2" 63 , cause an "I" 
error. Also, if the 8087 is in projective (unsigned infinity) mode, both infinite inputs 
give an "I" error. If "I" is masked, the value returned is the input if a NaN; the 
value INDEFINITE otherwise. If "I" is unmasked, the trap handler is called with 
the input still on the 8087 stack. 

Since the largest number which can be stored in the 8087's temporary real format is 
about e 11357 , an input greater than about 11357 causes an "O" overflow error. 
When "O" is masked, the value +INFINITY is returned. Likewise, an input less 
than about -11355 causes the "U" underflow error. If "U" is masked, the result is 
a gradual-underflow denormal if possible, zero otherwise. For unmasked "O" or 
"U" errors, the trap handler is called, with the input still on the 8087 stack. 

When the trap handler is called, mqerEXP first places the value 16B hex into the 
8087 opcode register. 



PL/M-86 usage example: 

mqerEXP: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerEXP; 

DECLARE Y_VALUE REAL; 

DECLARE X_VALUE REAL; 

DECLARE NORM CONSTANT LITERALLY '0.3989422803' ; 

/* 1 / <SQRT(2*PI>) */ 

X_VALUE = -2.0; /* Test value */ 

/* The following line gives the value for the graph of the 
normal distribution function, for mean and variance 
1. By plotting Y_VALUE against various different input 
X- values, one obtains the familiar '' bell-shaped 
curve 1 '. * / 
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Y VALUE = NORM_CONSTANT 

* mqerEXP( - X_VALUE * X_VALUE 12.)', 

/* For the test input, Y_VALUE is now approximately 
0.5399100 */ 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• n a i r s * 

EXTRN mqe.rEXP: FAR 



MINUS_2 

X_VALUE 
Y_VALUE 
NORM CONSTANT 



DQ -2. 



constant used in calculation 

below 

initialization is a test value 



DQ 0.4 

DQ ? 

DT 3FFDCC42299EA1B28468R 



; constant 

; 1/sqrt(2*PI) 



The following lines compute the normal distribution 
just as in the PL/M example, except with L0NG_REAL 
numbers . 



FLD X_VALUE 
FMUL ST, ST 
FIDIV MINUS_2 
CALL mqerEXP 
FLD N0RM_C0NSTANT 
FMUL 



load input onto 8087 stack 
input is squared 
negate and divide by 2 
exponent i ate 



divide by the square root of 

2 * PI 

store the result and pop the 

stack 



FSTP Y_VALUE 

; For the test case Y_VALUE is now about 0.3399562 
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mqerlA2 AX =■ roundaway(x) 

Input parameters (x) is the top number on the 8087 stack. 



Function: If x is not normal it is first normalized. Then x is rounded to the nearest 
integer. If there are two equally near integers (i.e., the fractional part of x is .5), then 
the integer farthest from zero is selected. The answer is given in 16-bit two's comple- 
ment form. For example, 3.1 returns hex 0003; 10.5 returns hex 000B; -2.5 returns 
hex FFFD, which is -3 in two's complement form. 



Output: The input is popped from the 8087 stack, and the answer is left in the AX 
register. 



Errors: The only numbers that will fit the 16-bit destination are those between, but 
not including, -32768.5 and +32767.5. Numbers not in this range, including 
infinities and NaN's, will cause an "I" error. If "I" is masked, the "indefinite 
integer" value 8000 hex is returned. If "I" is unmasked, the trap handler is called 
with the input still on the 8087 stack, and the 8087 opcode register set to 17E hex. 



PL/M-86 usage example: 

mqerIA2: PROCEDURE (X) INTEGER EXTERNAL; 

DECLARE X REAL; 
END mqerIA2; 

DECLARE REAL VAR REAL; 
DECLARE INTEGER_VAR INTEGER; 

REAL_VAR = 4„5; /* Test value */ 

INTEGERJ/AR = mqe r I A2 (REAL_VAR) ; 

/* Now in the test case, INTEGER_VAR = 0005 hex. */ 

ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• d a i r s * 

EXTRN mqerIA?.: FAR 

INTEGER_VAR DW ? 

REAL_VAR DQ -8.5 ; Initialization is a test value 

FLD REAL_VAR ; load the parameter 
CALLmqerIA2 ; round to nearest integer 
MOV INTEGER_VAR,AX ; store the answer 

; For the test case, I NTEG ER_ VAR now equals -9, which 
; is FFF7 hex. 
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mqerlA4 DXAX = roundaway(x) 

Input parameter: (x) is the top number on the 8087 stack. 



Function: If x is not normal it is first normalized. Then x is rounded to the nearest 
integer. If there are two equally near integers (i.e., the fractional part of x is .5), then 
the integer farthest from zero is selected. The answer is given in 32-bit two's comple- 
ment form. For example, 3.1 returns hex 00000003; 10.5 returns hex 0000000B; -2.5 
returns hex FFFFFFFD, which is -3 in two's complement form. 



Output: The input is popped from the 8087 stack, and the answer is left in the DX 
and AX registers, with DH the highest byte and AL the lowest byte. 



Errors: The only numbers that will fit the 32-bit destination are those between, but 
not including, -2147483648.5 and +2147483647.5. Numbers not in this range, 
including infinities and NaN's, will cause an "I" error. If "I" is masked, the 
"indefinite integer" value 80000000 hex is returned. If "I" is unmasked, the trap 
handler is called with the input still on the 8087 stack, and the 8087 opcode register 
set to 168 hex. 



PL/M-86 usage example: 

Since PL/M-86 does not support a 32-bit integer data type, mqerIA4 cannot be used 
by PL/M programs. You should use either mqerIA2 or mqerlAX. 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• n 3 i r s * 

EXTRN mqerIA4: FAR 

COUNT DD ? 

REAL_C0UNT DQ 100000.5 ; initialized to a test value 

; The following code assumes that the above variables are 
; in the DS segment. 



FLD REAL_C0UNT 

CALL mqerIA4 

MOV BX,0FFSET(C0UNT) 

MOV EBX] f AX 

MOV [BX+23,DX 



load the parameter onto the 
8087 stack 

convert the number to 32 bits 
point to the destination 
move lower 16 bits of answer 
move upper 16 bits of answer 



; With the test input, COUNT is now 100001, which is 
; 000186A1 hex. 
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mqerlAX x = roundaway(x) 

Input parameter: (x) is the top number on the 8087 stack. 

Function: If x is not normal it is first normalized. Then x is rounded to the nearest 
integer. If there are two equally near integers (i.e., the fractional part of x is .5), then 
the integer farthest from zero is selected. For example, 3.3 returns 3; 4.5 returns 5; 
-6.5 returns -7. 

Infinite values are returned unchanged, with no error. 
Output: The answer replaces the input on the 8087 stack. 



Errors: The "I" error is set if the input is a NaN. The NaN is left unchanged. If the 
"I" error is unmasked, the trap handler is called, with the 8087 opcode register set 
to 167 hex. 



PL/M-86 usage example: 

mqerlAX: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerlAX; 

DECLARE (HOURS, MINUTES) REAL; 

HOURS = 2.71; /* Test value */ 

/* The following statement converts HOURS into an integer 
number of MINUTES */ 

MINUTES = mqerIAX(HOURS*60.); 

/* Now MINUTES = 163. */ 

ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• d a i r s * 

EXTRN mqerlAX: FAR 

HOURS DO 2.71 ; initialized to test value 
MINUTES DC) ? 
SIXTY DD 60.0 

The following lines convert HOURS into an integer number 
of MINUTES „ It does the same thing as the above PL/M 
example, only with L0NG_REAL numbers. 



FLD HOURS 
FMUL SIXTY 
CALL mqerlAX 
FSTP MINUTES 



put HOURS onto 8087 stack 

convert to a real number of MINUTES 

round to the nearest integer 

store the integer in L0NG_REAL format 



; With the test case, MINUTES is now 163.0 
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mqerlC2 AX = chop(x) 

Input parameter: (x) is the top number on the 8087 stack. 



Function: If x is not normal it is first normalized. Then if x is an integer, it is 
returned unchanged. If x is not an integer, the fractional part of x is chopped. Thus 
the nearest integer in the direction of is returned. The answer is given in 16-bit 
two's complement form. For example, 4 returns hex 0004; 11.7 returns hex 000B; 
-6.9 returns hex FFFA, which is -6 in two's complement form. 



Output: The input is popped from the 8087 stack, and the answer is left in the AX 
register. 



Errors: The only numbers that will fit the 16-bit destination are those between, but 
not including, -32769 and +32768. Numbers not in this range, including infinities 
and NaN's, will cause an "I" error. If "I" is masked, the "indefinite integer" value 
8000 hex is returned. If "I" is unmasked, the trap handler is called with the input 
still on the 8087 stack. 

Note that "indefinite integer" can also represent the legal output value -32768. The 
two outputs can be distinguished only by the setting of the "I" exception bit of the 
8087. 

If a truncation does take place, the "P" error is set. If "P" is masked the answer is 
returned as usual, since this is not usually considered an error. If "P" is unmasked, 
the trap handler is called. Since the output in the AX register is likely to be lost 
before the trap handler can use it, we recommend that you do not unmask the "P" 
exception. 

If either trap handler is called, the 8087 opcode register is first set to 17E hex. 



PL/M-86 usage example: 

mqerIC2: PROCEDURE (X) INTEGER EXTERNAL; 

DECLARE X REAL; 
END mqerIC2; 

DECLARE C0NTR0L_SETTING INTEGER; 
DECLARE REAL_INPUT REAL; 

REAL_INPUT = 37.885; /* Test value */ 

/* The following line translates REAL_INPUT, which could 
have been calculated using floating point arithmetic, 
into an INTEGER value C0NTR0 L_S ETT I N G , which might be 
output as up to 16 logic lines to a physical device. */ 

C0NTR0L_SETTING = mqer I C2 (REAL_INPUT) ; 

/* For the test input, CONTROL SETTING is now 37; which is 
0025 hex. */ 



4-26 



IC2 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• d a i r s * 

EXTRN mqerIC2: FAR 



REAL_INPUT 
CONTROL SETTING 



DQ 65.7 
DW 



initialized to a test 
va Lue 



The following lines convert the REA L_ INPUT into the 

16-bit integer value CONTRO L_S ETT I NG , just as in the 

PL/M example above, except that REA L_I N PUT is a 
LONG REAL number. 



FLD REAL_INPUT 

CALL mqerIC2 

MOV CONTROL_SETTING,AX 



load the input onto the 8087 

stack 

chop to the integer part 

store the 16-bit answer 



For the input test value, CONTROL_SETTING is now 65, or 
41 hex. 
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mqerlC4 DXAX = chop(x) 

Input parameters: x is the top number on the 8087 stack. 



Function: If x is not normal it is first normalized. Then if x is an integer, it is 
returned unchanged. If x is not an integer, the fractional part of x is chopped. Thus 
the nearest integer in the direction of is returned. The answer is given in 32-bit 
two's complement form. For example, 4 returns hex 00000004; 11.7 returns hex 
0000000B; -6.9 returns hex FFFFFFFA, which is -6 in two's complement form. 



Output: The input is popped from the 8087 stack, and the answer is left in the DX 
and AX registers, with DH the highest byte and AL the lowest byte. 



Errors: The only numbers that will fit the 32-bit destination are those between, but 
not including, -2147483649 and +2147483648. Numbers not in this range, including 
infinities and NaN's, will cause an "I" eriror. If "I" is masked, the "indefinite 
integer" value 80000000 hex is returned. If "I" is unmasked, the trap handler is 
called with the input still on the 8087 stack. 

Note that "indefinite integer" can also represent the legal output value 
-2147483648. The two outputs can be distinguished only by the setting of the "I" 
exception bit of the 8087. 

If a truncation does take place, the "P" error is set. If "P" is masked the answer is 
returned as usual, since this is not usually considered an error. If "P" is unmasked, 
the trap handler is called. Since the output in the DXAX registers is likely to be lost 
before the trap handler can use it, we recommend that you do not unmask the "P" 
exception. 

If either trap handler is called, the 8087 opcode register is first set to 179 hex. 



PL/M-86 usage example: 

Since PL/M-86 does not support a 32-bit integer data type, mqerIC4 cannot be used 
by PL/M programs. You should use either mqerIC2 or mqerlCX. 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

■ d a i r s * 

EXTRN mqerIC4: FAR 

POPULATION DQ 5306279.0 ; total number of voters 

SUPP0RT_SHARE DQ .39 ; proportion of population 

; supporting us 

; above numbers are initialized to test 

; values 
SUPPORT VOTES DD ? 
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The following lines calculate the number of voters 
supporting an issue, given the total population, and the 
fractional share of the population which supports the 
issue. This is simply the share multiplied by the 
total, then chopped to a 32-bit integer SUPP0RT_V0TES . 



FLD POPULATION 

FMUL SUPPORT_SHARE 

CALL mqerIC4 

MOV WORD PTR SUP P0RT_V0T E S , AX 

MOV WORD PTR (SUPP0RT_V0TES+2) , DX 



load total onto 

8087 stack 

multiply by share 

chop to 32 bits 

store bottom 16 

bits 

store top 16 bits 



; With the test inputs, SUPPORT VOTES is now 2069448, 
; which is 001F93C8 hex. 
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mqerlCX x = chop(x) 

Input parameter: (x) is the top number on the 8087 stack. 

Function: If x is not normal it is first normalized. Then if x is an integer, it is 
returned unchanged. If x is not an integer, the fractional part of x is chopped. Thus 
the nearest integer in the direction of is returned. For example, 4 returns 4; -3.9 
returns -3; 1.7 returns 1 . 

Infinite values are returned unchanged, with no error. 
Output: The answer replaces the input on the 8087 stack. 



Errors: If a truncation does take place, the "P" error is set. The correct answer is 
on the 8087 stack. If "P" is unmasked (this is rarely done), the trap handler is then 
called. 

The "I" error is set if the input is a NaN. The NaN is left unchanged, and if the "I" 
error is unmasked, the trap handler is called. 

If either trap handler is called, the 8087 opcode register is first set to 166 hex. 



PL/M-86 usage example: 

mqerlCX: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerlCX; 

DECLARE PRICE_D0LLARS REAL; 

PRICE_DOLLARS = 37.596; /* Test value */ 

/* The following statement chops PR I C E_D0 L L AR S to an even 
number of pennies. First, PR I C E_D0 LLARS is multiplied 
by 100 to get a number of pennies; second, the pennies 
are chopped by mqerlCX; third, the answer is divided by 
100 to convert back into dollars. */ 

PRICE_D0LLARS = ( mqerlCX (PRICE_D0LLARS * 100.) / 100. ) ; 

/* Now PRICE DOLLARS = 37.59 */ 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

pai rs : 

EXTRN mqerlCX: FAR 

PRICE_D0LLARS DQ 37.596 ; initialized to a test 

; v a I u e 
0NE_HUNDRED DD 100.00 ; constant which is used 

; twice below 
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The following lines chop PR I C E_D0 LLARS to an even number 
of pennies, just as in the PL|M example above, except 
that PRICE DOLLARS is here a LONG REAL variable. 



FLD PRICE_.DOLLARS 

FMUL ONE_HUNDRED 
CALL mqerlCX 
FDIV ONE HUNDRED 



amount is loaded onto 8087 

stack 

amount is converted to pennies 

pennies are chopped 

amount is converted back to 

dollars 

converted number is stored, 

8087 stack is popped 



FSTP PRICE_DOLLARS 
; With the above test input, PR I CE_DOLLARS is now 37.59 
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mqerlE2 AX = roundeven(x) 

Input parameter: (x) is the top number on the 8087 stack. 



Function: If x is not normal it is first normalized. Then x is rounded to the nearest 
integer. If the are two integers equally near (i.e., the fractional part of x is .5), then 
the even integer is selected. The answer is given in 16-bit two's complement form. 
For example, 3.1 returns hex 0003; 10.5 returns hex 000A; -2.5 returns hex FFFE, 
which is -2 in two's complement form. 



Output: The input is popped from the 8087 stack, and the answer is left in the AX 
register. 



Errors: The only numbers that will fit the 16-bit destination are those greater than 
or equal to -32768.5 and less than +32767.5. Numbers not in this range, including 
infinities and NaN's, will cause an "I" error. If "I" is masked, the "indefinite 
integer" value 8000 hex is returned. If "I" is unmasked, the trap handler is called 
with the input still on the 8087 stack. 

If a rounding does take place, the "P" error is set. If "P" is masked the answer is 
returned as usual, since this is not usually considered an error. If "P" is unmasked, 
the trap handler is called. Since the output in the AX register is likely to be lost 
before the trap handler can use it, we recommend that you do not unmask the "P" 
exception. 

When the trap handler is called, the 8087 opcode register is first set to 180 hex. 



PL/M-86 usage example: 

mqerIE2: PROCEDURE (X) INTEGER EXTERNAL; 

DECLARE X REAL; 
END mqerIE2; 

DECLARE REAL_VAR REAL; 
DECLARE INTEGER_VAR INTEGER; 

REAL_VAR = 4.5; /* Test value */ 

INTEGER_VAR = mqer I E2 (REAL_VAR) ; 

/* Now in the test case, INTEGER_VAR = 0004 hex. */ 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

* d a i r s * 

EXTRN mqerIE2: FAR 
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INTEGER_VAR 
REAL VAR 



DW ? 
DQ -8.5 



FLD REAL_VAR 
CALL mqerIE2 
MOV INTEGER_VAR,AX 



IE2 



Initialization is a test value 

load the parameter 
round to nearest integer 
store the answer 



; For the test case, I NTEG ER_VAR now equals -8, which 
: is FFF8 hex . 
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mqerlE4 DXAX = roundeven(x) 

Input parameter: (x) is the top number on the 8087 stack. 



Function: If x is not normal it is first normalized. Then x is rounded to the nearest 
integer. If the are two integers equally near (i.e., the fractional part of x is .5), then 
the even integer is selected. The answer is given in 32-bit two's complement form. 
For example, 3.1 returns hex 00000003; 10.5 returns hex 0000000A; -2.5 returns hex 
FFFFFFFE, which is -2 in two's complement form. 



Output: The input is popped from the 8087 stack, and the answer is left in the DX 
and AX registers, with DH the highest byte and AL the lowest byte. 



Errors: The only numbers that will fit the 32-bit destination are those greater than 
or equal to -2147483648.5 and less than +2147483647.5. Numbers not in this range, 
including infinities and NaN's, will cause an "I" error. If "I" is masked, the 
"indefinite integer" value 80000000 hex is returned. If "I" is unmasked, the trap 
handler is called with the input still on the 8087 stack. 

If a rounding does take place, the "P" error is set. If "P" is masked the answer is 
returned as usual, since this is not usually considered an error. If "P" is unmasked, 
the trap handler is called. Since the output in the DXAX registers is likely to be lost 
before the trap handler can use it, we recommend that you do not unmask the "P" 
exception. 

When the trap handler is called, the 8087 opcode register is first set to 17B hex. 

PL/M-86 usage example: 

Since PL/M-86 does not support a 32-bit integer data type, mqerIE4 cannot be used 
by PL/M programs. You should use either mqerIE2 or mqerlEX. 

ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

pai rs : 

EXTRN mqerIE4: FAR 

COUNT DD ? 

REAL_C0UNT DQ 100000.5 ; initialized to a test value 

; The following code assumes that the above variables are 
; in the DS segment. 



FLD REAL_C0UNT 

CALL mqerIE4 

MOV WORD PTR COUNT, AX 

MOV WORD PTR (C0UNT+2), DX 



load the parameter onto 

the 8087 stack 

convert the number to 32 

bits. 

move lower 16 bits of 

answer 

move upper 16 bits of 

answer 



; With the test input, COUNT is now 100000, which is 
; 000186A0 hex. 
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mqerlEX x = roundeven(x) 

Input parameter: (x) is the top number on the 8087 stack. 

Function: If x is not normal it is first normalized. Then x is rounded to the nearest 
integer. If the are two integers equally near (i.e., the fractional part of x is .5), then 
the even integer is selected. For example, 3.3 returns 3; 4.5 returns 4; -6.5 returns 
-6. 

Infinite values are returned unchanged, with no error. 
Output: The answer replaces the input on the 8087 stack. 



Errors: If a rounding does take place, the "P" error is set. The correct answer is on 
the 8087 stack. If "P" is unmasked (this is rarely done), the trap handler is then 
called. 

The "I" error is set if the input is a NaN. The NaN is left unchanged, and if the "I" 
error is unmasked, the trap handler is called. 

If either trap handler is called, the 8087 opcode register is first set to 178 hex. 



PL/M-86 usage example: 

mqerlEX: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerlEX; 

DECLARE UNITS REAL; 
DECLARE THOUSANDS REAL; 

UNITS = 4500.00; /* Test value */ 

/* The follow i rig line computes an integer number of 
thousands, given an input number of units. */ 

THOUSANDS = mqerlEX (UNITS/1 000 .) ; 

/* With the test input value, THOUSANDS now equals 4.00 */ 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• n 3 i p S * 

EXTRN mqerlEX: FAR 

THOUSANDS DQ ? 

UNITS DQ 5890.14 ; initialization is a test 

; value 

0NE_GRAND DD 1000.00 ; constant for the division 

; below 
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The following lines compute an integer number of 
thousands, just as in the PL|M example above, except the 
LONG REAL variables are used. 



FLD UNITS 
FDIV 0NE_GRAND 
CALL mqerlEX 
FSTP THOUSANDS 



load the input onto the 8087 stack 

convert to thousands 

round to an integer value 

store it into 8086 memory and pop 

the 8087 stack 



With the test value, THOUSANDS is now 6.00 
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mqerLGD x = common log(x) 

Input parameter: (x) is the top number on the 8087 stack. 



Function: mqerLGD returns the number L such that 10 L equals the input number. 
For positive inputs this is a well-defined number. For example, an input 10 gives 1; 
100 gives 2; .001 gives -3. Inputs between 10 and 100 give outputs between 1 and 2. 

If the 8087 processor is in affine (signed infinity) mode, the input +INFINITY is 
also valid, returning itself as an answer. 



Output: The answer replaces the input on the 8087 stack. 



Errors: Negative numbers (including -INFINITY), all NaN's, and all unnormals 
give an "I" error. Also, when the 8087 is in projective (unsigned infinity) mode, 
both values of INFINITY give an "I" error. If "I" is masked, the result is the input 
for NaN's; the result is the value INDEFINITE for other invalid inputs. If "I" is 
unmasked, the trap handler is called with the input number still on the 8087 stack. 

Zero input, of either sign, gives a "Z" error. If "Z" is masked, the result is 
-INFINITY. If "Z" is unmasked, the trap handler is called with the input number 
still on the 8087 stack. 

If the input is a denormal, mqerLGD tests the 8087 "D" exception bit to see if the 
8087 is in normalizing mode. If in normalizing mode ("D" unmasked), the input is 
treated as a zero, and a "Z" error is issued. If not in normalizing mode ("D" 
masked), the input is treated as an unnormal, and an "I" error is issued. Note that 
even though the "D" masking bit is tested, the "D" error is never issued by 
mqerLGD, and the "D" trap handler is never called during mqerLGD. 

Whenever the trap handler is called, the number 16D hex is first placed into the 8087 
opcode register. 



PL/M-86 usage example: 

mqerLGD: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerLGD; 

DECLARE QUANTITY REAL; 
DECLARE TENS._POWER REAL; 

QUANTITY = 1900.00; /* Test value */ 

TENS_P0WER = mqerLGD(QUANTITY); 

/* Since the test value QUANTITY is between 10 ** 3 and 
10 ** 4, the answer TEST_P0WER is between 3 and 4. It 
is about 3.27875. */ 



4-37 



LGD 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

■ n 3 i r s * 

EXTRN mqerLGD: FAR 



QUANTITY 
TENS POWER 



DQ .0001 
DQ ? 



initialized to a test value 



The following code implements the above PL/M call in 
assembly language, with L0NG_REAL variables. 



FLD QUANTITY 
CALL mqerLGD 
FSTP TENS POWER 



load input onto 8087 stack 
take the base 10 logarithm 
store the answer and pop the 
8087 stack 



Since the test input was 10 * * - 4 , the output should be 
-4.00. Due to accumulated rounding errors, it may not 
be the exact integer. 
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mqerLGE x = natural log(x) 

Input parameter: (x) is the top number on the 8087 stack. 



Function: mqerLGE returns the number L such that e L equals the input number. 
The constant e is 2.718281828459+. This logarithm is called "natural" because it 
occurs in calculus as the inverse derivative of 1/X; and on many computers it is 
slightly easier to compute than other logarithms. For positive inputs, the logarithm 
is a well-defined number. If the 8087 processor is in affine (signed infinity) mode, 
the input +INFINITY is also valid, returning itself as an answer. 



Output: The answer replaces the input on the 8087 stack. 



Errors: Negative numbers (including -INFINITY), all NaN's, and all unnormals 
give an "I" error. Also, when the 8087 is in projective (unsigned infinity) mode, 
both values of INFINITY give an "I" error. If "I" is masked, the result is the input 
for NaN's; the result is the value INDEFINITE for other invalid inputs. If "I" is 
unmasked, the trap handler is called with the input number still on the 8087 stack. 

Zero input, of either sign, gives a "Z" error. If "Z" is masked, the result is 
-INFINITY. If "Z" is unmasked, the trap handler is called with the input number 
still on the 8087 stack. 

If the input is a denormal, mqerLGE tests the 8087 "D" exception bit to see if the 
8087 is in normalizing mode. If in normalizing mode ("D" unmasked), the input is 
treated as a zero, and a "Z" error is issued. If not in normalizing mode ("D" 
masked), the input is treated as an unnormal, and an "I" error is issued. Note that 
even though the "D" masking bit is tested, the "D" error is never issued by 
mqerLGE, and the "D" trap handler is never called during mqerLGE. 

Whenever the trap handler is called, the number 16C hex is first placed into the 8087 
opcode register. 



PL/M-86 usage example: 

mqerLGE: PROCEDURE (X) REAL EXTERNAL; 

DECLARE X REAL; 
END mqerLGE; 

DECLARE X REAL; 
DECLARE THETA REAL; 

X = 3.0; /* Test value */ 

/* The following code calculates the inverse hyperbolic 
sine of X. That is, THETA is set to the number whose 
hyperbo lie sine is X . */ 

THETA = mqerLGEC X + FSQRT87(X*X + 1.)); 

/* For the test input, THETA now equals about 1.81845 */ 
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ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• d a i r s * 

EXTRN mqerLGE: FAR 



X 

THETA 

ONE 



DQ -3.0 
DQ ? 
DD 1.00 



initialized to a test value 
constant used below 



The following code calculates the inverse hyperbolic 
sine, just as in the PL/M example above, except with 
L0NG_REAL inputs. 



FLD X 

FMUL ST, ST 

FADD ONE 

FSQRT 

FADD X 

CALL mqerLGE 



input parameter onto 8087 stack 

X squared 

X squared + 1 

X + SQRTQ squared + 1) 

take the natural logarithm -- this 

is the answer 

store it and pop the 8087 stack 



FSTP THETA 
; With the test input THETA is now about -1.81845 
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mqerMAX x = max(x,y) 



Input parameters: (x) is the top number on the 8087 stack; (y) is the next number 
on the 8087 stack. 



Function: mqerMAX returns the greater of the numbers (x) and (y). That is, if the 
8087 FCOM instruction indicates that (x) > (y), then (x) is returned. If (x) < (y), then 
(y) is returned. When (x) and (y) test as equal even though they may be given in 
different formats (for example, +0 and -0), (x) is returned. 

If the 8087 chip is in affine (signed infinity) mode, either or both inputs can be 
infinite. If the 8087 is in projective (unsigned infinity) mode, then if either input is 
infinite, both must be infinite. In that case, the values test as equal, and (x) is the 
answer. 



Output: The 8087 stack pops once, with the answer replacing the two inputs. 

Errors: The first error which is detected is the "D" error, for unnormal or denor- 
mal inputs. If the 8087 is in warning mode ("D" is masked), the calculation con- 
tinues; but the "D" error bit remains set. If the 8087 is in normalizing mode ("D" is 
unmasked), the "D" trap handler is called. The trap handler is called directly from 
the interior of mqerMAX; the 8087 opcode register contains the code for the FCOM 
instruction which caused the "D" error. The (y) and (x) inputs are left on the 8087 
stack. Most trap handlers will replace the denormal input(s) with normalized 
numbers, reexecute the FCOM or FSUB instruction, and continue through 
mqerMAX. The trap handler provided by EH87.LIB, described in Chapter 5, will 
do this and replace denormals with zero. 

mqerMAX next checks for NaN inputs. If either input is a NaN, an "I" error is 
given. If "I" is masked, the answer returned is the input NaN (the larger NaN if 
both inputs are NaN's). 

An "I" error is also given if the 8087 chip is in projective (unsigned infinity) mode, 
and exactly one of the inputs is infinite. In this case a masked 'I' yields the value 
INDEFINITE for an answer. 

For all "I" errors., when "I" is unmasked, the trap handler is called with the inputs 
still on the 8087 stack, and the 8087 opcode register set to 265 hex. 

PL/M-86 usage example: 

mqerMAX: PROCEDURE (Y,X) REAL EXTERNAL; 

DECLARE (Y,X) REAL; 
END mqerMAX; 

DECLARE POPULATIONSdO) REAL; 
DECLARE LARGEST REAL; 
DECLARE N BYTE; 

/* The following code sets LARGEST to the greatest of the 
ten values in the array POPULATIONS. */ 

LARGEST = POPULATIONS(O); 
DO N = 1 TO 9; 

LARGEST = mqerMAX(LARGEST,P0PULATI0NS(N)); 
END; 
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ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

* d a i r s * 

EXTRN mqerMAX: FAR 



VAR1 
VAR2 



LARGEST 



DQ 
DQ 



DQ 



-5.3 
-7.9 



the above initializations are test 
va Lues 



The following code sets LARGEST to the maximum of VAR1 
and VAR2. 



FLD VAR1 

FLD VAR2 
CALL mqerMAX 
FSTP LARGEST 



load the first parameter onto 
the 8087 stack 
load the second parameter 
stack now contains the maximum 
maximum is stored, and 8087 
stack is popped 



; With the test inputs, LARGEST is now -5.3 
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mqerMIN x = min(x,y) 



Input parameters: (x) is the top number on the 8087 stack; (y) is the next number 
on the 8087 stack. 



Function: mqerMIN returns the lesser of the numbers (x) and (y). That is, if the 
8087 FCOM instruction indicates that (x) < (y), then (x) is returned. If (x) > (y), then 
(y) is returned. When (x) and (y) test as equal even though they may be given in 
different formats (for example, +0 and -0), (x) is returned. 

If the 8087 chip is in affine (signed infinity) mode, either or both inputs can be 
infinite. If the 8087 is in projective (unsigned infinity) mode, then if either input is 
infinite, both must be infinite. In that case, the values test as equal, and (x) is the 
answer. 



Output: The 8087 stack pops once, with the answer replacing the two inputs. 

Errors: The first error which is detected is the "D" error, for unnormal or denor- 
mal inputs. If the 8087 is in warning mode ("D" is masked), the calculation con- 
tinues; but the "D" error bit remains set. If the 8087 is in normalizing mode ("D" is 
unmasked), the "D" trap handler is called. The trap handler is called directly from 
the interior of mqerMIN; the 8087 opcode register contains the code for the FCOM 
instruction which caused the "D" error. The (y) and (x) inputs are left on the 8087 
stack. Most trap handlers will replace the denormal input(s) with normalized 
numbers, reexecute the FCOM or FSUB instruction, and continue through 
mqerMIN. The trap handler provided by EH87.LIB, described in Chapter 5, will do 
this and replace denormals with zero. 

mqerMIN next checks for NaN inputs. If either input is a NaN, an "I" error is 
given. If "I" is masked, the answer returned is the input NaN (the larger NaN if 
both inputs are NaN's). 

An "I" error is also given if the 8087 chip is in projective (unsigned infinity) mode, 
and exactly one of the inputs is infinite. In this case a masked T yields the value 
INDEFINITE for an answer. 

For all "I" errors, when "I" is unmasked, the trap handler is called with the inputs 
still on the 8087 stack, and the 8087 opcode register set to 265 hex. 

PL/M-86 usage example: 

mqerMIN: PROCEDURE (Y,X) REAL EXTERNAL; 

DECLARE (Y f X) REAL; 
END mqerMIN; 

DECLARE POPULATIONSC10) REAL; 
DECLARE SMALLEST REAL; 
DECLARE N BYTE; 

/* The following code sets SMALLEST to the smallest of the 
ten values in the array POPULATIONS. */ 

SMALLEST = POPULATIONS (0) ; 
DO N = 1 TO 9; 

SMALLEST = mqerMIN (SMALLEST , POPULATIONS (N) ) ; 
END; 
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ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• d a i r s • 

EXTRN mqerMIN: FAR 



VAR1 
VAR2 



SMALLEST 



DQ 
DQ 



-5.3 
-7.9 



the above initializations are test 
values 



DQ 



The following code sets SMALLEST to the minimum of VAR1 
and VAR2. 



FLD VAR1 

FLD VAR2 
CALL mqerMIN 
FSTP SMALLEST 



load the first parameter onto 
the 8087 stack 
load the second parameter 
stack now contains the minimum 
minimum is stored, and 8087 
stack is popped 



; With the test inputs, SMALLEST is now -7.9 
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mqerMOD x = (y mod x), same sign as (y) 



Input parameters: (x) is the top number on the 8087 stack; (y) is the next number 
on the 8087 stack. 



Function: mqerMOD returns the answer (y - (x * mqerICX(y/x)). In more intuitive 
terms, this is the "remainder" left when (y) is divided by (x). The answer is always 
exact; there is no roundoff error. 

mqerMOD always returns the same values when (x) is negative as for the corre- 
sponding positive input (-x). Whenever (x) is mentioned in the following 
paragraphs, we are really talking about the absolute value of (x). 

mqerMOD is calculated by subtracting an integer multiple of (x) from the input (y), 
to bring it down to within (x) units of zero. The choice of which integer multiple to 
subtract determines the range of possible values the function can yield. For 
mqerMOD, the integer is mqerICX(y/x), which is the value (y/x) chopped towards 
zero. If (y) is postive, the answer is greater than or equal to 0, and less than (x). If (y) 
is negative, the answer is greater than (-x), and less than or equal to 0. 

For example, suppose (x) equals either -5 or 5. Then for negative y mqerMOD(y,5) 
gives values in the range between -5 and 0. mqerMOD(-7,5) is -2; 
mqerMOD(-10,5) is 0; mqerMOD(-19.99,5) is -4.99. For positive y, 
mqerMOD(y,5) gives values in the range between and 5. mqerMOD(2,5) is 2; 
mqerMOD(45,5) is 0; mqerMOD(44.75,5) is 4,75. 

It is legal to have infinite (x) inputs. In that case, mqerMOD simply returns (y), if (y) 
is finite and normal. If (y) is unnormal, the normalized (y) is returned. If (y) is 
denormal, the result depends on the setting of the 8087's normalization mode, as 
determined by the "D" error masking flag. If in normalizing mode ("D" 
unmasked), the result is with no error. If in warning mode ("D" masked), the 
result is the unchanged denormal (y), also with no error. 

It is often legal to have unnormal and denormal (y) inputs. The cases with infinite (x) 
are discussed in the above paragraph. If (y) is unnormal, and if the normalization of 
(y) does not produce a denormal, then (y) is legal, and the normalized input is used. 



Output: The 8087 stack pops once, with the answer replacing the two inputs. 



Errors: First, the inputs are checked to see if either is a NaN. If so, an "I" error is 
given. If "I" is masked, the input NaN is returned (if both inputs are NaN's the 
larger NaN is returned). 

If (x) is unnormal, denormal, or any zero value, an "I" error is given. Also, if (y) is 
infinite, an "I" error is given. If "I" is masked, the value INDEFINITE is returned. 

If "I" is unmasked for any of the above errors, the trap handler is called with the 
inputs still on the 8087 stack, and the number 269 hex in the 8087 opcode register. 

A "U" error is given when (y) is unnormal, and the normalization of (y) produces a 
denormal. A "U" error is also given if (y) is denormal. If "U" is masked, an answer 
of is returned. If "U" is unmasked, the trap handler is called, with the value 169 
hex placed in the 8087 opcode register; but the inputs are not on the 8087 stack. 
Instead, the correct answer is given, with a "wrapped" exponent. To obtain the 
correct exponent, subtract the decimal number 24576 from the given exponent. 
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PL/M-86 usage example: 

mqerMOD: PROCEDURE (Y,X) REAL EXTERNAL; 

DECLARE (Y f X) REAL; 
END mqerMOD; 

DECLARE Y REAL; 

DECLARE LAST_THREE_DIGITS REAL; 

Y = 456789.00; /* Test value * / 

/* The following line sets LAST_TH R E E_D I G I TS to Y MOD 
1000. If Y is a positive integer, then the answer is 
the number formed by the last three decimal digits of 
Y. */ 

LAST_THREE_DIGITS = mqerMOD (y , 1 000 .) ; 

/* With the test value, LAST_T HR E E_D I G I T S is now 789.00 */ 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

; pairs: 

EXTRN mqerMOD: FAR 



LAST_THREE_DIGITS 
Y 

ONE GRAND 



DQ ? 

DQ 181137.00 

DD 1000.00 



initialization is a 
test value 
constant for 
ca leu lat i on below 



The following code calculates Y MOD 1000, as in the PL/M 
example above, except with L0NG_REAL variables 



FLD Y ; load first parameter onto 8087 

stack 

load modulus 1000 onto 8087 stack 

take the modu lus 
FSTP LAST_THREE_DIG ITS ; store answer and pop the 

; 8087 stack 



FLD 0NE_GRAND 
CALL mqerMOD 



; With the test value, LAST_THRE E_D I G ITS is 137.00 
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mqerRMD x = (y mod x), close to 



Input parameters: (x) is the top number on the 8087 stack; (y) is the next number 
on the 8087 stack. 



Function: mqerRMD returns the answer (y - (x * mqerIEX(y/x)). In more intuitive 
terms, this is the "remainder" left when (y) is divided by (x). The answer is always 
exact; there is no roundoff error. 

mqerRMD always returns the same values when (x) is negative as for the corre- 
sponding positive input (~x). Whenever (x) is mentioned in the following 
paragraphs, we are really talking about the absolute value of (x). 

mqerRMD is calculated by subtracting an integer multiple of (x) from the input (y), 
to bring it down to within (x) units of zero. The choice of which integer multiple to 
subtract determines the range of possible values the function can yield. For 
mqerRMD, the integer is mqerIEX(y/x), which is the value (y/x) rounded to the 
nearest integer. Thus the range of answers is from (-x/2) to (x/2). 

For example, suppose (x) equals either -5 or 5. Then the value of mqerRMD(y,5) 
ranges from -2.5 to 2.5. mqerRMD(-7,5) is -2; mqerRMD(-10,5) is 0; 
mqerRMD(-19.99,5) is +0.01; mqerRMD(2,5) is 2; mqerRMD(4,5) is -1; 
mqerRMD(44.75,5) is -0.25. 

When the input (y) is an odd integer multiple of (x/2), the answer returned by 
mqerRMD can be either (x/2) or (-x/2). The number chosen is determined by 
the convention of mqerlEX, which rounds to the even integer in case of a tie. This 
results in values alternating between (-x/2) and (x/2). Inputs (y) which yield (x/2) 
form the series -C... -7x/2, -3x/2, x/2, 5x/2, 9x/2, ... >. Inputs (y) which yield 
(-x/2) form the series {..,. -5x/2, -x/2, 3x/2, 7x/2, llx/2, ... >. For example, 
mqerRMD(2.5,5) is 2.5; mqerRMD(7.5,5) is -2.5. 

It is legal to have infinite (x) inputs. In that case, mqerRMD simply returns (y), if (y) 
is finite and normal. If (y) is unnormal, the normalized (y) is returned. If (y) is 
denormal, the result depends on the setting of the 8087's normalization mode, as 
determined by the "D" error masking flag. If in normalizing mode ("D" 
unmasked), the result is with no error. If in warning mode ("D" masked), the 
result is the unchanged denormal (y), also with no error. 

It is often legal to have unnormal and denormal (y) inputs. The cases with infinite (x) 
are discussed in the above paragraph. If (y) is unnormal, and if the normalization of 
(y) does not produce a denormal, then (y) is legal, and the normalized input is used. 



Output: The 8087 stack pops once, with the answer replacing the two inputs. 



Errors: First, the inputs are checked to see if either is a NaN. If so, an "I" error is 
given. If "I" is masked, the input NaN is returned (if both inputs are NaN's the 
larger NaN is returned). 

If (x) is unnormal, denormal, or any zero value, an "I" error is given. Also, if (y) is 
infinite, an "I" error is given. If "I" is masked, the value INDEFINITE is returned. 

If "I" is unmasked for any of the above errors, the trap handler is called with the 
inputs still on the 8087 stack, and the number 27A hex in the 8087 opcode register. 
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A "U" error is given when (y) is unnormal, and the normalization of (y) produces a 
denormal. A "U" error is also given if (y) is already denormal. If "U" is masked, 
an answer of is returned. If "U" is unmasked, the trap handler is called, with the 
value 17A hex placed in the 8087 opcode register; but the inputs are not on the 8087 
stack. Instead, the correct answer is given, with a "wrapped" exponent. To obtain 
the correct exponent, subtract the decimal number 24576 from the given exponent. 



PL/M-86 usage example: 

mqerRMD: PROCEDURE (Y,X) REAL EXTERNAL; 

DECLARE (Y,X) REAL; 
END mqerRMD; 

DECLARE TW0_PI LITERALLY ■6.283185307179586476925'; 

/* 2 * PI -- a full circle expressed in radians * / 

DECLARE THETA REAL; /* angle to be reduced */ 

THETA = 6. ; /* Test value */ 

/* The following line reduces THETA to a principal value 
-- a value between -PI and PI. */ 

THETA = mqerRMD(THETA, TW0_PI); 

/* Now THETA is 6 radians, reduced to the principal value: 
about -0.2831853 */ 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• n a i r s * 

EXTRN mqerRMD: FAR 



THETA 



DQ -6.00 



initialization is a test value 



The following code performs the same reduction of an 
angle to a principal value as the PL/M code above, 
except with a L0NG_REAL variable. 



FLD THETA 
FLDPI 

FADD ST, ST 
CALL mqerRMD 
FSTP THETA 



angle parameter onto 8087 stack 

constant PI onto stack 

2 * PI 

modulus is taken 

principal value is stored, 8087 

stack is popped 



; With the test value, THETA is now about 0.2831853 
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mqerSGN x = (y with x's sign) 



Input parameters: (x) is the top number on the 8087 stack; (y) is the next number 
on the 8087 stack. 

Function: If (x) is greater than or equal to zero, mqerSGN returns the absolute 
value of (y). If (x) is less than zero, mqerSGN returns the negative of the absolute 
value of (y). 

The positive absolute value of (y) is returned for all values of (x) which are zeroes or 
pseudo-zeroes; even if (x) is equivalent to -0. 

Unnormal values of (x) are legal. If (x) is not a pseudo-zero, only the sign of (x) is 
relevant to the final answer. 

Infinite values of (x) are allowed. The sign of the infinity determines the sign on the 
answer, even when the 8087 is in projective (unsigned infinity) mode. 

Any input (y) is legal, including NaN's, unnormals, denormals, and infinities. The 
only part of (y) which might be changed upon output is the sign. 



Output: The 8087 stack pops once, with the answer replacing the two inputs. 



Errors: If (x) is a denormal, the "D" error is given by an FTST instruction within 
the interior of mqerSGN. If the 8087 is in warning mode ("D" is masked), 
mqerSGN will use the denormal to determine the sign of the answer. If the 8087 is in 
normalizing mode ("D" is unmasked), the "D" trap handler will be called with the 
input still on the 8087 stack. Most trap handlers will normalize the argument, reper- 
form the FTST instruction, and continue with the computation of mqerSGN. The 
trap handler provided by EH87.LIB, described in Chapter 5, will replace the denor- 
mal with 0. Thus, the absolute value of (y) will be returned by mqerSGN. 

If (x) is a NaN, an "I" error results. If "I" is masked, (x) is returned. If "I" is 
unmasked, the trap handler is called with the inputs still on the 8087 stack; and the 
8087 opcode register set to 264 hex. 



PL/M-86 usage example: 

mqerSGN: PROCEDURE (Y f X) REAL EXTERNAL; 

DECLARE (Y f X) REAL; 
END mqerSGN; 

DECLARE THETA REAL; 

DECLARE Y_C00R REAL; 

DECLARE PI LITERALLY '3.14159265358979323' ; 

Y_C00R = -0.0000001; /* Test value */ 

/* The following code returns either the value PI or -PI. 
If Y_C00R is positive, it returns PI. If Y_C00R is 
negative, it returns -PI. */ 

THETA = mqerSGN(PI f Y_C00R); 

/* With the test value, THETA now is -PI. */ 
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ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

* d a i r s * 

EXTRN mqerSGN: FAR 



THETA 
Y COOR 



DQ ? 

DQ -0.001 



The following code -0.001 
va lue 



initialized to a test value 
; initialized to a test 



The following code returns PI with the sign of Y_C00R, 
just as in the PL/M example above. 



FLDPI 

FLD Y COOR 



CALL mqerSGN 
FSTP THETA 
; With the test case, THETA is now PI 



first parameter PI onto 8087 stack 

second parameter Y_C00R onto 8087 

stack 

combine sign of Y_C00R with magnitude 

of PI 

store answer and pop the 8087 stack 
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mqerSIN x = sine(x) 

Input parameter: (x) is the top number on the 8087 stack. 



Function: mqerSIN returns the trigonometric sine of x, where x is an angle 
expressed in radians. All input zeroes, pseudo-zeroes, and denormals return the 
input value. Also, unnormals whose value is less than 2~ 63 return the input value. 



Output: The answer replaces the input on the 8087 stack. 



Errors: An "I" error is given for input infinities and NaN's. An "I" error is also 
given for unnormals which do not represent values less than 2~ 63 . 

If "I" is unmasked, the trap handler is called with the input still on the stack, and 
the 8087 opcode register set to 171 hex. If "I" is masked, the answer is the input for 
NaNs; the answer is the value INDEFINITE for other invalid inputs. 



PL/M-86 usage example: 

mqerSIN: PROCEDURE (THETA) REAL EXTERNAL; 

DECLARE THETA REAL; 
END mqerSIN; 

DECLARE (P0LAR_R f POLAR_THETA) REAL; 
DECLARE REC_Y REAL; 

DECLARE PI LITERALLY '3.141 59265358979 ' ; 
DECLARE DEG_T0_RAD LITERALLY °PI/180.'; 

P0LAR_R = 2.; POLAR_THETA = 30.; /* Test values */ 

/* The following line computes the Y-coordinate of a 

polar-to-rectangular conversion. The input angle is in 
degrees, so it must be converted to radians. */ 

REC_Y = P0LAR_R * mqerSIN <P0LAR_THETA * DEG_T0_RAD); 

/* Now in the test case, REC_Y = 1. */ 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• d a i r s * 

EXTRN mqerSIN: FAR 

P0LAR_THETA DQ 30.0 
P0LAR_R DQ 2.0 

; the above initializations are test 

; values. 
REC_Y DQ ? 

DEG_T0_RAD DT 3F F98EFA351 294E9C8AER ; the constant 

; PI/180. 
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The following lines c omp ute the Y-coordinate of 
polar-to-rectangular conversion, as in the PL/M 
example above; except that the variables are 
LONG REAL. 



FLD POLAR_THETA 

FLD DEG_T0_RAD 

FMUL 

CALL mqerSIN 

FMUL P0LAR_R 

FSTP REC Y 



degrees angle onto 8087 stack 

convert ed to radians 

sine is taken 

answer scaled to correct radius 

Y-coordinate stored and stack is 

popped 



With the test case, REC_Y is now 1 
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mqerSNH x = hyperbolic sine(x) 

Input parameter: (x) is the top number on the 8087 stack. 

Function: mqerSNH returns the hyperbolic sine of x, where x is an angle expressed 
in radians. All input zeroes, pseudo-zeroes, and denormals return the input value. 
Also, unnormals whose value is less than 2~ 63 return the input value. 

Infinite inputs are legal, and return the input value. 
Output: The answer replaces the input on the 8087 stack. 



Errors: An "I" error is given for input NaN's. An "I" error is also given for 
unnormals which do not represent values less than 2~ 63 . 

If "I" is unmasked, the trap handler is called with the input still on the stack. If "I" 
is masked, the answer is the input for NaNs; the answer is the value INDEFINITE 
for other invalid inputs. 

mqerSNH will give an "O" overflow error if the input is greater than about 1 1355. 
When "O" is masked, the value +INFINITY is returned. Likewise, "O" is given for 
inputs less than about -11355, with -INFINITY returned for masked "O". When 
"O" is unmasked, the trap handler is called with the input still on the 8087 stack. 

When either trap handler is called, mqerSNH first sets the 8087 opcode register to 
16E hex. 



PL/M-86 usage example: 

mqerSNH: PROCEDURE (THETA) REAL EXTERNAL; 

DECLARE THETA REAL; 
END mqerSNH; 

DECLARE (INPUT_VALUE, OUTPUTJ/ALUE) REAL; 

INPUT_VALUE = 2.7; /* Test value */ 

OUTPUTJ/ALUE = mqe rSNH ( INPUTJ/ALUE) ; 

/* Now with the test input, 0UTPUT_VALUE is about 
14.812526 */ 

ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• d a i r s * 

EXTRN mqerSNH: FAR 

INPUT_VALUE DQ -2.7 ; initialization is a test 

; value 
OUTPUT VALUE DQ ? 



4-53 



SNH 



The following code duplicates the above PL|M 
assignment statement, 

except with L0NG_REAL variables. 



FLD INPUT_VALUE 

CALL mqerSNH 
FSTP OUTPUT VALUE 



load the parameter onto the 8087 

stack 

take the hyperbolic sine 

store the answer and pop the 

8087 stack 



; With the test input, OUTPUT_VALUE is now about 
; -14.812526 
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mqerTAN x = tangent(x) 

Input parameter: (x) is the top number on the 8087 stack. 

Function: mqerTAN returns the trigonometric tangent of x, where x is an angle 
expressed in radians. All input zeroes, pseudo-zeroes, and denormals return the 
input value. Also, unnormals whose value is less than 2~ 63 return the input value. 

Output: The answer replaces the input on the 8087 stack. 

Errors: An "I" error is given for input infinities and NaN's. An "I" error is also 
given for unnormals which do not represent values less than 2~ 63 . 

If "I" is unmasked, the trap handler is called with the input still on the stack, and 
the 8087 opcode register set to 173 hex. If "I" is masked, the answer is the input for 
NaNs; the answer is the value INDEFINITE for other invalid inputs. 

A "Z" error is given when the input number is an exact odd multiple of the closest 

TEMP REAL number to n/2. When "Z" is masked, the answer 

+INFINITY is returned. When "Z" is unmasked, the trap handler is called with the 
input still on the 8087 stack, and the 8087 opcode register set to 173 hex. 

PL/M-86 usage example: 

mqerTAN: PROCEDURE (THETA) REAL EXTERNAL; 

DECLARE THETA REAL; 
END mqerTAN; 

DECLARE PI LITERALLY '3.141 59265358979 * ; 

DECLARE DE6_T0_RAD LITERALLY 'PI/180.'; 
DECLARE THETA_DEGREES REAL; 
DECLARE SLOPE REAL; 

THETA_DE6REES = 135.0; /* Test value */ 

/* The following line computes the tangent of the angle 
THETA_DEGREES. The answer is called SLOPE because it is 
the slope of a line which is displaced by TH ET A_DE6 RE ES 
from the X-axis. */ 

SLOPE = mqerTAN(THETA_DEGREES * DEG_T0_RAD); 

/* Now with the test value, SLOPE = -1. */ 

ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

" d a i r s * 

EXTRN mqerTAN: FAR 

THETA_DEGREES DQ 45.00 ; initialization is a test 

; value 
SLOPE DQ ? 

DEG_T0_RAD DT 3 F F98EFA351 294E9C8AER ; the constant 

; PI/180. 
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The following code computes the tangent just as in 
the PL/M example above, except with LON G_R E A L 
variables. 



FLD THETA_DEGREES 

FLD DEG_T0_RAD 

FMUL 

CALL mqerTAN 



FSTP SLOPE 
; With the test input, SLOPE is now 1.00 



load the first parameter onto 
the 8087 stack 

convert from degrees to radians 

take the tangent of the radians 

value; 

store the answer and pop the 

8087 stack 
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mqerTNH x = hyperbolic tangent(x) 

Input parameter: (x) is the top number on the 8087 stack. 



Function: mqerTNH returns the hyperbolic tangent of x, where x is an angle 
expressed in radians. All input zeroes, pseudo-zeroes, and denormals return the 
input value. Also, unnormals whose value is less than 2~ 63 return the input value. 

Infinite inputs are allowed. Input +INFINITY gives an answer of +1; -INFINITY 
gives -1. The sign of the infinity is significant even if the 8087 is in projective 
(unsigned infinity) mode. 



Output: The answer replaces the input on the 8087 stack. 



Errors: An "I" error is given for input NaN's. An "I" error is also given for 
unnormals which do not represent values less than 2" 63 . 

If "I" is unmasked, the trap handler is called, with the input still on the stack, and 
the 8087 opcode register set to 170 hex. If "I" is masked, the answer is the input for 
NaNs; the answer is the value INDEFINITE for illegal unnormals. 



PL/M-86 usage example: 

mqerTNH: PROCEDURE (THETA) REAL EXTERNAL; 

DECLARE THETA REAL; 
END mqerTNH; 

DECLARE (INPUT_VALUE, OUTP.UT_VALUE) REAL; 

INPUT_VALUE = 0.62; /* Test value */ 

0UTPUT_VALUE = mqerTNH ( INPUT_VALUE) ; 

/* Now with the test input, 0UTPUT_VALUE is about 
0.55112803 */ 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• d a i r s * 

EXTRN mqerTNH: FAR 

INPUT_VALUE DQ -0.62 ; initialization is a test 

; value 
0UTPUT_VALUE DQ ? 

; The following code duplicates the above PL/M 
; assignment statement, except with L0NG_REAL 
.•variables. 
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FLD INPUT_VALUE 

CALL mqerTNH 
FSTP OUTPUT VALUE 



Load the parameter onto the 8087 

stack 

take the hyperbolic tangent 

store the answer and pop the 

8087 stack 



; With the test input, OUTPUT_VALUE is now about 
; -0.55112803 
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mqerY2X x = y x 



Input parameters: (x) is the top number on the 8087 stack; (y) is the next number 
on the 8087 stack. 



Function: mqerY2X computes (y) to the (x) power; neither (y) nor (x) is required to 
be an integer. For most inputs, the formula used is 2( X * LG2 M). 

The base (y) must be positive for the logarithmic formula to have any meaning. 
There are some cases, however, in which (y) is allowed to be non-positive. 

• If (x) is positive, a zero (y) gives a zero answer, with no error. 

• If (x) is zero, (y) can be negative; the answer is 1 . 

• If (x) is an integer, (y) can be negative, The function is evaluated using the 
absolute value of (y). The result is the answer if (x) is even; it is the negative of 
the answer if (x) is odd. 

Inputs which are zero, infinite, unnormal, or denormal are accepted under certain 
conditions. 

When either input is denormal, it is replaced with an alternate value. The value 
selected depends on the setting of the 8087 "D" masking bit. If the 8087 is in nor- 
malizing mode ("D" is unmasked), the input is replaced by 0. If the 8087 is in warn- 
ing mode ("D" is masked), the input is replaced by the equivalent unnormal. Note 
that even though mqerY2X references the "D" masking bit, it never gives a "D" 
error, and it never calls the "D" trap handler. 

An unnormal (y) input is legal only if (x) is a normal integer or infinite value. If (x) is 
infinite, the function is evaluated as if (y) were the equivalent normal value. If (x) is 
zero, (y) must be nonzero; in that case, the answer is 1. If (x) is any other integer, it 
must fit into 32 bits; in that case, the function mqerYI4 is called to obtain the 
answer. 

An unnormal (x) input is legal only if it is non-zero, and if (y) is infinite or zero. In 
those cases, (x) is replaced by its normal equivalent. 

When the 8087 is in affine (signed infinity) mode, there are a number of cases in 
which mqerY2X allows infinite inputs: 

• If (y) is -INFINITY, then (x) must be a non-zero integer. The magnitude of the 
answer is then INFINITY if (x) is positive; zero if (x) is negative. The sign of the 
answer is positive if (x) is even; negative if (x) is odd. 

• If (y) is +INFINITY, then any non-zero (x) is legal. The answer is +INFINITY 
if (x) is positive; zero if (x) is negative. 

• If (x) is +INFINITY, then (y) must be positive or any zero, and not equal to 1. 
The answer is zero if (y) is less than 1 , and +INFINITY if (y) is greater than 1 . 

• If (x) is -INFINITY, then (y) must likewise be positive or any zero, and not 
equal to 1. The answer is +INFINITY if (y) is less than 1, and -INFINITY if (y) 
is greater than 1 . 

When the 8087 chip is in projective (unsigned infinity) mode, there is only one case 
in which any infinite input is allowed. This occurs when (y) is infinite, and (x) is a 
non-zero integer. The result is the same as if the 8087 were in affine mode. 



Output: The 8087 stack pops once, with the answer replacing the two inputs. 
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Errors: mqerY2X first checks for NaN inputs. If either input is a NaN, an "I" 
error is given. If "I" is masked, the input NaN is returned (the larger NaN is 
returned if both inputs are NaN's). 

The legal cases involving unnormal inputs, infinite inputs, and negative (y) inputs 
are described above. Illegal cases yield an "I" error. If "I" is masked, the value 
INDEFINITE is returned. The case (y = 1) and (x infinite), among others, falls into 
this category. 

It is illegal for both (x) and (y) to have zero values. This too gives an "I" error, with 
an INDEFINITE answer if "I" is masked. 

If (y) is any zero and (x) is any negative value (including negative infinity), then a 
"Z" error is given. If "Z" is masked, the value +INFINITY is returned. 

The "O" overflow or "U" underflow error occurs when (y x ) cannot be represented 

by the TEMP REAL exponent. If "O" is masked, an overflow will return the 

answer +INFINITY. In underflow cases, if "U" is masked, the correct answer is 
given if it can be represented by a denormal; otherwise, is given. 

All of the errors, "I", "Z", "O", and "U", will cause the trap handler to be called 
when the corresponding exception bit is unmasked. mqerY2X leaves the input 
numbers on the 8087 stack, and places the value 26A hex into the 8087 opcode 
register before calling the trap handler. 



PL/M-86 usage example: 

mqerY2X: PROCEDURE (Y,X) REAL EXTERNAL; 

DECLARE (Y,X) REAL; 
END mqerY2X; 

DECLARE CUBE_R00T REAL; 
DECLARE INPUT_VALUE REAL; 

INPUT_VALUE = 17.00; /* Test value */ 

/* The following line takes the cube root of the positive 
value INPUTJ/AL.UE. */ 

CUBE_R00T = mqerY2X(INPUT_VALUE f 1./3.); 

/* With the test input, INPUT_VALUE is now about 
2.5712816 */ 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

• d a i r s * 

EXTRN mqerY2X: FAR 

INPUT_VALUE DQ 64.00 ; initialization is a test 

; value 
CUBE_R00T DQ ? 
0NE_THIRD DT 0.333333333333333333333 ; constant 

; used below 
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; The following lines take the cube root just as in the 
; PL/M example above, except with L0NG_REAL variables. 



FLD INPUT_.VALUE 

FLD 0NE_THIRD 

CALL mqerY2X 
FSTP CUBE ROOT 



load first parameter onto 

8087 stack 

load second parameter onto 

8087 stack 

exponent i ate 

store the answer and pop the 

8087 stack 



Y2X 



; With the test input, CUBE_R00T is now about 4.00 
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mqerYI2 x = x**AX 



Input parameters: (x) is the top number on the 8087 stack. The power to which (x) 
is raised is the 8086 AX register, interpreted as a twos-complement signed integer. 



Function: mqerYI2 raises the real input (x) to an integer power. If the integer is 
zero, the answer is 1 . In this case, no error is given, even if (x) is a NaN. 

If AX is not zero, the input (x) is first checked for unusual values. 

If (x) is +INFINITY, the answer is +INFINITY for positive AX; +0 for negative 
AX. There is no error. 

If (x) is -INFINITY, the magnitude of the answer is INFINITY for positive AX; 
for negative AX. The sign of the answer is positive for even AX; negative for odd 
AX. 

Zero (x) input is legal if AX is positive. The answer is -0 if (x) is -0 and AX is odd; 
the answer is +0 in all other cases. 

If (x) is denormal, no error is given. However, the 8087 "D" error masking bit is 
checked to see what action to take. If the 8087 is in normalizing mode ("D" is 
unmasked), (x) is replaced by zero. If the 8087 is in warning mode ("D" is masked), 
(x) is replaced by the unnormal number with the same numeric value as the denor- 
mal. The evaluation of mqerYI2 proceeds with the new (x). 

If (x) is unnormal and AX is negative, then (1/x) is computed, preserving the 
number of bits of unnormalization. Then the positive power is computed by 
successive squaring and multiplying by (x). 

If (x) is unnormal and AX is positive, then the power is computed by successive 
squaring and multiplying by (x). 

For normal, non-zero values of (x), computation of the power proceeds according to 
the value of the integer power AX. 

If the integer power is 64 or greater, or -64 or less, the answer is computed with 
logarithms. The answer is 2 ** (AX * LG2(x)). 

If the integer power is from 1 to 63, the answer is computed by successively squaring 
and multiplying by (x) to achieve the correct power. 

If the integer power is from -63 to -1, mqerYI2 determines if any exceptions would 
occur if the expression 1 / (x * x * ... x) were evaluated. If not, the expression is 
evaluated and the answer is returned. If so, then the expression (1/x) * (1/x) * ... * 
(1/x) is evaluated. If the second expression causes exceptions, the trap handler is 
called. 

The maximum number of multiplications performed for any of the above squaring- 
and-multiplying algorithms is 9. 



Output: The answer replaces the input (x) on the 8087 stack. The AX input is 
destroyed, as allowed by PL/M-86 procedure conventions. 
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Errors: As stated above, there can be no errors if AX is 0. Otherwise, errors occur 
in the following cases: 

If (x) is a NaN, an "I" error is given. If "I" is masked, the input NaN is returned as 
the answer. 

A "U" underflow error occurs when the answer computed is too close to zero to be 
represented by the exponent. If "U" is masked, the answer is replaced by the 
equivalent denormal if it exists; otherwise. 

An "O" overflow error occurs when the magnitude of the answer is too great to be 
represented by the exponent. If "O" is masked, the answer is INFINITY with the 
appropriate sign. 

If any of the "I", "O", or "U" errors occurs with the error unmasked, the trap 
handler is called. Before calling the trap handler, mqerYI2 sets the 8087 opcode 
register to 27C hex and leaves the inputs on the 8087 stack. The integer power is con- 
verted to TEMP REAL and pushed onto the top of the 8087 stack. The base (the 

original (x)) becomes the second stack element. 



PL/M-86 usage example: 

For PL/M-86, the 16-bit input parameter should be on the 8086 stack rather than in 
the AX register. Therefore, use mqerYIS (instead of mqerYI2) in PL/M-86 
programs. 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

; pairs: 

EXTRN mqerYI2: FAR 



POWER 
REAL BASE 



REAL OUTPUT 



DW 9 ; exponent which will be used 
DQ 1.3 ; number which will be raised to 
; POWER 

; the above initializations are test 

; values 
DQ ? 



The following code multiplies POWER copies of REAL_BASE 
together, and stores the answer in REAL_OUTPUT. 



FLD REAL_BASE 
MOV AX, POWER 
CALL mqerYl2 
FSTP REAL OUTPUT 



base parameter goes onto the 8087 

stack 

exponent goes into the AX 

register 

REAL_BASE ** POWER is now on 8087 

stack 

store the answer and pop the 8087 

stack 



; With the test inputs, REAL_0UTPUT is now about 10.604499 
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mqerYI4 x = x**DXAX 



Input parameters: (x) is the top number on the 8087 stack. The power to which (x) 
is raised is a 32-bit twos complement value in the 8086 DX and AX registers. DX is 
the most significant half, and AX is the least significant half. 

Function: mqerYI2 raises the real input (x) to an integer power. The input integer is 
presented in a 32-bit format. Note, however, that 32 bits are rarely necessary to 
represent an integer exponent. Raising a number to a power greater than 32768 
rarely gives meaningful results. Thus mqerYI2 is sufficient for almost every applica- 
tion in which mqerYI4 might be used. We have provided mqerYI4 mainly for com- 
patibility with the 32-bit integer types found in Pascal-86 and FORTRAN-86. 

Ff the input integer power is zero, the answer is 1, no matter what (x) is. In this case, 
no error is given, even if (x) is a NaN. 

If DXAX is not zero, the input (x) is first checked for unusual values. 

If (x) is +INFINITY, the answer is +INFINITY for positive DXAX; +0 for negative 
DXAX. There is no error. 

If (x) is -INFINITY, the magnitude of the answer is INFINITY for positive DXAX; 
for negative DXAX. The sign of the answer is positive for even DXAX; negative 
for odd DXAX. 

Zero (x) input is legal if DXAX is positive. The answer is -0 if (x) is -0 and DXAX 
is odd; the answer is +0 in all other cases. 

If (x) is denormal, no error is given. However, the 8087 "D" error masking bit is 
checked to see what action to take. If the 8087 is in normalizing mode ("D" is 
unmasked), (x) is replaced by zero. If the 8087 is in warning mode ("D" is masked), 
(x) is replaced by the unnormal number with the same numeric value as the denor- 
mal. The evaluation of mqerYI2 proceeds with the new (x). 

If (x) is unnormal and DXAX is negative, then (1/x) is computed, preserving the 
number of bits of unnormalization. Then the positive power is computed by 
successive squaring and multiplying by (x). 

If (x) is unnormal and DXAX is positive, then the power is computed by successive 
squaring and multiplying by (x). 

For normal, non-zero values of (x), computation of the power proceeds according to 
the value of the integer power DXAX. 

If the integer power is 64 or greater, or -64 or less, the answer is computed with 
logarithms. The answer is 2 ** (DXAX * LG2(x)). 

If the integer power is from 1 to 63, the answer is computed by successively squaring 
and multiplying by (x) to achieve the correct power. 

If the integer power is from -63 to -1, mqerYI4 determines if any exceptions would 
occur if the expression 1 / (x * x * ... x) were evaluated. If not, the expression is 
evaluated and the answer is returned. If so, then the expression (1/x) * (1/x) * ... * 
(1/x) is evaluated. If the second expression causes exceptions, the trap handler is 
called. 

The maximum number of multiplications performed for any of the above squaring- 
and-multiplying algorithms is 9. 
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Output: The answer replaces the input (x) on the 8087 stack. The DXAX input is 
destroyed, as allowed by PL/M-86 procedure conventions. 



Errors: As stated above, there can be no errors if DXAX is 0. Otherwise, errors 
occur in the following cases: 

If (x) is a NaN, an "I" error is given. If "I" is masked, the input NaN is returned as 
the answer. 

A "U" underflow error occurs when the answer computed is too close to zero to be 
represented by the exponent. If "U" is masked, the answer is replaced by the 
equivalent denormal if it exists; otherwise. 

A "O" overflow error occurs when the magnitude of the answer is too great to be 
represented by the exponent. If "O" is masked, the answer is INFINITY with the 
appropriate sign. 

If any of the "I", "O", or "U" errors occurs with the error unmasked, the trap 
handler is called. Before calling the trap handler, mqerYI2 sets the 8087 opcode 
register to 27C hex and leaves the inputs on the 8087 stack. The integer power is con- 
verted to TEMP__REAL and pushed onto the top of the 8087 stack. The base (the 
original (x)) becomes the second stack element. 



PL/M-86 usage example: 

Since PL/M-86 does not support a 32-bit integer data type, mqerYI4 cannot be used 
by PL/M programs. You should use either mqerYIS or mqerY2X. 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

* d a i r s * 

EXTRN mqerYU: FAR 

POWER DD 9 ; exponent which wi LL be used 

REAL_BASE DQ 1.3 ; number which- will be raised to 

; POWER 

; the above initializations are 

; test values 
REAL_0UTPUT DQ ? 

; The following code multiplies POWER copies of RE A L_B AS E 
; together, and stores the answer in REAL_0UTPUT. 



FLD REAL_BASE 

MOV AX, WORD PTR POWER 

MOV DX, WORD PTR (POWER+2) 
CALL mqerYH 



base parameter goes onto 

the 8087 stack 

low 16 bits of exponent 

to AX 

---high 16 bits to DX 

REAL_BASE ** POWER is 

now on 8087 stack 

store the answer and pop 

the 8087 stack 



FSTP REAL_.0UTPUT 
; With the test inputs, REAL_0UTPUT is now about 10.604499 
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mqerYIS x = x**STK 



Input parameters: (x) is the top number on the 8087 stack. The power STK to 
which (x) is to be raised is a 16-bit twos-complement integer that is pushed onto the 
stack before mqerYIS is called. 



Function: mqerYIS raises the real input (x) to an integer power. If the integer is 
zero, the answer is 1 . In this case, no error is given, even if (x) is a NaN. 

If STK is not zero, the input (x) is first checked for unusual values. 

If (x) is +INFINITY, the answer is +INFINITY for positive STK; +0 for negative 
STK. There is no error. 

If (x) is -INFINITY, the magnitude of the answer is INFINITY for positive STK; 
for negative STK. The sign of the answer is positive for even STK; negative for odd 
STK. 

Zero (x) input is legal if STK is positive. The answer is -0 if (x) is -0 and STK is odd; 
the answer is +0 in all other cases. 

If (x) is denormal, no error is given. However, the 8087 "D" error masking bit is 
checked to see what action to take. If the 8087 is in normalizing mode ("D" is 
unmasked), (x) is replaced by zero. If the 8087 is in warning mode ("D" is masked), 
(x) is replaced by the unnormal number with the same numeric value as the denor- 
mal. The evaluation of mqerYIS proceeds with the new (x). 

If (x) is unnormal and STK is negative, then (1/x) is computed, preserving the 
number of bits of unnormalization. Then the positive power is computed by suc- 
cessive squaring and multiplying by (x). 

If (x) is unnormal and STK is positive, then the power is computed by successive 
squaring and multiplying by (x). 

For normal, non-zero values of (x), computation of the power proceeds according to 
the value of the integer power STK. 

If the integer power is 64 or greater, or -64 or less, the answer is computed with 
logarithms. The answer is 2 ** (STK * LG2(x)). 

If the integer power is from 1 to 63, the answer is computed by successively squaring 
and multiplying by (x) to achieve the correct power. 

If the integer power is from -63 to -1, mqerYIS determines if any exceptions would 
occur if the expression 1 / (x * x * ... x) were evaluated. If not, the expression is 
evaluated and the answer is returned. If so, then the expression (1/x) * (1/x) * ... * 
(1/x) is evaluated. If the second expression causes exceptions, the trap handler is 
called. 

The maximum number of multiplications performed for any of the above squaring- 
and-multiplying algorithms is 9. 



Output: The answer replaces the input (x) on the 8087 stack. mqerYIS returns with 
the value STK popped off the 8086 stack, so it no longer exists. 
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Errors: As stated above, there can be no errors if STK is 0. Otherwise, errors occur 
in the following cases: 

If (x) is a NaN, an "I" error is given. If "I" is masked, the input NaN is returned as 
the answer. 

A "U" underflow error occurs when the answer computed is too close to zero to be 
represented by the exponent. If "U" is masked, the answer is replaced by the 
equivalent denormal if it exists; the answer is otherwise. 

An "O" overflow error occurs when the magnitude of the answer is too great to be 
represented by the exponent. If "O" is masked, the answer is INFINITY with the 
appropriate sign. 

If any of the "I", "O", or "U" errors occurs with the error unmasked, the trap 
handler is called. Before calling the trap handler, mqerYIS sets the 8087 opcode 
register to 27C hex and leaves the inputs on the 8087 stack. The integer power is con- 
verted to TEMP REAL and pushed onto the top of the 8087 stack. The base (the 

original (x)) becomes the second stack element. 



PL/M-86 usage example: 

mqerYIS: PROCEDURE (Y,I) REAL EXTERNAL; 

DECLARE Y REAL, I INTEGER; 
END mqerYIS; 

DECLARE INTEREST_RATE REAL; 

DECLARE NUMBER_OF_PERIODS INTEGER; 

DECLARE START_AMOUNT REAL; 

DECLARE FINISH_AM0UNT REAL; 

INTEREST_RATE = 0.015; /* Test value */ 
NUMBER_0F_PERI0DS = 12; /* Test value */ 
START_AM0UNT = 1000.00; /* Test value */ 

/* The following line calculates compound interest for the 
given NUMB ER_0 F_PER I 0DS , given the rate I NTEREST_R ATE 
for each period. I NTEREST_R AT E is presented as a 
fraction of 1; for example, the value 0.015 represents 
1.5 percent interest for each time period. */ 

FINISH_AM0UNT = START_AMOUNT* mqerYIS 

(1 . + INTEREST_RATE, NUMBER_0F_PER I0DS ) ; 

/* With the test inputs, FINISH_AM0UNT is now about 
1195.62. This is the balance of an unpaid loan of 
$1000.00 after one year, if the loan accumulates 1.5 
percent interest every month. */ 



ASM-86 usage example: 

; This EXTRN must appear outside of all SEGMENT-ENDS 

; pairs: 

EXTRN mqerYIS: FAR 



4-67 



YIS 



INTEREST_RATE 
NUMBER_OF_PERIODS 
START AMOUNT 



FINISH AMOUNT 



DQ 0.015 

DW 12 

DQ 1000.00 

; the above initializations are test 

; va Lues 

DQ ? 



The following code implements the above PL/M example in 
assembly language, with LONG_REAL variables. 



FLD1 

FADD INTEREST_RATE 

PUSH NUMBER OF PERIODS 



1 onto 8087 stack 
(1 + I) is on stack 
exponent parameter goes onto 
8086 stack 

(1 + I) ** N is on 8087 stack 
scaled up by the amount of 
money 

store result and pop the 8087 
stack 

; NOTE: Do not explicitly POP the NUMB E R_0 F_PE R I ODS from 
; the 8086 stack -- mqerYIS does that for you. 



CALL mqerYIS 
FMUL START_AMOUNT 

FSTP FINISH AMOUNT 
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The Common Elementary Function Library 



Linkage of CEL87.LIB to Your Program Modules 

The final action you must take to use CEL87.LIB in your programs is to include the 
file name CEL87.LIB into an appropriate LINK86 command. 

CEL87 requires either the 8087 component or the 8087 emulator. If the component 
is present, you must also link in 8087. LIB. If the emulator is used, you must link in 
both E8087 and E8087. LIB. 

If you are also using EH87.LIB, you must give the name EH87.LIB after 
CEL87.LIB. If you put EH87.LIB before CEL87.LIB, the program will link with no 
error messages, but it will halt after the first CEL87 function is called. 

Following is the suggested order for object modules in your LINK86 statement. 

Your object modules 

DCON87.LIB if you are using it 

CEL87.LIB 

EH87.LIB if you are using it 

8087. LIB if you are using the component, or 

E8087, E8087.LIB if you are using the emulator 

As an example, if you are linking your PL/M-86 modules MYMODl.OBJ and 
MYMOD2.0BJ into a program using the 8087 emulator and the error handler, issue 
the command 



> 
>> 



LINK86 :F1 : MYM0D1 .OBJ, : F1 : MYMqDZ . OBJ , &<cr> 

3 : FO : CEL87- LIB , : FO : EH87 . LIB , i : F : E 8 8 7 , &<cr> 
:F0:E8087.:LIB TO :F1 : MYPROG . UN K 



If you have a single ASM-86-generated object module :Fl:MYPROG.OBJ to be 
executed in a system with an 8087 chip, issue the command 



>> 



LINK86 : F1 :MYPR0G .OBJ , : FO : CEL87. LIB, &<cr> 
: F : 8 8 7 . L I B TO : F 1 : M Y P R G . L N K 
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CHAPTER 5 
THE ERROR HANDLER MODULE 



Overview 

This chapter describes EH87.LIB, a library of five utility procedures which you can 
use to write trap handlers. Trap handlers are procedures which are called when an 
unmasked 8087 error occurs. 

EH87.LIB also contains a set of "dummy" public symbols for all the functions of 
CEL87.LIB. These symbols save you code when you do not use all the functions of 
CEL87.LIB. Their presence, however, makes it absolutely necessary that you link 
EH87.LIB and CEL87.LIB in the correct order. The last section of this chapter tells 
you how to do so. 

EH87.LIB also contains a set of alternate public names for its procedures, which are 
used by some Intel translators. They are listed in Appendix F. 

The 8087 error reporting mechanism can be used not only to report error conditions, 
but also to let software implement modes and functions not directly supported by 
the chip. This chapter defines three such extensions to the 8087: normalizing mode, 
non-trapping NaN's, and non-ordered comparison. The utility procedures support 
these extra features. 

DECODE is called near the beginning of the trap handler. It preserves the complete 
state of the 8087, and also identifies what function called the trap handler, with what 
arguments and/or results. DECODE eliminates much of the effort needed to deter- 
mine what error caused the trap handler to be called. 

NORMAL provides the "normalizing mode" capability for handling the "D" 
exception (described in the following section). By calling NORMAL in your trap 
handler, you eliminate the need to write code in your application program which 
tests for non-normal inputs. 

SIEVE provides two capabilities for handling the "I" exception. It implements non- 
trapping NaN's and non-ordered comparisons (both described in the following 
sections). These two IEEE standard features reduce the incidence of multiple error 
reports for a single bad input. 

ENCODE is called near the end of the trap handler. It restores the state of the 8087 
saved by DECODE, and performs a choice of concluding actions, by either retrying 
the offending function or returning a specified result. ENCODE provides a common 
path for exiting the trap handler and resuming execution of the user program. 

FILTER calls each of the above four procedures. If your error handler does nothing 
more than detect fatal errors and implement the features supported by SIEVE and 
NORMAL, then your interface to EH87.LIB can be accomplished with a single call 
to FILTER. 

Normalizing Mode 

Normalizing mode allows you to perform floating point operations without having 
to worry about whether the operands are in normal form. All denormal inputs will 
be normalized before the operation, without any user intervention. 

The 8087 provides the "D" error, which warns you that an operand is not normal. 
You can implement normalizing mode in software by unmasking the "D" error, and 
providing a "D" trap handler. The handler should perform the needed normaliza- 
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tion, and then retry the operation. We have provided NORMAL, which gives nor- 
malizations adequate for a majority of applications. 

Non-Trapping NaN's 

The large number of representations for NaN's gives you the chance to put 
diagnostic information into a NaN. The information can be passed along as the NaN 
undergoes multiple floating point operations. The information can not only tell 
where the NaN came from, but also control further error action. 

EH87.L1B adopts the convention that the top bit of the fractional part of a NaN is a 
control bit, to be used in the following way: if the bit is 1, the NaN is to be con- 
sidered a "non-trapping NaN", for which no further error need be explicitly 
reported. You can return a non-trapping NaN as the result of an invalid operation. 
Then when the NaN passes through more arithmetic, there will be no more errors 
reported. You thus avoid multiple error messages which really come from only one 
error. 

The 8087 does not distinguish between trapping and non-trapping NaN's. All NaN's 
will generate an "I" error when they are used. However, you can provide an "I" 
trap handler which distinguishes between trapping and non-trapping NaN's. The 
procedure SIEVE does this for you. 



Non-Ordered Comparisons 



When you are testing two floating point numbers for equality, you may or may not 
want an error to be reported if those numbers are NaN's. The 8087 provides the 
FCOM and FTST instructions, which report an "I" error if they are given a NaN 
input. 

To suppress error reporting for NaN's in FCOM and FTST, we recommend the 
following convention: if either FCOM or FTST is followed by a MOV AX,AX 
instruction (8BC0 hex), then the "I" trap handler should treat non-trapping NaN's 
as legal inputs. It should return the answer NON-ORDERED (C3 = CO = 1), even if 
the two inputs are the same NaN, and act as if no "I" error had occurred. The pro- 
cedure SIEVE follows this suggested convention. 

Note that comparisons coded in PL/M-86 generate FCOM and FTST instructions 
which are not followed by a MOV AX, AX instruction. Therefore, according to the 
EH87.LIB convention, PL/M-86 comparisons of non-trapping NaN's are not con- 
sidered legal. There is no way to cause PL/M-86 to insert a MOV AX, AX instruc- 
tion after a comparison. 

The ESTATE87 Data Structure 

ESTATE87 is a 144-byte data structure created by DECODE, and used by the other 
EH87.LIB utility procedures. It contains most of the information your trap handler 
will need to provide customized error recovery: the state of the 8087, the identity of 
the offending operation, the values and formats of the operands, and possible 
already-calculated results. 

You will typically receive an ESTATE87 structure from DECODE, and pass it back 
to ENCODE mostly unchanged. You do not need to become familiar with those 
parts of ESTATE87 which you do not change. 

Following is a description of each of the fields of the structure ESTATE87. The off- 
sets mentioned are decimal numbers which give the number of bytes from the begin- 
ning of ESTATE87 to the beginning of the field. 
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OPERATION is a WORD, offset 0, which contains an error code identifying which 
procedure or 8087 instruction caused the error. The error codes for 8087 instructions 
are given in Appendix D. The error codes for CEL87 functions are the last two digits 
of the codes given in Appendix E. The error code for all DCON87 functions is 0C8 
hex. 

ARGUMENT is a BYTE, offset 2, which identifies the types and locations of the 
arguments of the interrupted operation. See figure 5-1 for the layout of the bit fields 
within the ARGUMENT byte. The 3-bit fields ATYPE1 and ATYPE2 indicate the 
types of ARG1 and ARG2, according to the following codes: 

no operand 

1 top of 8087 stack: ST(0) 

2 next element on 8087 stack: ST(1) 

3 the element of the 8087 stack specified by REGISTER 

4 a number in 8086 memory of a type given by FORMAT 

5 a TEMP_REAL operand 

6 a 64-bit integer operand 

7 a binary-coded decimal operand 

The PUSH ONCE bit is 1 if the result is pushed onto the 8087 stack (rather than 
replacing one or two of the input arguments). 

For example, ARGUMENT would equal 21 hex if the instruction FPREM caused 
the error, since the arguments to FPREM are the top two elements on the 8087 
stack, and the result replaces the inputs rather than causing the 8087 stack to be 
pushed. 

ARG1(5) is a WORD array, offset 3, which gives the first argument of the opera- 
tion, in the format specified by ARGUMENT. 

ARG1 FULL is a boolean BYTE, offset 13, whose bottom bit is 1 if ARG1 is pres- 
ent. If the error handler is called after the offending operation is done, the argument 
may no longer exist. 

ARG2(5) is a WORD array, offset 14, which gives the second argument of the 
operation, in the format specified by ARGUMENT. 

ARG2 FULL is a boolean BYTE, offset 24, whose bottom bit is 1 if ARG2 is pres- 
ent. If there is only one argument, or if the error handler is called after the offending 
operation is done, then this argument may not exist. 

RESULT is a BYTE, offset 25, which identifies the types and locations of the results 
of the interrupted operation. See figure 5-2. The 3-bit fields RTYPE1 and RTYPE2 
use the same codes as the corresponding fields in the ARGUMENT byte. The POP 
ONCE bit is set if the operation causes the 8087 stack to pop exactly once; the POP 
TWICE bit is set if the operation causes the 8087 stack to pop twice (i.e., the opera- 
tion is FCOMPP.) 



NOT 
USED 


1 1 

ATYPE2 

i 1 


PUSH 
ONCE 


1 1 

ATYPE1 

I I 



Figure 5-1 . Bit Fields of the ARGUMENT Byte in ESTATE87 121725-2 
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POP 
ONCE 


I I 

RTYPE1 
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Figure 5-2. Bit Fields of the RESULT Byte in ESTATE87 



121725-3 



5-3 



The Error Handler Module 8087 Support Library 



RES1(5) is a WORD array, offset 26, which gives the first result of the operation, in 
the format specified by RESULT. 

RES1 FULL is a boolean BYTE, offset 36, whose bottom bit is 1 if RES1 is pres- 
ent. If the error handler is called before the offending operation is done, the result 
does not yet exist. 

RES2(5) is a WORD array, offset 37, which gives the second result of the operation, 
in the format specified by RESULT. 

RES2 FULL is a boolean BYTE, offset 47, whose bottom bit is 1 if RES2 is pres- 
ent. If there is only one result, or if the error handler is called before the offending 
operation is done, then this result does not exist. 

FORMAT is a BYTE, offset 48, which specifies the memory data type when a field 
of ARGUMENT or REiSULT has value 4. (There is never more than one such field.) 
The possible values of FORMAT are as follows: 

forSHORT_REAL(32bits) 

1 for a 32-bit integer 

2 for LONG^REAL (64 bits) 

3 for a 16-bit integer 

REGISTER is a BYTEi, offset 49, which specifies the 8087 stack element number 
when a field of ARGUMENT or RESULT has the value 3. (No more than one such 
field can have this value.) The values of REGISTER range from for the top stack 
element to 7 from the bottom-most stack element. 

SAVE87(47) is a WORD array, offset 50, which contains the state of the 8087, as 
defined by the 8087 FSAVE instruction. Since DECODE is called after the 8087 
exceptions have been cleared, the 8087 status word stored into SAVE87 differs from 
the status word as it existed when the trap handler was invoked. This former value 
must be maintained separately from ESTATE87; it appears as the parameter 
ERRORS87 in all the EH87.LIB routines. 



How To Write an Exception Handler in ASM-86 
Using EH87.LIB 

By using EH87.LIB, you eliminate the difficult aspects of interfacing to the 8087 for 
your error recovery. However, there remains a strict protocol which must be fol- 
lowed. Following is a template of coding structure for an 8087 error handler written 
in ASM-86. If your customized error recovery is limited to errors which do not 
examine the ESTATE87 structure, you can follow the simpler template given under 
FILTER. 

1 . The handler procedure must be an INTERRUPT procedure. The 4-byte pointer 
is placed into the fixed slot for interrupt 16 (00040 hex). 

2. If you expect to return to the code which caused the handler to be called, you 
must preserve all 8086 registers. The 8086 flags are automatically saved and 
restored because this is an interrupt routine. 

3. The first floating point instruction of the handler should be an FNSTSW 
instruction, which stores the 8087 status word into a 16-bit memory location. 
This status word is the parameter ERRORS87 given to the EH87.LIB 
procedures. 

4. To insure synchronization with the 8087 at this point, there must be an arbitrary 
access to the 8086 memory space. We recommend a PUSH AX instruction 
followed by a POP AX instruction. 
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5. There should follow an FNCLEX instruction, to clear the 8087 exceptions. The 
exceptions must be cleared for the following EH87.L1B calls to work properly. 

6. The parameters to DECODE should be pushed onto the 8086 stack, and 
DECODE should be called. 

7. If you intend to use NORMAL and SIEVE, their calls should come immediately 
after the call to DECODE. You should not have intervening code that alters 
ESTATE87. NORMAL and SIEVE require ESTATE87 to be defined by 
DECODE. If ESTATE87 does not have values which could be output by 
DECODE, the results of NORMAL and SIEVE are undefined. 

8. If you have any customized code in your error handler, it should appear here. 

9. If the handler is returning to the calling environment, the parameters to 
ENCODE should be pushed onto the 8086 stack, and ENCODE should be 
called. 

10. The 8086 registers which were saved at the beginning of the exception handler 
should now be restored. 

11. The exception handler should be exited using an IRET instruction. A simple 
RET instruction will not do because the 8086 pushes its flags onto its stack when 
the interrupt is executed. 

You should also remember that if there is the possibility of a recursive call to the 
error handler (which can happen if ENCODED retries an instruction with exceptions 
unmasked), then all data storage used by your handler must be allocated on the 8086 
stack. This includes the ERRORS87 word, and the entire ESTATE87 structure. 



An Example of an Exception Handler 
Written in ASM-86 

NAME HANDLER_.87 

EXTRN DECODE: FAR 

EXTRN SIEVE: FAR 

EXTRN ENCODE: FAR 

EXTRN NORMAL: FAR 

INTERRUPT_TABLE SEGMENT AT 

ORG 4 * 16 

DD TRAP_HANDLER 
INTERRUPT_TABLE ENDS 

CODE SEGMENT PUBLIC 
ASSUME CS : CODE 



I MASK 



EQU 0001H 



STACK_LAY0UT STRUC 
OPERATION 



ARGUMENT 

ARG1 

ARG1 FULL 



DW ? 



DB ? 

DW 5 DUP(?) 

DB ? 



Position of "I" 
error bit and "I" 
mask 

Pointed at by BP 
during TRAP_HANDLER 

ESTATE87 template 
begi ns with this 
I i n e 
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AR62 


DW 


5 


DUP(?) 


ARG2 FULL 


DB 


? 




RESULT 


DB 


? 




RES1 


DW 


5 


DUP(?) 


RES1 FULL 


DB 


? 




RES2 


DW 


5 


DUP(?) 


RES2 FULL 


DB 


? 




FORMAT 


DB 


? 




REGISTER 


DB 


? 




CONTROL WORD 


DW 


? 





STATUS_WORD 

TAG_WORD 

ERROR_POINTERS 

STACK87 

RETRY_CONTROL 

RETRY_F LAG 

ERR0RS87 



DW ? 



DW 


? 




DW 


4 


DUP (?) 


DT 


8 


DUP (?) 


DW 


? 




DB 


? 




DW 


? 





Start of 94-byte 
FSAVE template 
As it exists after 
clearing exceptions 



Last line of FSAVE 
and of ESTATE87 
Retry 8087 Control 
setting for ENCODE 
Boolean parameter to 
ENCODE 

8087 Status Word 
before clearing 



You can place additional stack-allocated variables here, 

for your custom code. You refer to them as 

[BP] . your_var_name . If the number of bytes of inserted 

space is odd, you should eliminate the following dummy 

variable. 



DUMMY_EVEN_ALIGN 
REG ISTERS_86 
FLAGS_86 
RET_ADDRESS 
STACK LAYOUT ENDS 



DB ? 

DW 8 DUP (?) 

DW ? 

DD ? 



Not used; filler to 
keep SP even 
Pushed by 
PUSH_REG ISTERS 
Pushed by 8086 when 
interrupt called 
From TRAP HANDLER 



TRAP_HANDLER is a functioning 8087 exception handler 
written in ASM-86, using the procedures of EH87.LIB. It 
assumes that the only unmasked exceptions are "I" and 
"D", though we indicate below where you could insert 
code to handle other exceptions. 



TRAP_HANDLER PR0C FAR 

CALL PUSH_REGISTERS 

SUB SP, OFFSET REG ISTERS_86 

MOV BP , SP 

FNSTSW [BP3 . ERR0RS87 

PUSH AX 

POP AX 
FNCLEX 



This interrupt will 

preserve 8086 registers 

A I locate room for 

STACK_LAYOUT 

Set up indexing into 

STACK_LAYOUT 

Save the errors that 

caused the exception 

Arbitrary memory access 

to synch roni ze 

with the 8087 
Clear the except i ons so 
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CALL PUSH_IESTATE__ERRORS 

CALL DECODIE 

MOV AX, [BP1 .CONTROL_WORD 

MOV [BP] .RETRY_CONTROL,AX 

MOV [BP] .RETRY_FLAG,OFFH 

CALL PUSH_ESTATE_ERRORS 

CALL NORMAL 
RCR A L , 1 

JC ENCODE EXIT 



OR [BP] .RETRY_CONTROL,I_MASK 
CALL PUSH_ESTATE ERRORS 
CALL SIEVE 



RCR AL,1 

JC ENCODE EXIT 



TEST [BP] .ERR0RS87, I_MASK 
JNZ ACTUAL I ERROR 



hand L e r can use 8087 
Push parameters to 
DECODE onto stack 
ESTATE87 is now f i I Led 
with valid data 
Ex i st i ng cont ro L word is 
the def au It for 

any ENCODE retry 

attempt 
Default setting: retry 
will take place 
Push parameters to 
NORMAL onto stack 
Handle "D" errors 
Test boolean answer: was 
"D M the only error? 
If so then no other 
checking is necessary 

Mask "I " except i on 
Push parameters to SIEVE 
Check for lega I non- 
Trapping N a N 

Test boolean answer: was 
there such a NaN? 
If so then we can exit; 
"•I" has been cleared 

Was "I" error bit set? 



Here you could insert code to examine [ BP ] . ERR0RS87 to 

detect exceptions other than "I" or "D". If those 
other exceptions are detected, you can provide your 
own code to handle them. 

JMP ENCODE EXIT 



ACTUAL I ERROR: 



Here you may place your own customized code to deal 
with "I" exceptions other than legal non-trapping 
NaN's and denormalized inputs. The following lines set 
the " I" error bit, mask the "I" exception, and drop 
through to EN C0DE_EX I T . This simulates the masked "I" 
exception, but with non-trapping NaN's implemented. 
The user program must unmask the "I" exception when it 
tests and clears the "I" error bit. 



OR [BP] .STATUS_WORD, I_MASK 
OR CBP] .C0NTR0L_W0RD, I_MASK 

ENC0DE_EX IT : 

CALL PUSH_ESTATE_ERRORS 

PUSH [BP] .RETRY_C0NTR0L 

PUSH WORD PTR [ BP] . RETRY_F LAG 

CALL ENCODE 



Set the "I" error bi t 
Mask the "I" exception 



Push first two ENCODE 

pa ramete rs 

Push third ENCODE 

pa ramet er 

Push fourth ENCODE 

parameter 

Restore post-exc< t i on 
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ADD SP, OFFSET REG ISTERS_86 

CALL P0P_RE6ISTERS 

1RET 

TRAP HANDLER ENDP 



8087 state 

Release the STACK_LAY0UT 

storage 

Restore the eight 8086 

regi sters 

Restore 8086 flags; 

long-return to caller 



PUSH_ESTATE_ERRORS causes the two parameters 
ESTATE87_PTR and ERR0RS87 to be pushed onto the 8086 
stack, prior to a call to one of the EH87.LIB procedures 
which call for these parameters. 



PUSH_ESTATE_ERRORS PROC NEAR 
POP DX 
PUSH SS 

LEA AX, [BP] .OPERATION 

PUSH AX 

PUSH [BP] .ERR0RS87 

JMP DX 

PUSH ESTATE ERRORS ENDP 



save the return address 

push the segment half of 

ESTATE87_PTR 

push the offset half of 

ESTATE87_PTR 

--push complete 
ER"R0RS87 is the bottom 
half of this byte 
this is the RETURN from 
PUSH ESTATE ERRORS 



PUSH_REG I STERS causes the eight 8086 registers 
SI,DI ,ES,BP,DX,CX,BX,AX to be pushed onto the 8086 
stack. The registers can be popped off with a call to 
POP REGISTERS. 



PUSH_REG I STERS PROC NEAR 
PUSH DI 
PUSH ES 
PUSH BP 
PUSH DX 
PUSH CX 
PUSH BX 
PUSH AX 
MOV BP f SP 



XCHG SI, [BP+14] 
JMP SI 
PUSH REGISTERS ENDP 



Get stack pointer into an index 

regi ster 

Exchange the return address with 

SI which is being saved 

This is the RETURN of the 

procedure 



P0P_REG I STERS causes the eight registers which were 
pushed by PUS H_REG I STERS to be popped, restoring them to 
their original values. 



P0P_REGISTERS PROC NEAR 
POP SI 

MOV BP,SP 

XCHG SI, [BP + 14] 



Hold this call's return address 

temprar i ly in SI 

Get the stack pointer into an 

i ndex register 

Restore SI, and position the 

return address 
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POP AX 
POP BX 
POP CX 
POP DX 
POP BP 
POP ES 
POP DI 

RET ; Returnaddress is inpositiondue 

; to above XCHG 
POP_REG ISTERS ENDP 

CODE ENDS 

END 

How To Write an 8087 Exception Handler in 
PL/M-86 Using EH87. LIB 

Following is a template of coding structure for any 8087 error handler written in 
PL/M-86. Note that many of the protocols required for an ASM-86 error handler 
are not needed in PL/M-86. They are provided automatically by PL/M-86 and its 
built-in functions. 

If you customized error recovery is limited to errors which do not examine the 
ESTATE87 structure, you can follow the simpler template given under FILTER. 

1 . Your error handler procedure must have the attribute INTERRUPT 16. 

2. You must declare the structure ESTATE87, and the WORD variable 
ERRORS87. 

3. If you use any other variables in your error handler, they should also be 
declared within the procedure. All data must be declared within the procedure 
for the error handler to be reentrant. 

4. You should make the assignment ERRORS87 = GET$REAL$ERROR. 

5. You should call DECODE with the approriate parameters. 

6. If you intend to use NORMAL and SIEVE, their calls should come immeditely 
after the call to DECODE. NORMAL and SIEVE require ESTATE87 to be 
defined by DEiCODE. If ESTATE87 does not have values which could be out- 
put by DECODE, the results of NORMAL and SIEVE are undefined. 

7. If you have any customized code in your error handler, it should appear here. 

8. If the handler is returning to the calling environment, you should call ENCODE 
with the appropriate parameters. 

An Example of an 8087 Exception Handler Written in 
PL/M-86 

HANDLER_MODULE: DC- 
DECODE: PROCEDURE (ESTATE87_PTR,ERR0RS87) EXTERNAL; 

DECLARE (ESTATE87_PTR POINTER, ERR0RS87) WORD; 
END; 

ENCODE: PROCEDURE (ESTATE87_PTR , ERR0RS87 , C0NTR0L87 , 

RETRY_FLAG) EXTERNAL; 
DECLARE ESTATE87 PTR POINTER; 
DECLARE ERR0RS87 WORD; 
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DECLARE C0NTR0L87 WORD; 
DECLARE RETRY_FLA6 BYTE; 
END; 

NORMAL: PROCEDURE (ESTATE87_PTR , ERR0RS87) BYTE EXTERNAL; 

DECLARE (ESTATE87_PTR POINTER, ERR0RS87) WORD; 
END; 

SIEVE: PROCEDURE ( ESTATE87_PTR , ERR0RS87) BYTE EXTERNAL; 

DECLARE <ESTATE87_PTR POINTER, ERR0RS87) WORD; 
END; 

/* TRAP_HANDLER is a functioning 8087 exception handler 
written in PL/M-86, using the procedures of EH87.LIB. 
It assumes that the only unmasked exceptions are "I" 
and "D", though we indicate below where you could 
insert code to handle other exceptions. */ 

DECLARE I$ERROR$BIT LITERALLY • 1 H ' ; /* To set the "I" 

error bi t */ 
DECLARE TRUE LITERALLY ' F F H ' ; 

TRAP_HANDLER: PROCEDURE INTERRUPT 16; 

DECLARE ESTATE87 STRUCTURE ( 

OPERATION WORD, 

ARGUMENT BYTE, 

ARGK5) WORD, ARG1_FULL BYTE, 

ARG2(5) WORD, ARG2_FULL BYTE, 

RESULT BYTE, 

RESK5) WORD, RES1_FULL BYTE, 

RES2C5) WORD, RES2_FULL BYTE, 

FORMAT BYTE, 

REGISTER BYTE, 

CONTROL_WORD WORD, 

STATUS_WORD WORD, 

TAG_WORD WORD, 

ERR0R_P0INTERS(4) WORD, 

STACK_87(40) WORD); 
DECLARE ERR0RS87 BYTE; 
DECLARE C0NTR0L87 WORD; 
DECLARE RETRY_FLAG BYTE; 

ERR0RS87 = GET$REAL$ERROR; 
CALL DECODE O ESTATE87 , ERR0RS87 ); 
C0NTR0L87 = ESTATE87 . CONTROL_WORD; 
RETRY_FLAG = TRUE; 

IF NOT N0RMALOESTATE87,ERR0RS87) 
THEN DO; 
C0NTR0L87 = C0NTR0L87 OR I $ ER ROR$ B I T ; 
IF NOT SIEVEOESTATE87, ERR0RS87) 

THEN DO; 
IF ERR0RS87 
THEN DO; 
/* Here you may place your own customized code to 
deal with "I" exceptions other than non- 
trapping NaN's and denorma li z.ed inputs. The 
following lines set the "I" error bit, mask the 
"I" exception, and drop through to the RETRY 
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exit. This simulates the masked "I" exception, 
but with non-trapping NaN's implemented. The 
user program must unmask the "I" exception when 
it tests and clears the "I" error bit. */ 
ESTATE87.STATUS_W0RD = ESTATE87 . STATUS_WORD OR 

I$ERROR$BIT ■ 
ESTATE87. CONTROL_WORD = C0NTR0L87;' 
END; 
ELSE DO; 
/* Here you could insert code to examine ERR0RS87 
to detect exceptions other than "I" and "D". If 
those other exceptions are detected, you can 
provide your own code. */ 
END; 
END; 
END; 

CALL EN CODE OESTATE87 , ERR0RS87 , C0NTR0L87 , RETRY_F LAG); 

END TRAP_HANDLER; 

END HANDLER_MODULE; 
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DECODE Decode trapped operation and save 

8087 status 



Input: Two parameters, totalling six bytes, must be pushed onto the 8086 stack 
before calling DECODE. 

First, the four-byte pointer ESTATE87 PTR is pushed. As usual, the segment of 

ESTATE87^PTR is pushed first, followed by the offset. DECODE will send its 
output to the 144-byte memory buffer pointed to by ESTATE87 PTR. 

Second, a two-byte quantity whose bottom byte is ERRORS87 is pushed onto the 
stack. ERRORS87 is the bottom byte of the 8087 status word as it existed when the 
trap handler was called, before the exceptions were cleared. The top byte of the two- 
byte quantity is ignored. 



Function: DECODE provides, in a standardized format, all the information an 
exception handler might need to deal with the error condition. This includes the 
complete state of the 8087, the identity of the offending operation, and the 
existence, formats, and values of any arguments or results. 

We have programmed into DECODE the error calling specifications of all 8087 
instructions and all CEL87.LIB functions. Once DECODE has identified the inter- 
rupted operation, it uses these specifications to fill ESTATE87 with the correct 
arguments and results. 

An exception to the programming specifications occurs for the CEL87 functions 
that return values in 8086 registers. These functions are mqerIA2, mqerIA4, 
mqerlC2, mqerIC4, mqerIE2, and mqerlE4. The results are specified by DECODE 
to go to the 8087 stack, not to 8086 registers. 

DECODE identifies DCON87 errors by setting OPERATION to 0C8 hex, but no 
information about DCON87 arguments or results is available, because they reside in 
a location on the 8086 stack unknown to DECODE. 

Note that for FCOMP and FCOMPP instructions with denormal arguments, the 
trap handler is called with the inputs already popped off the 8087 stack. DECODE 
recovers these arguments from the bottom of the 8087 stack, and pushes them back 
onto the stack top. 



Output: All output is performed through the 144-byte memory buffer pointed at by 
the input parameter ESTATE87_PTR. The fields of the ESTATE87 structure are 
described in the section "The ESTATE87 Data Structure", at the beginning of this 
chapter. 

The 8087 itself is left completely cleared to its initialized state. This is the input state 
expected by the other procedures of EH87.LIB. 



PL/M-86 declaration and calling sequence:: 

DECODE: PROCEDURE ( EST AT E87_PT R , ERR0RS87) EXTERNAL; 

DECLARE ESTATE87_PTR POINTER, ERR0RS87 WORD; 
END; 

/* Assume ESTATE87 is already declared as presented 
earlier in this chapter */ 
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DECLARE ERR0RS87 WORD; 



/* Assume ERR0RS87 has been set to the value of the bottom 
part of the 8087 status word before exceptions were 
cleared. * / 

CALL DECODE ( 3ESTATE87, ERR0RS87 ); 



ASM-86 declaration and calling sequence: 

; the following line must occur outside of all SEGEMENT- 
; ENDS pairs 

EXTRN DECODE: FAR 



ESTATE_TEMPLATE STRUC 
DB 144 DUP(?) 



ESTATE TEMPLATE ENDS 



The individual fields of ESTATE87 
can be dec lared 
as in STACK_LAYOUT at the 
beginning of this chapter. 



DATA SEGMENT 

ESTATE_87 E ST AT E_T EMP L AT E 

ERR0RS87 DW ? 
DATA ENDS 

; the following code assumes that DS contains the segment 
; of ESTATE87 



ASSUME DS : DATA 
PUSH DS 

MOV AX,0I : FSET(ESTATE87) 
PUSH AX 

MOV AL,ERR0RS87 
PUSH AX 
CALL DECODE 



push the segment of 
ESTATE87 onto the stack 

4-byte pointer is now 
pushed onto stack 
second parameter 

-- pushed onto stack 
ESTATE87 is now fi lied 
with information 
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ENCODE Restore 8087 status and operation 

environment 



Input: ENCODE expects the 8087 to be completely cleared to its initialized state. 
Information about the 8087 is obtained from ESTATE87, not from the entry state 
of the 8087 itself. 

Four parameters, totalling ten bytes, must be pushed onto the 8086 stack before call- 
ing ENCODE. The first parameter is four bytes, and each of the other three 
parameters is two bytes. 

The first parameter, ESTATE87_PTR, points to the 144-byte ESTATE87 structure 
which was filled by a previous call to DECODE, and possibly modified by interven- 
ing code. ESTATE87 PTR is a four-byte pointer, with the segment pushed first 

and the offset pushed next. 

The second parameter pushed onto the stack is ERRORS87. ERRORS87 is a 
WORD parameter whose bottom byte is the bottom byte of the 8087 status word as 
it existed when the trap handler was called, before the exceptions were cleared. The 
top byte of ERRORS87 is ignored. 

The third parameter, CONTROL87, is a WORD parameter. It contains the value 
which the caller wants to place into the 8087 control word before retrying the func- 
tion when RETRY__FLAG is true. If RETRY_FLAG is false, CONTROL87 is 
ignored. 

The fourth and last parameter pushed onto the stack is RETRY FLAG. 

RETRY__FLAG is a BYTE parameter, the bottom half of the WORD which is 

pushed onto the 8086 stack (the top half is ignored). RETRY FLAG is a boolean 

control input which tells ENCODE what action to take after restoring the 8087 
environment, as described below. 

Function: ENCODE restores the 8087 to a state representing the completion of the 
operation which caused an error trap. This is accomplished by restoring the 8087 to 
the state stored in the ESTATE87 structure, after performing one of two actions, 
depending on RETRY __FLAG. 

If the bottom bit of RETRY FLAG is 0, the operation is assumed to have already 

been performed, and the results are assumed to be in ESTATE87. The results are 
copied to their destination. If insufficient results are available in ESTATE87, the 
value INDEFINITE may be used. 

If the bottom bit of RETRY FLAG is 1, the operation is identified from 

ESTATE87, and reperformed using operands as specified in ESTATE87. The 
parameter CONTROL87 is used as the 8087 control word for the retry; you can thus 
select which exceptions will be unmasked. After the retry is complete, the 8087 con- 
trol word as it exists in ESTATE87 is restored. 

If the operation being retried is a load or store that merely moves data without 
transforming it, ENCODE will perform the operation with exceptions masked, 
using the ARGUMENT and RESULT information of ESTATE87. This allows you 
to call ENCODE with true RETRY ___FLAG in all cases when NORMAL returns 
TRUE, without fear that the retry will repeat the "D" error. 

Note that ENCODE cannot be called with RETRY JLAG true if the error came 
from DCON87, since ESTATE87 does not contain enough information about 
DCON87. For the same reason, ENCODE cannot retry the CEL87 functions 
mqerIC2, mqerIC4, mqerIE2, or mqerIE4 after a "P" precision error. 
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ENCODE can perform the retry for ''I" errors coming from the above four CEL87 
functions, as well as mqerIA2 and mqerIA4, but the results will go to the 8087 stack 
and not to the correct 8086 registers. Your code must identify these functions and 
pop the answer from the 8087 to the correct destination in order to resume execution 
of the trapped program. 

If your CONTROL87 parameter contains unmasked exceptions, it is possible for an 
exception handler to be called during the retry, and for ENCODE to be invoked 
recursively. If you unmask the exception that caused the error, you should have 
modified the original arguments or results to avoid an infinite recursion. 



Output: ENCODE returns with the 8087 restored to the appropriate post-exception 
state, as indicated by ESTATE87 and the effects of either the retry or of the move- 
ment of results to their destination. Also, the input parameters are popped from the 
8086 stack. 



PL/M-86 declaration and calling sequence: 

ENCODE: PROCEDURE ( ESTATE87_PTR , ERR0RS87 , C0NTR0L87 , 
RETRY_F LAG) EXTERNAL; 

DECLARE ESTATE87_PTR POINTER; 

DECLARE ERR0RS87 WORD; 

DECLARE C0NTR0L87 WORD; 

DECLARE RETRY_FLAG BYTE; 
END; 

/* Assume ESTATE87 is already declared as presented 
earlier in this chapter. */ 

DECLARE ERR0RS87 WORD; 
DECLARE C0NTR0L87 WORD; 
DECLARE RETRY_FLAG BYTE; 

/* Assume that all the parameters have been filled with 
va lues . */ 

CALL ENCODE ( 3ESTATE87, ERR0RS87, C0NTR0L87, RETRY_FLAG ); 

/* The 8087 has now been returned to an appropriate post- 
exception state. * / 



ASM-86 declaration and calling sequence: 

; the following line must occur outside of all SEGEMENT- 
; ENDS pairs 

EXTRN ENCODE: FAR 

DATA SEGMENT 

ESTATE_87 ESTATE_TEMPLATE 

ERR0RS87 DW ? 
DATA ENDS 

The following code assumes that DS contains the segment 
of ESTATE87. It also assumes that the parameters to 
ENCODE have been set to appropriate values. 
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ASSUME DS : DATA 
PUSH DS 

MOV AX f OFFSET (ESTATE87) 

PUSH AX 

MOV AX,ERR0RS87 
PUSH AX 

MOV AX f C0NTR0L87 

PUSH AX 

MOV AL, RETRY_FLA6 

PUSH AX 

CALL ENCODE 



push the segment of 
ESTATE87 onto the stack 

4-byte pointer is now 
pushed onto stack 

second parameter is 
pushed onto stack 
load third parameter 
push onto stack 
last parameter 
-- pushed onto stack 



; the 8087 is now restored to an appropriate post- 
; exception state 
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FILTER Filter denormals and non-trapping NaN's from 

user error handling 



Input: The two-byte quantity whose bottom byte is ERRORS87 is pushed onto the 
8086 stack before calling FILTER. ERRORS87 is the bottom byte of the 8087 status 
word as it existed when the trap handler was called, before the exceptions were 
cleared. The top byte of the two-byte quantity is ignored. 

Function: FILTER provides a single interface to EH87.LIB which calls all of the 
other four functions DECODE, SIEVE, NORMAL, and ENCODE. If ERRORS87 
indicates that a "D" error caused the trap handler to be called, the denormal argu- 
ment is found and normalized. If an "I" error caused the trap, FILTER indicates 
whether the exception was caused by a legal non-trapping NaN. 

FILTER is a boolean function. It returns TRUE if the trap was indeed caused by one 
of the above-mentioned "I" or "D" cases. TRUE means that you may return im- 
mediately from the trap handler. FILTER returns FALSE if it has not handled the 
error. If FILTER returns FALSE, your trap handler should execute appropriate 
customized error recovery. 

Output: The boolean value of FILTER is output as the bottom bit of the AL 
register. The input word ERRORS87 is popped from the 8086 stack upon return. 

If FILTER returns FALSE, the 8087 is returned to the same state that it was prior to 
the call to FILTER. 

If FILTER returns TRUE, the 8087 is changed to reflect a retry of the instruction 
with normal operands (in case of "D" error) or "I" masked (in case of "I" error). 

PL/M-86 programming example: 

FILTER_MODULE: DC- 
FILTER: PROCEDURE (ERR0RS87) BYTE EXTERNAL; 

DECLARE ERR0RS87 WORD; 
END; 

/* SIMPLE_TRAP_HANDLER is an floating point exception 
handler which provides a bare minimum interface to 
EH87.LIB. If the error which caused this interrupt has 
been handled by FILTER, then the interrupt returns to 
the user program immediately. Otherwise, your inserted 
custom code is executed. */ 

SIMPLE_TRAP_HANDLER: PROCEDURE INTERRUPT 16; 

DECLARE ERR0RS87 WORD; 

IF FILTER (ERR0RS87 := GET$REAL$ERROR) THEN RETURN; 

/* Here you could place- code to handle exceptions other 
than the denormals and non- trapping NaN's handled by 
FILTER. */ 

END SIMPLE_TRAP_HANDLER; 

END FILTER_M0DULE; 
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ASM-86 programming example: 



NAME FILTER_MODULE 
EXTRN FILTER: FAR 

INTERRUPT_TABLE SEGMENT AT 

ORG 4 * 16 

DD SIMPLE_TRAP_HANDLER 
INTERRUPT_TABLE ENDS 

CODE SEGMENT PUBLIC 
ASSUME CS : CODE 

STACK LAYOUT STRUC 



P dinted at by BP during 
SIMPLE TRAP HANDLER 



You can place additional Itack-a I located variables here, 
for your custom code. You refer to them as 
BP.your var name. 



ERR0RS87 



DW ? 



REGISTERS_86 DW 8 DUP(?) 

FLAGS_86 DW ? 

RET_ADDRESS DD ? 

STACK LAYOUT ENDS 



8087 status word before 

clearing 

Pushed by PUSH_REG ISTERS 

Pushed by 8086 when interrupt 

called 

From SIMPLE TRAP HANDLER 



SIMPLE_TRAP_HANDLER is an floating point exception 
handler which provides a bare minimum interface to 
EH87.LIB. If the error which caused this interrupt has 
been handled by FILTER, then the interrupt returns to 
the user program immediately. Otherwise, your inserted 
custom code is executed. 



SIMPLE_TRAP_HANDLER PROC FAR 

CALL PUSH_REG ISTERS 

SUB SP, OFFSET REG ISTERS_86 

MOV BP,SP 

FNSTSW CBP] . ERR0RS87 

MOV AX, [BP] .ERR0RS87 
PUSH AX 

FNCLEX 

CALL FILTER 

RCR AL,1 

JC TRAP EXIT 



Save the 8086 registers 
Allocate room for stac k 
layout 

Set up indexing into 
STACK_LAYOUT 

Save the errors that caused 
the exception 
Fetch the status word ... 
.... and store it as FILTER 
parameter 
CI (?ar the exceptions so 
handler can use 8087 
Check for de norma I or non- 
trapping NaN 
Was error handled? 
If so then do nothing more 



Here you could place code to handle exceptions other 
than the denormals and non-trapping NaN's handled by 
FILTER. 
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TRAP_EXIT: 

CALL POP_REGISTERS 
IRET 



; Restore the 8086 registers 



SIMPLE_TRAP_HANDLER ENDP 

PUSH_REG I STERS causes the eight 8086 registers 
SI, DI,ES, BP f DX, CX,BX,AX to be pushed onto the 8086 
stack. The registers can be popped off with a call to 
POP_REGISTERS. 

PUSH_REGISTERS PROC NEAR 
PUSH DI 
PUSH ES 
PUSH BP 
PUSH DX 
PUSH CX 
PUSH BX 
PUSH AX 
MOV BP f SP 



XCHG SI, [BP + 14] 

JMP SI 
PUSH REGISTERS ENDP 



Get stack pointer into an index 

regi ster 

Exchange the return address with SI 

which is being saved 

This is the RETURN of the procedure 



P0P_R EGISTERS causes the eight registers which were 
pushed by PUSH_REG ISTERS to be popped, restoring them to 
their original values. 



POP REGISTERS PROC NEAR 

Hold this call's return address 

temporari ly in SI 

Get the stack pointer into an index 

regi ster 

Restore SI, and position the return 

address 



POP SI 

MOV BP,SP 

XCHG SI, [BP+14] 



POP AX 
POP BX 
POP CX 
POP DX 
POP BP 
POP ES 
POP DI 

RET ; Return address is in position due to above 
; XCHG 
P0P_REG ISTERS ENDP 

CODE ENDS 

END 
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NORMAL Detect and normalize "D" error non-normal 

arguments 



Input: NORMAL expects the 808Tto be completely cleared to its initialized state. 
Information about the 8087 is obtained from ESTATE87, not from the entry state 
of the 8087 itself. 

Two parameters, totalling six bytes, must be pushed onto the 8086 stack before call- 
ing NORMAL. 

First, the four-byte pointer ESTATE87 PTR is pushed. The segment of 

ESTATE87„PTR is pushed first, followed by the offset. ESTATE87JTR points 
to the 144-byte ESTATE87 structure which was filled by a previous call to 
DECODE. 

Second, a two-byte quantity whose bottom byte is ERRORS87 is pushed onto the 
stack. ERRORS87 is the bottom byte of the 8087 status word as it existed when the 
trap handler was called, before the exceptions were cleared. The top byte of the two- 
byte quantity is ignored. 



Function: NORMAL first checks ERRORS87 to see if the "D" error bit is the only 
unmasked error bit which is set. If this is not the case, NORMAL immediately 
returns FALSE, indicating that no normalization retry is to take place. 

If "D" is set, and no other unmasked error bits are set, NORMAL returns TRUE. If 
the operation that caused the "D" error was a load operation, the non-normal 
arguments are left unchanged. If the operation was not a load operation, NORMAL 

modifies the denormal arguments. Denormal SHORT REAL and LONG REAL 

arguments are normalized; denormal TEMP REAL arguments are converted to 

zero. Only the copies of arguments in ESTATE87's ARG1 and ARG2 fields are 
modified. 

Whenever NORMAL returns TRUE, you should call ENCODE with 
RETRY FLAG set to TRUE. You may leave "D" unmasked, using the same con- 
trol word in effect when the trap handler was called. If the operation was a load 
instruction, ENCODE will perform the load with "D" masked, and not cause a 
repeat of the "D" error. 

Note that NORMAL always returns FALSE if the operation that caused the trap is 
DCON87 or a CEL87 function. Only individual 8087 instructions can cause a "D" 
error. 



Output: NORMAL returns normalized arguments by modifying ESTATE87. In 
addition, it returns a boolean value in the All register indicating further action to be 
taken. If the bottom bit of AL is 1, then a normalization has taken place, and 
NORMAL is requesting the caller to retry the offending operation by calling 
ENCODE with a true RETRY__FLAG. If the bottom bit of AL is 0, no such retry is 
requested. 

NORMAL returns the 8087 itself to the same cleared state it had upon entry. 
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PL/M-86 declaration and calling sequence: 

NORMAL: PROCEDURE ( ESTATE87_PTR , ERR0RS87) BYTE EXTERNAL; 

DECLARE ESTATE87 PTR POINTER, ERR0RS87 WORD; 
END; 

/* Assume ESTATE87 is already declared as presented 
earlier in this chapter */ 

DECLARE ERR0RS87 WORD; 
DECLARE NORM__RETRY BYTE; 

NORM_RETRY = NORMALC 3ESTATE87, ERR0RS87 ); 

/* Now NORM_RETRY is true if a retry of the operation 
should be made */ ■ 



ASM-86 declaration and calling sequence: 

; the following line must occur outside of all SEGEMENT- 
; ENDS pairs 

EXTRN NORMAL: FAR 



DATA SEGMENT 
ESTATE_87 
ERR0RS87 
NORM_RETRY 

DATA ENDS 



ESTATE_TEMPLATE 
DW ? 
DB ? 



; the following code assumes that DS contains the segment 
; of ESTATE87 

ASSUME DS : DATA 

PUSHDS ;pushthe segment of ESTATE87 onto 

• the stack 
MOV AX,0FFSET(ESTATE87) 
PUSH AX 



MOV AL f ERR0RS87 

PUSH AX 

CALL NORMAL 

MOV NORM RETRY, AL 



4-byte pointer is now pushed onto 

stack 

second parameter 

-- pushed onto stack 
AL now tells whether to retry 
boolean answer stored in NORM RETRY 
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SIEVE Detect non-trapping NaN's which should 

be ignored 



Input: SIEVE expects the 8087 to be completely cleared to its initialized state. 
Information about the 8087 is obtained from ESTATE87, not from the entry state 
of the 8087 itself. 

Two parameters, totalling six bytes, must be pushed onto the 8086 stack before call- 
ing SIEVE. 

First, the four-byte pointer ESTATE87 PTR is pushed. The segment of 

ESTATE87__PTR is pushed first, followed by the offset. ESTATE87_PTR points 
to the 144-byte ESTATE87 structure which was filled by a previous call to 
DECODE. 

Second, a two-byte quantity whose bottom byte is ERRORS87 is pushed onto the 
stack. ERRORS87 is the bottom byte of the 8087 status word as it existed when the 
trap handler was called, before the exceptions were cleared. The top byte of the two- 
byte quantity is ignored. 



Function: SIEVE signals those cases in which the "I" exception should not have 
been given because the argument was a legal non-trapping NaN. This detection 
applies to all arithmetic operations which check for NaN inputs. 

SIEVE follows the conventions described in the sections titled "Non-Trapping 
NaN's" and "Non-Ordered Comparisons" towards the beginning of this chapter. If 
it is determined that this "I" error should not have taken place because the NaN 
arguments are non-trapping and the operation is not an ordered comparison, then 
SIEVE returns TRUE. In this case, the caller should retry the offending operation 
by calling ENCODE with a true RETRY_FLAG, and CONTROL87 modified so 
that "I" is masked. ENCODE will restore the original control word, with "F 
unmasked, after the retry is executed. 



kid 



If it is determined that this "I" error should still be flagged (for example, a stack 
error has occurred), SIEVE returns FALSE. Some appropriate "I" error recovery 
action should be performed. 

Note that if the operation which caused the trap is DCON87 or a CEL87 function, 
SIEVE will always return FALSE. All legal non-trapping NaN situations arise from 
individual 8087 instructions. 



Output: The boolean answer is returned as the bottom bit of the AL register. If this 
bit is 1, the "I" error should not have taken place, according to the non-trapping 
NaN conventions described earlier in this chapter. 

SIEVE leaves the 8087 in the same cleared state it had upon entry. The structure 
ESTATE87 is also unchanged by SIEVE. 



PL/M-86 declaration and calling sequence: 

SIEVE: PROCEDURE ( E S T AT E 87_PT R , ERR0RS87) BYTE EXTERNAL; 

DECLARE ESTATE87_PTR POINTER, ERR0RS87 WORD; 
END; 

/* Assume EST AT E87 is already declared as presented 
earlier in this chapter */ 



5-22 



SIEVE 



DECLARE ERR0RS87 WORD; 
DECLARE LEGAL_NAN BYTE; 

LE6AL_NAN = S I EVE ( 3ESTATE87, ERR0RS87 ); 

/* Now LEG A L_ N AN is true if there should not have been an 
■ ' I ' • error */ 



ASM-86 declaration and calling sequence: 

; the following line must occur outside of all SEGEMENT- 
; ENDS pairs 

EXTRN SIEVE: FAR 

DATA SEGMENT 

ESTATE_87 E ST ATE_TEMP LATE 

ERR0RS87 DW ? 
DATA ENDS 

; the following code assumes that DS contains the segment 
; of ESTATE87 

ASSUME DS : DATA 

PUSHDS ;pushthe segment of ESTATE87 onto 

• the stack 
MOV AX,0FFSET(ESTATE87) 
PUSH AX ; 4-byte pointer is now pushed onto 

; stack 
MOV AL f ERR0RS87 ; second parameter 
PUSH AX ; --pushed onto stack 

CALL SIEVE ; AL now contains boolean answer 
MOV LEGAL_NAN,AL ; LEGAL_NAN true if there is no 

: ' ' I ' ' error 
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Linkage of EH87.LIB to Your Program Modules 

This section tells how to link EH87.LIB into your program using a LINK86 
command. 

If you are linking both CEL87.LIB and EH87.LIB, it is absolutely necessary that 
EH87.LIB appear after CEL87.LIB in the LINK86 command for this reason: 
EH87.LIB contains references to all the functions of CEL87.L1B. However, if there 
are any unused CEL87 functions, they do not need to be linked for EH87.LIB to 
work correctly. We have therefore provided public symbols for CEL87.LIB in 
EH87.L1B, all of which point to an 8086 HLT instruction. CEL87.LIB must appear 
first in the LINK86 invocation to pick up the actual references to CEL87 functions; 
then EH87.LIB will supply HLT references (which do not add to the code size) for 
all unused CEL7 symbols. 

If you mistakenly put CEL87.LIB after EHS7.LIB in your link invocation, the link 
will be performed without any visible problems. But all your calls to CEL87.LIB will 
go to a HLT instruction. When you attempt to execute your program, the NDP will 
halt the first time a CEL87 function is called, 

EH87.LIB requires either the 8087 component or the 8087 emulator. If the com- 
ponent is present, you must link in 8087. LIB. If the emulator is used, you must 
instead link in E8087 and E8087.LIB. 

Following is the recommended ordering of modules in the LINK86 statement: 

Your object modules 

DCON87.LIB if your program uses it 

CEL87.LIB if your program uses it 

EH87.LIB 

8087. LIB if you are using the component, or 

E8087, E8087.LIB if you are using the emulator 



Examples: 

If your modules are MYMODl.OBJ and MYMOD2.0BJ and you are using the 
8087 emulator, issue the command 



LINK86 :F1 : MYM0D1 .OBJ, : F 1 : M Y M|0 D 2 . B J f &<cr> 

3 : FO : DC0N87 .LIB, : FO : CE L87 . L I|B , : FO : EH87 . L I B , &<cr> 

] :F0:E8087, : FO : E8087 .LIB TO I: F1 : MYPROG . LNK 



If you have a single module MYPROG. OBJ, are using the 8087 component, and do 
not use any libraries other than EH87.LIB, issue the command 



>> 



LINK86 : F1 : MYPROG . OB J f : FO : EH87 . L I B , &<cr> 
: FO : 8087 . L I B TO :F1 : MYPROG . LNK 
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CHAPTER 6 
IMPLEMENTING THE IEEE STANDARD 



The NDP, together with the 8087 Support Library, provides an implementation of 
"A Proposed Standard for Binary Floating Point Arithmetic", Draft 8.0 of IEEE 
TaskP754, Computer, March 1981, pp. 51-62. 

This chapter describes the relationship between the NDP and the IEEE Standard. It 
gives the choices we have made in the places where the Standard has options. It pin- 
points two areas in which we do not conform to the Standard. If your application 
requires it, you can provide software to meet the Standard in these areas. We 
indicate how to write this software. 

This chapter contains many terms with precise technical meanings, as specified by 
the Standard. We follow the Standard's convention of emphasizing the precision of 
our meaning by capitalizing those terms. The glossary in Appendix A provides the 
definitions for all capitalized phrases in this chapter. 



Options We Have Chosen 



Our SHORT_REAL and LONG_REAL formats conform perfectly to the 
Standard's Single and Double Floating Point Numbers, respectively. Our 

TEMP REAL format is the same as the Standard's Double Extended format. The 

Standard allows a choice of Bias in representing the exponent; we have chosen the 
Bias 16383 decimal. 

For the Double Extended format, the Standard contains an option for the meaning 
of the minimum exponent combined with a non-zero significand. The Bias for this 
special case can be either 16383, as in all the other cases, or 16382, making the 
smallest exponent equivalent to the second-smallest exponent. We have chosen the 
Bias 16382 for this case. This allows us to distinguish between Denormal numbers 
(integer part is zero, fraction is non-zero, Biased Exponent is 0) and Unnormal 
numbers of the same value (same as the denormal except the Biased Exponent is 1 .) 

The Standard endorses the support of only one extended format, and that is all we 
provide. 

The Standard allows us to specify which NaN's are Trapping and which are non- 
Trapping. EH87.LIB, which provides a software implementation of non-Trapping 
NaN's, defines the distinction. If the most significant bit of the fractional part of a 
NaN is 1 , the NaN is non-Trapping. If it is 0, the NaN is Trapping. 

When a masked "I" error involves two NaN inputs, the Standard lets us set the rule 
for which NaN is output. We choose the NaN whose absolute value is greatest. 



Areas Which Need the Support Library to 
Meet the Standard 

There are five features of the Standard which are not implemented directly on the 
8087. 

1 . The Standard recommends (but does not absolutely require) that a Normalizing 
Mode be provided, in which any non-normal operands to functions are 
automatically normalized before the function is performed. The NPX instead 
has a "D" exception, not mentioned in the Standard, which gives the exception 
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handler the opportunity to perform the normalization specified by the 
Standard. The "D" exception handler provided by EH87.LIB implements the 
Standard's Normalizing Mode completely for Single and Double precision 
arguments. Normalizing mode for Double Extended operands is implemented 
with one non-Standard feature, mentioned in the next section. 

2. The Standard specifies that the "=" comparison test yield an answer of FALSE, 
with no "I" error, if the relation between the input operands is "unordered". 
The 8087 FCOM and FTST instructions issue an "I" error for this case. The 
error handler EH87.LIB filters out the "I" error, using the following conven- 
tion: whenever an FCOM or FTST instruction is followed by a MOV AX, AX 
instruction (8BC0 hex), and neither argument is a trapping NaN, the error 
handler will assume that a Standard "=" comparison was intended, and return 
the correct answer with "I" erased. Note that "I" must be unmasked for this 
action to occur. 

3. The Standard requires that two kinds of NaN's be provided, Trapping and 
non-Trapping. Non-Trapping NaN's will not cause further "I" errors when 
they occur as operands to calculations. The NPX directly supports only Trapp- 
ing NaN's. The error handler module EH87.LIB implements non-Trapping 
NaN's by returning the correct answer with "I" erased. Note that "I" must be 
unmasked for this action to occur. 

4. The Standard requires that all functions which convert real numbers to integer 
formats automatically normalize the inputs if necessary. The 8087 FIST instruc- 
tion does not do so. CEL87.LIB's integer conversion functions fully meet the 
Standard in this aspect. 

5. The Standard specifies the remainder function which is provided by mqerRMD 
in CEL87.LIB. The 8087 FPREM instruction returns answers in a different 
range. 



Further Software Required to Fully 
Meet the Standard 

There are two cases in which you will need to provide further software to meet the 
Standard. We have not provided this software because we feel that a vast majority 
of applications will never encounter these cases. 

1. Non-Trapping NaN's are not implemented when "I" is masked. Likewise, the 
Standard's "=" function is not implemented when "I" is masked. You can 
simulate the Standard's concept of a masked "I" exception by unmasking the 
8087 "I" bit, and providing an "1" trap handler which supports non-Trapping 
NaN's and the "=" function, but otherwise acts just as if "I" were masked. 
Chapter 5 contains examples in both ASM-86 and PL/M-86 for doing this. 

2. Denormal operands in the TEMP REAL format are converted to by the 

EH87.LIB Normalizing Mode, giving sharp Underflow to 0. The Standard 
specifies that the operation be performed on the real numbers represented by the 
denormals, giving gradual Underflow. To correctly perform such arithmetic, 
you will have to normalize the operands into a format identical to 

TEMP REAL except for two extra exponent bits, then perform the operation 

on those numbers. Thus your software must manage the 17-bit exponent 
explicitly. 

We felt that it would be disadvantageous to most users to increase the size of the 
Normalizing routine by the amount necessary to provide this expanded arithmetic. 

Remember that the TEMP REAL exponent field is so much larger than the 

LONG REAL exponent field that it is extremely unlikely that you will encounter 

TEMP REAL Underflow. Remember further that the Standard does not require 
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Normalizing mode. If meeting the Standard is a more important criterion to you 
than the choice between Normalizing and warning modes, you can select warning 
mode ("D" masked), which fully meets the Standard. 

If you do wish to implement the TEMP REAL arithmetic with extra exponent bits, 

we list below some useful pointers about when the "D" error occurs: 

a. TEMP REAL numbers are considered Denormal by the NPX whenever the 

Biased Exponent is (minimum exponent). This is true even if the explicit 
integer bit of the significand is 1. Such numbers can occur as the result of 
Underflow. 

b. The 8087 FLD instruction can cause a "D" error if a number is being loaded 
from memory. It cannot cause a "D" error if it is coming from elsewhere in the 
8087 stack. 

c. The 8087 FCOM and FTST instructions will cause a "D" error for unnormal 
operands as well as denormal operands. 

d. In cases where both "D" and 'T" errors occur, you will want to know which is 
signalled first. When a comparison instruction is issued between a non-existent 
stack element and a denormal number in 8086 memory, the "D" and "I" errors 
are issued simultaneously. In all other situations, a stack "I" error takes 
precedence over a "D" error, and a "D" error takes precedence over a non- 
stack "I" error. 
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APPENDIX A 

GLOSSARY OF 8087 AND 

FLOATING-POINT TERMINOLOGY 



We continue in this appendix the convention of Chapter 6, capitalizing terms which 
have precise technical meanings. Such terms appear as entries in this glossary. Thus 
you may interpret any non-standard capitalization as a cross-reference. 

Affine Mode: a state of the 8087, selected in the 8087 Control Word, in which 
infinities are treated as having a sign. Thus the values +INFINITY and -INFINITY 
are considered different; and they can be compared with finite numbers and with 
each other. 

Base: (1) a term used in logarithms and exponentials. In both contexts, it is a 
number which is being raised to a power. The two equations (y = log base b of x) and 
(by = x) are the same. 

Base: (2) a number which defines the representation being used for a string of 
digits. Base 2 is the binary representation; Base 10 is the decimal representation; 
Base 16 is the hexadecimal representation. In each case, the Base is the factor of 
increased significance for each succeeding digit (working up from the bottom). 

Bias: the difference between the unsigned Integer which appears in the Exponent 
field of a Floating Point Number and the true Exponent that it represents. To obtain 
the true Exponent, you must subtract the Bias from the given Exponent. For 
example, the Short Real format has a Bias of 127 whenever the given Exponent is 
non-zero. If the 8-bit Exponent field contains 10000011, which is 131, the true 
Exponent is 131-127, or +4. 

Biased Exponent: the Exponent as it appears in a Floating Point Number, inter- 
preted as an unsigned, positive number. In the above example, 131 is the Biased 
Exponent. 

Binary Coded Decimal: a method of storing numbers which retains a base 10 
representation. Each decimal digit occupies 4 full bits (one hexadecimal digit). The 
hex values A through F (1010 through 1111) are not used. The 8087 supports a 
"Packed Decimal" format which consists of 9 bytes of Binary Coded Decimal (18 
decimal digits), and one sign byte. 

Binary Point: an entity just like a decimal point, except that it exists in binary 
numbers. Each binary digit to the right of the Binary Point is multiplied by an 
increasing negative power of two. 

C3 — CO: the four "condition code" bits of the 8087 Status Word. These bits are set 
to certain values by the compare, test, examine, and remainder functions of the 
8087. 

Characteristic: a term used for some non-Intel computers, meaning the Exponent 
field of a Floating Point Number. 

Chop: to set the fractional part of a real number to zero, yielding the nearest 
integer in the direction of zero. 

Control Word: a 16-bit 8087 register that the user can set, to determine the modes 
of computation the 8087 will use, and the error interrupts that will be enabled. 

Denormal: a special form of Floating Point Number, produced when an 
Underflow occurs. On the 8087, a Denormal is defined as a number with a Biased 
Exponent which is zero. By providing a Significand with leading zeros, the range of 
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possible negative Exponents can be extended by the number of bits in the Signifi- 
cand. Each leading zero is a bit of lost accuracy, so the extended Exponent range is 
obtained by reducing significance. 

Double Extended: the Standard's term for the 8087 Temporary Real format, with 
more Exponent and Significand bits than the Double (Long Real) format, and an 
explicit Integer bit in the Significand. 

Double Floating Point Number: the Standard's term for the 8087's 64-bit Long 
Real format. 

Environment: the 14 bytes of 8087 registers affected by the FSTENV and 
FLDENV instructions. It encompasses the entire state of the 8087, except for the 8 
Temporary Real numbers of the 8087 stack. Included are the Control Word, Status 
Word, Tag Word, and the instruction, opcode, and operand information provided 
by interrupts. 

Exception: any of the six error conditions (I, D, O, U, Z, P) signalled by the 8087. 

Exponent: (1) any power which is raised by an exponential function. For example, 
the operand to the function mqerEXP is an Exponent. The Integer operand to 
mqerY12 is an Exponent. 

Exponent: (2) the field of a Floating Point Number which indicates the magnitude 
of the number. This would fall under the above more general definition (1), except 
that a Bias sometimes needs to be subtracted to obtain the correct power. 

Floating Point Number: a sequence of data bytes which, when interpreted in a 
standardized way, represents a Real number. Floating Point Numbers are more ver- 
satile than Integer representations in two ways. First, they include fractions. Second, 
their Exponent parts allow a much wider range of magnitude than possible with 
fixed-length Integer representations. 

Gradual Underflow: a method of handling the Underflow error condition which 
minimizes the loss of accuracy in the result. If there is a Denormal number which 
represents the correct result, that Denormal is returned. Thus digits are lost only to 
the extent of denormalization. Most computers return zero when Underflow occurs, 
losing all significant digits. 

Implicit Integer Bit: a part of the Significand in the Short Real and Long Real for- 
mats which is not explicitly given. In these formats, the entire given Significand is 
considered to be to the right of the Binary Point. A single Implicit Integer Bit to the 
left of the Binary Point is always 1, except in one case. When the Exponent is the 
minimum (Biased Exponent is 0), the Implicit Integer Bit is 0. 

Indefinite: a special value that is returned by functions when the inputs are such 
that no other sensible answer is possible. For each Floating Point format there exists 
one non-Trapping NaN which is designated as the Indefinite value. For binary 
Integer formats, the negative number furthest from zero is often considered the 
Indefinite value. For the 8087 Packed Decimal format, the Indefinite value contains 
all 1 's in the sign byte and the uppermost digits byte. 

Infinity: a value which has greater magnitude than any Integer, or any Real 
number. The existence of Infinity is subject to heated philosophical debate. 
However, it is often useful to consider Infinity as another number, subject to special 
rules of arithmetic. All three Intel Floating Point formats provide representations 
for +INFINITY and -INFINITY. They support two ways of dealing with Infinity: 
Projective (unsigned) and Affine (signed). 
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Integer: a number (positive, negative, or zero) which is finite and has no fractional 
part. "Integer" can also mean the computer representation for such a number: a 
sequence of data bytes, interpreted in a standard way. It is perfectly reasonable for 
Integers to be represented in a Floating Point format; this is what the 8087 does 
whenever an Integer is pushed onto the 8087 stack. 

Invalid Operation: the error condition for the 8087 which covers all cases not 
covered by other errors. Included are 8087 stack overflow and underflow, NaN 
inputs, illegal infinite inputs, out-of-range inputs, and illegal unnormal inputs. 

Long Integer: an Integer format supported by the 8087 which consists of a 64-bit 
Two's Complement quantity. 

Long Real: a Floating Point Format supported by the 8087, which consists of a 
sign, an 11-bit Biased Exponent, an Implicit Integer Bit, and a 52-bit Significand; a 
total of 64 explicit bits. 

Mantissa: a term used for some non-Intel computers, meaning the Significand of a 
Floating Point Number. 

Masked: a term which applies to each of the six 8087 Exceptions I,D,Z,0,U,P. An 
exception is Masked if a corresponding bit in the 8087 Control Word is set to 1 . If an 
exception is Masked, the 8087 will not generate an interrupt when the error condi- 
tion occurs; it will instead provide its own error recovery. 

NaN: an abbreviation for "Not a Number"; a Floating Point quantity which does 
not represent any numeric or infinite quantity. NaN's should be returned by func- 
tions which encounter serious errors. If created during a sequence of calculations, 
they are transmitted to the final answer, and can contain information about where 
the error occurred. 

NDP: Numeric Data Processor. This is any iAPX86 or iAPX88 system that con- 
tains an 8087 or the full 8087 emulator. 

Non-Trapping NaN: a NaN in which the most significant bit of the fractional part 
of the Significand is 1. By convention, these NaN's can undergo certain operations 
without visible error. Non-Trapping NaN's are implemented for the 8087 via the 
software in EH87. LIB. 

Normal: the representation of a number in a Floating Point format in which the 
Significand has an Integer bit 1 (either explicit or Implicit). 

Normalizing Mode: a state, recommended by the Standard, in which non-normal 
inputs are automatically converted to normal inputs whenever they are used in 
arithmetic. Normalizing Mode is implemented for the 8087 via the software in 
EH87.LIB. 

NPX: Numeric Processor Extension. This is either the 8087 or the full 8087 
emulator. 

Overflow: an error condition in which the correct answer is finite, but has 
magnitude too great to be represented in the destination format. 

Packed Decimal: an Integer format supported by the 8087. A Packed Decimal 
number is a 10-byte quantity, with nine bytes of 18 Binary Coded Decimal digits, 
and one byte for the sign. 

Pop: to remove from a stack the last item that was placed on the stack. 
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Precision Control: an option, programmed through the 8087 Control Word, which 
allows all 8087 arithmetic to be performed with reduced precision. Since no speed 
advantage results from this option, its only use is for strict compatibility with the 
Standard, and with other computer systems. 

Precision Exception: an 8087 error condition which results when a calculation does 
not return an exact answer. This exception is usually Masked and ignored; it is used 
only in extremely critical applications, when the user must know if the results are 
exact. 

Projective Mode: a state of the 8087, selected in the 8087 Control Word, in which 
infinities are treated as not having a sign. Thus the values -(-INFINITY and 
-INFINITY are considered the same. Certain operations, such as comparison to 
finite numbers, are illegal in Projective Mode but legal in Affine Mode. Thus Pro- 
jective Mode gives you a greater degree of error control over infinite inputs. 

Pseudo Zero: a special value of the Temporary Real format. It is a number with a 
zero significand and an Exponent which is neither all zeroes or all ones. Pseudo- 
zeroes can come about as the result of multiplication of two Unnormal numbers; but 
they are very rare. 

Real: any finite value (negative, positive, or zero), which can be represented by a 
decimal expansion. The fractional part of the decimal expansion can contain an 
infinite number of digits. Reals can be represented as the points of a line marked off 
like a ruler. The term "Real" can also refer to a Floating Point Number which 
represents a Real value. 

Short Integer: an Integer format supported by the 8087 which consists of a 32-bit 
Two's Complement quantity. Short Integer is not the shortest 8087 Integer format- 
there is the 16-bit Word. Integer. 

Short Real: a Floating Point Format supported by the 8087, which consists of a 
sign, an 8-bit Biased Exponent, an Implicit Integer Bit, and a 23-bit Significand; a 
total of 32 explicit bits. 

Significand: the part of a Floating Point Number which consists of the most 
significant non-zero bits of the number, if the number were written out in an 
unlimited binary format. The Significand alone is considered to have a Binary Point 
after the first (possibly Implicit) bit; the Binary Point is then moved according to the 
value of the Exponent. 

Single Extended: a Floating Point format, required by the Standard, which pro- 
vides greater precision than Single; it also provides an explicit Integer Significand 
bit. The 8087's Temporary Real format meets the Single Extended requirement as 
well as the Double Extended requirement. 

Single Floating Point Number: the Standard's term for the 8087's 32-bit Short 
Real format. 

Standard: "a Proposed Standard for Binary Floating-Point Arithmetic," Draft 8.0 
of IEEE Task P754, Computer, March 1981, pp. 51-62. 

Status Word: A 16-bit 8087 register which can be manually set, but which is usually 
controlled by side effects to 8087 instructions. It contains condition codes, the 8087 
stack pointer, busy and interrupt bits, and error flags. 

Tag Word: a 16-bit 8087 register which is automatically maintained by the 8087. 
For each space in the 8087 stack, it tells if the space is occupied by a number; if so, it 
gives information about what kind of number. 
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Temporary Real: the main Floating Point Format used by the 8087. It consists of a 
sign, a 15-bit Biased Exponent, and a Significand with an explicit Integer bit and 63 
fractional-part bits. 

Transcendental: one of a class of functions for which polynomial formulas are 
always approximate, never exact for more than isolated values. The 8087 supports 
trigonometric, exponential, and logarithmic functions; all are Transcendental. 

Trapping NaN: a NaN which causes an "I" error whenever it enters into a calcula- 
tion or comparison, even a non-ordered comparison. 

Two's Complement: a method of representing Integers. If the uppermost bit is 0, 
the number is considered positive, with the value given by the rest of the bits. If the 
uppermost bit is 1, the number is negative, with the value obtained by subtracting 
(2 bit count ) from all the given bits. For example, the 8-bit number 1 1 1 1 1 100 is -4, ob- 
tained by subtracting 2 8 from 252. 

Unbiased Exponent: the true value that tells how far and in which direction to 
move the Binary Point of the Significand of a Floating Point Number. For example, 
if a Short Real Exponent is 131, we subtract the Bias 127 to obtain the Unbiased 
Exponent +4. Thus the Real number being represented is the Significand with the 
Binary Point shifted 4 bits to the right. 

Underflow: an error condition in which the correct answer is non-zero, but has a 
magnitude too small to be represented as a Normal number in the destination 
Floating Point format. The Standard specifies that an attempt be made to represent 
the number as a Denormal. 

Unmasked: a term which applies to each of the six 8087 Exceptions I,D,Z,0,U,P. 
An exception is Unmasked if a corresponding bit in the 8087 Control Word is set to 
0. If an exception is Unmasked, the 8087 will generate an interrupt when the error 
condition occurs. You can provide an interrupt routine that customizes your error 
recovery. 

Unnormal: a Temporary Real representation in which the explicit Integer bit of the 
Significand is zero, and the exponent is non-zero. We consider Unnormal numbers 
as distinct from Denormal numbers. 

Word Integer: an Integer format supported by both the 8086 and the 8087 that con- 
sists of a 16-bit Two's Complement quantity. 

Zerodivide: an error condition in which the inputs are finite, but the correct 
answer, even with an unlimited exponent, has infinite magnitude. 
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EXTERNAL SYMBOLS USED 

BY THE 8087 EMULATOR 



Intel translators put into generated object modules a number of special external 
symbols when the NPX is used by programs. These symbols are explained and listed 
in this appendix, so that you can recognize them in listings issued by LINK86 and 
LOC86. 

The machine code generated by an 8087 instruction for the 8087 component is 
usually the WAIT opcode, followed by an ESCAPE opcode, followed by further 
code which identifies the operands and/or the specific instruction. (See the 
8086/8087/8088 Macro Assembly Language Reference Manual for the details of 
these opcodes). 

The emulator can be called using the same amount of code. Instead of the WAIT 
and ESCAPE opcodes, a two-byte call to an 8086 interrupt routine is used. The 
interrupt routines are written to decode the remaining instruction bytes in the same 
way that the 8086/8087 decodes them. 

Thus, for each 8087 instruction, there are two opcode bytes which differ when the 
emulator is used instead of the chip. Each two-byte pattern for the 8087 chip is 
stored under a special name in 8087. LIB. The corresponding two-byte pattern for 
the emulator is stored under the same name in E8087.LIB. The special names are 
built into all Intel translators that support the 8087. Instead of issuing the opcode 
for an 8087 instruction, the translator issues the appropriate external name. LINK86 
then supplies the correct opcodes when the externals are satisfied, either by 
8087.LIBorE8087.LIB. 

87NULL.LIB satisfies the external names with the same opcodes as 8087. LIB. It is 
up to the user program to insure that any resulting 8087 instructions are never 
executed. 

Here are the special names. The colon given is part of the name; it exists to make it 
impossible for the name to be generated by any Intel translators. 



M: 


NCS 


M: 


NDS 


M: 


NES 


M: 


NSS 


M: 


NST 


M: 


WCS 


M: 


WDS 


M: 


WES 


M: 


WSS 


M: 


WST 


M: 


WT 
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POINT FORMATS 



The 8087 supports three formats for real numbers. All formats are in reverse order 
when stored in 8086 memory. That is, the first byte (lowest memory address) is the 
least significant byte of the Significand; the last byte (highest memory address) con- 
tains the sign and the top seven bits of the Exponent. 

The three formats supported are the 32-bit Short Real, the 64-bit Long Real, and the 
80-bit Temporary Real. 



Short Real: 

1 -bit Sign 

8-bit Exponent: 

0-bit Implicit Integer Bit: 

23-bit Fractional Part: 

Long Real: 

1-bit Sign 

11 -bit Exponent: 

0-bit Implicit Integer Bit: 

52-bit Fractional Part: 

Temporary Real: 

1 -bit Sign 
15-bit Exponent: 
1-bit Explicit Integer Biit 
63-bit Fractional Part: 



Bias 126 if Exponent zero; 127 if non-zero 
if Exponent zero; 1 if non-zero 
digits to the right of the Binary Point 



Bias 1022 if Exponent zero; 1023 if non-zero 
if Exponent zero; 1 if non-zero 
digits to the right of the Binary Point 



Bias 16382 if Exponent zero; 16383 if non-zero 
digits to the right of the Binary Point 



Special Values in All Three Formats: 

Infinity: maximum exponent and zero Significand 
NaN: maximum exponent and non-zero Significand 

Zero: minimum exponent and zero Significand 

Denormai: minimum exponent and non-zero Significand 
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SUMMARY OF 8087 INSTRUCTIONS 



On the following pages you will find a reference chart of NPX instructions. Here is 
an explanation of each column of the chart. 

Opcode: This column gives the machine codes generated by ASM86 and LINK86 
for the given instruction. Digits and upper-case letters "A" through "F" are hexa- 
decimal digits. In addition, the following special codes are used: 

• "i" denotes a 4-bit quantity whose top bit is 0, and whose bottom three bits give 
the 8087 stack element number for the instruction. 

• "j" denotes "i" plus 8. It is a 4-bit quantity whose top bit is 1, and whose 
bottom three bits give the 8087 stack element number for the instruction. 

• "/" followed by a digit denotes a MODRM byte, as described in the ASM86 
language manual. The digit gives the value of the middle REG field of the byte 
(bits 5,4,3). The MOD field (bits 7,6) can be any of the values 00, 01, or 10. The 
R/M field (bits 2,1,0) can be any of the 8 possible values. For some values of 
MOD and R/M, there will be one or two immediate displacement bytes follow- 
ing the MODRM byte, as defined by the 8086 architecture. 

The machine codes are those for the 8087 component, produced by linking 8087. LIB 
to your ASM-86 object modules. If there is a segment override byte, it goes between 
the first (WAIT or NOP) byte and the second (ESCAPE) byte. 

The code for the 8087 emulator, produced by linking E8087.LIB, differs as follows: 
The FWAIT instruction produces a NOP (90 hex) byte instead of the 8086 WAIT 
(9B hex) byte. All other instructions produce a first byte of CD hex, the 8086 INT 
instruction, instead of the WAIT or NOP byte produced for the 8087. If there is no 
segment override byte, the second byte's top hexadecimal digit is "1" instead of 
"D", giving an 8086 interrupt number between 18 hex and IF hex. The remaining 
parts of the instruction are the same. If there is a segment override byte, the byte is 
replaced by an interrupt number between 14 hex and 17 hex, as shown below. The 
third ESCAPE byte remains the same, and is interpreted by the emulator. 

Interrupt Override 

14 ES 

15 CS 

16 SS 

17 DS 

Instruction: This column gives the 8087 instruction just as it appears in an ASM-86 
program. Operands which begin with "mem" can be replaced with any memory 
operand denoting a data area of the correct number of bytes. The number following 
"mem" gives the decimal number of memory bytes acted upon. A trailing "r" in the 
name denotes a REAL format, "i" denotes an INTEGER format, and "d" denotes 
a Binary Coded Decimal format. 

Function: This column gives a concise description of the function performed by the 
instruction. 

Clocks: This column gives the typical number of clock cycles used by the 8087 chip 
to execute the instruction. It is not an exact number. If a MODRM byte is involved, 
a typical time of 9 cycles is added for calculating the effective address of the memory 
operand. 
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ErrlD: This column gives the hexadecimal value returned by the procedure 
DECODE in the library EH87.LIB, described in Chapter 5. The value indicates the 
type of instruction which caused an error, and is returned in the OPERATION field 
of the structure ESTATE87. 

Errors: This column lists the possible exceptions which can occur if the instruction 
is executed. 



Opcoc 


le 


Instruction 


Function 


Clocks 


ErrlD 


Errors 


9B D9 


FO 


F2XM1 


ST-(2**ST)-1 


500 


19 


UP 


9B D9 


E1 


FABS 


ST - | ST | 


14 


01 


I 


9B DE 


C1 


FADD 


ST(1)-ST(1l + ST, pop 


90 


05 


IDOUP 


9B DC 


Ci 


FADD ST(i),ST 


ST(i)-ST(i) + ST 


85 


05 


IDOUP 


9B D8 


Ci 


FADDST,ST(i) 


ST -ST(i) + ST 


85 


05 


IDOUP 


9B D8 


/0 


FADD mem4r 


ST -ST + mem4r 


114 


05 


IDOUP 


9B DC 


/0 


FADDmem8r 


ST -ST + m em8r 


119 


05 


IDOUP 


9B DE 


Ci 


FADDPST(i),ST 


ST(i) -ST(i) + ST, pop 


90 


05 


IDOUP 


9B DF 


IA 


FBLDmemlOd 


push, ST— memlOd 


309 


10 


I 


9B DF 


/6 


FBSTPmemlOd 


memlOd —ST, pop 


539 


OF 


I 


9B D9 


EO 


FCHS 


ST - - ST 


15 


02 


I 


9B DB 


E2 


FCLEX 


clear exceptions 


5 






9B D8 


D1 


FCOM 


compare ST- ST(1) 


45 


03 


ID 


9B D8 


Di 


FCOM ST(i) 


compare ST- ST(i) 


45 


03 


ID 


9B D8 


12 


FCOM mem4r 


compare ST - mem4r 


74 


03 


ID 


9B DC 


12 


FCOM mem8r 


compare ST - mem8r 


79 


03 


ID 


9B D8 


D9 


FCOMP 


compare ST - ST(1 ), pop 


47 


03 


ID 


9B D8 


Dj 


FCOMPST(i) 


compare ST - ST(i), pop 


47 


03 


ID 


9B D8 


/3 


FCOMP mem4r 


compare ST - mem4r , pop 


77 


03 


ID 


9B DC 


/3 


FCOMP mem8r 


compare ST - mem8r , pop 


81 


03 


ID 


9B DE 


D9 


FCOMPP 


compare ST - ST(1 ), pop 2 


50 


03 


ID 


9B D9 


F6 


FDECSTP 


decrement stack pointer 


9 






9B DB 


E1 


FDISI 


set interrupt mask 


5 






9B DE 


F1 


FDIV 


ST(1)-ST(1i/ST, pop 


202 


09 


IDZOUP 


9B DC 


Fj 


FDIVST(i),ST 


ST(i) - ST(i) / ST 


198 


09 


IDZOUP 


9B D8 


Fi 


FDIV St.ST(i) 


ST - ST / ST(i) 


198 


09 


IDZOUP 


9B D8 


/6 


FDIVmem4r 


ST -ST / mem4r 


229 


09 


IDZOUP 


9B DC 


/6 


FDIVmem8r 


ST -ST / mom8r 


234 


09 


IDZOUP 


9B DE 


Fj 


FDIVPST(i),ST 


ST(i) -ST(i) / ST, pop 


202 


09 


IDZOUP 


9B DE 


F1 


FDIVR 


ST(1)-ST/ST(1),pop 


203 


0A 


IDZOUP 


9B DC 


Fj 


FDIVRST(i),ST 


ST(i) -ST / ST(I) 


199 


0A 


IDZOUP 


9B D8 


Fj 


FDIVR ST, ST(i) 


ST -ST(i) / ST 


199 


0A 


IDZOUP 


9B D8 


11 


FDIVR mem4r 


ST-mem4r/ ST 


230 


0A 


IDZOUP 


9B DC 


11 


FDIVR mem8r 


ST-mem8r/ ST 


235 


0A 


IDZOUP 


9B DE 


Fi 


FDIVRPST(i),ST 


ST(i) -ST / ST(i), pop 


203 


0A 


IDZOUP 


9B DB 


EO 


FENI 


clear interrupt mask 


5 






9B DD 


Ci 


FFREEST(i) 


empty ST(i) 


11 






9B DE 


/O 


FIADD mem2i 


ST-ST + mem2i 


129 


05 


IDOP 


9B DA 


/0 


FIADDmem4i 


ST-ST + mem4i 


134 


05 


IDOP 


9B DE 


12 


FICOM mem2i 


compare ST - mem2i 


89 


03 


ID 


9B DA 


12 


FICOM mem4i 


compare ST - mem4i 


94 


03 


ID 


9B DE 


13 


FICOMPmem2i 


compare ST - mem2i , pop 


91 


03 


ID 


9B DA 


13 


FICOMPmem4i 


compare ST - mem4i , pop 


96 


03 


ID 


9B DE 


16 


FIDIVmem2i 


ST -ST - mem2i 


239 


09 


IDZOUP 


9B DA 


IB 


FIDIV mem4i 


ST -ST - mem4i 


245 


09 


IDZOUP 


9B DE 


11 


FIDIVRmem2i 


ST - mem2i / ST 


239 


0A 


IDZOUP 


9B DA 


11 


FIDIVR mem4i 


ST-mem4i / ST 


246 


0A 


IDZOUP 


9B DF 


10 


FILDmem2i 


push, ST — mem2i 


59 


10 


I 


9B DB 


10 


FILDmem4i 


push, ST — mem4i 


65 


10 


I 


9B DF 


/5 


FILDmem8i 


push, ST — mem8i 


73 


10 


I 


9B DE 


/1 


FIMULmem2i 


ST - ST * mem2i 


139 


08 


IDOP 


9B DA 


/1 


FIMULmem4i 


ST - ST * mem4i 


145 


08 


IDOP 


9B D9 


F7 


FINCSTP 


increment stack pointer 


9 






9B DB 


E3 


FINIT 


initialize 8087 


5 






9B DF 


12 


FIST mem2i 


mem2i — ST 


95 


OF 


IP 


9B DB 


12 


FISTmem4i 


mem4i — ST 


97 


OF 


IP 


9B DF 


13 


FISTPmem2i 


mem2i — ST, pop 


97 


OF 


IP 


9B DB 


13 


FISTP mem4i 


mem4i — ST, pop 


99 


OF 


IP 


9B DF 


11 


FISTPmem8i 


mem8i — ST, pop 


109 


OF 


IP 


9B DE 


IA 


FISUBmem2i 


ST-ST-mem2i 


129 


06 


IDOP 



D-2 
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Summary of 8087 Instructions 



Opcode 



Instruction 



Function 



Clocks ErrlD Errors 



9B DA IA 

9B DE /5 

9B DA 

9B D9 

9B DB 

9B D9 

9B DD 

9B D9 

9B D9 

9B D9 

9B D9 

9B D9 

9B D9 

9B D9 

9B D9 

9B D9 

9B DE C9 

9B DC Cj 

9B D8 

9B D8 

9B DC /1 

9B DE Cj 

90 

90 

DB E0 
DB E3 



90 
90 



/5 

Ci 

/5 

/0 

10 

E8 

/5 

IA 

EA 

E9 

EC 

ED 

EB 

EE 



Cj 

/1 



DB E2 
DB E1 



9B D9 

90 

90 



DO 
DD /6 
D9 17 



90 
90 

9B D9 
9B D9 
9B D9 
9B D9 



D9 /6 
DD 77 
F3 
F8 
F2 
FC 



9B DD IA 



9B DD 
9B D9 
9B D9 
9B DD Di 
9B D9 12 
9B DD 
9B D9 
9B D9 



/6 
FD 
FA 



12 
II 
16 



9B DD Dj 



9B DB 
9B D9 
9B DD 
9B DD 
9B DE 
9B DC Ej 
9B D8 Ei 
9B D8 IA 
9B DC IA 
9B DE Ej 
9B DE E1 
9B DC Ei 
9B D8 
9B D8 
9B DC 
9B DE 
9B D9 
9B 

9B D9 
9B D9 
9B D9 
9B D9 
9B D9 
9B D9 



Ej 
/5 
/5 
Ei 
E4 

E5 
C9 
Cj 
F4 
F1 
F9 



FISUBmem4i 

FISUBRmem2i 

FISUBRmem4i 

FLDST(i) 

FLDmemlOr 

FLD mem4r 

FLD mem8r 

FLD1 

FLDCW mem2i 

FLDENVmenn14 

FLDL2E 

FLDL2T 

FLDLG2 

FLDLN2 

FLDPi 

FLDZ 

FMUL 

FMULST(i),ST 

FMULST,ST(i) 

FMULmem4r 

FMUL mem8ir 

FMULPST(i),ST 

FNCLEX 

FNDISI 

FNENI 

FNINIT 

FNOP 

FNSAVE mem94 

FNSTCW mem2i 

FNSTENVmem14 

FNSTSWmem2i 

FPATAN 

FPREM 

FPTAN 

FRNDINT 

FRSTOR mem94 

FSAVEmem94 

FSCALE 

FSQRT 

FST ST(i) 

FSTmem4r 

FST mem8r 

FSTCW mem2i 

FSTENVmern14 

FSTPST(i) 

FSTPmemlOr 

FSTP mem4r 

FSTP mem8r 

FSTSW mem2i 

FSUB 

FSUBST(i),ST 

FSUBST,ST(i) 

FSUBmem4r 

FSUBmem8r 

FSUBPST(i),ST 

FSUBR 

FSUBRST(i),ST 

FSUBR ST,ST(i) 

FSUBR mem4r 

FSUBR mem8r 

FSUBRPST(i),ST 

FTST 

FWAIT 

FXAM 

FXCH 

FXCH ST(i) 

FXTRACT 

FYL2X 

FYL2XP1 



ST*-ST-mem4i 134 

ST *- mem2i - ST 129 

ST *- mem4i - ST 134 

push,ST*-oldST(i) 20 

push,ST*-mem10r 66 

push, ST *- mem4r 52 

push, ST*- mem8r 55 

push,ST*-1.0 18 

control word *- mem2i 19 

environment *- mem14 49 

push, ST *- log base 2 of e 18 

push, ST*- log base 2 of 10 19 

push, ST*- log base 10 of 2 21 

push, ST*- log base e of 2 20 

push, ST*- Pi 19 

push, ST*- +0.0 14 

ST(1)*-ST(1)*ST, pop 142 

ST(i) *- ST(i) * ST 138 

ST*-ST*ST(i) 138 

ST*-ST*mem4r 127 

ST*-ST*mem8r 170 

ST(i) *- ST(i) * ST, pop 142 

nowait clear exceptions 5 

nowaitset interrupt mask 5 

nowait clear interrupt mask 5 

nowait initialize 8087 5 

no operation 13 

nowait mem94 *- 8087 state 219 

nowait mem2i *- control word 24 

nowait mem14 *- environment 54 

nowait mem2i *- status word 24 

ST *- arctan(ST(1)°ST), pop 650 

ST*-REPEAT(ST-ST(1)) 125 

push,ST(1)/ST*-tan(oldST) 400 

ST*-round(ST) 45 

8087 state *- mem94 219 

mem94 *- 8087 state 219 

ST*-ST*2**ST(1) 35 

ST*- square root of ST 183 

ST(i)*-ST 18 

mem4r*-ST 96 

mem8r*-ST 109 

mem2i *- control word 24 

mem14 *- environment 54 

ST(i) *- ST, pop 20 

mem10r*- ST, pop 64 

mem4r*- ST, pop 98 

mem8r*- ST, pop 111 

mem2i *- status word 24 

ST(1)*-ST(1)-ST, pop 90 

ST(i) *- ST(i) - ST 85 

ST *- ST - ST(i) 85 

ST*-ST-mem4r 114 

ST*-ST-mem8r 119 

ST(i) *- ST(i) - ST, pop 90 

ST(1)*-ST-ST(1),pop 90 

ST(i) *- ST - ST(i) 87 

ST - ST(i) - ST 87 

ST*-mem4r-ST 114 

ST*-mem8r-ST 119 

ST(i) *- ST - ST(i), pop 90 

compare ST -0.0 42 
wait for 8087 ready 

C3-C0 *- type of ST 17 

exchange ST and ST(1) 12 

exchange ST and ST(i) 1 2 

push, ST(1) — expo, ST*- sig 50 

ST *- ST(1 ) * log2(ST), pop 950 

ST*-ST(1)*log2(ST+1),pop 850 



06 
07 
07 
10 
10 
10 



06 
06 
06 
06 
06 
06 
07 
07 
07 
07 
07 
07 
04 



0E 
0E 
0B 
1A 
1B 



IDOP 

IDOP 

IDOP 

I 

ID 

ID 



10 ID 

11 I 



13 
12 
15 
16 
14 
17 
08 
08 
08 
08 
08 
08 



I 

I 

IDOUP 

IDOUP 

IDOUP 

IDOUP 

IDOUP 

IDOUP 



1D UP 

1E IDU 

1C IP 

1F IP 



18 IOU 

0C IDP 

OF I 

OF IOUP 

OF IOUP 



OF I 

OF I 

OF IOUP 

OF IOUP 



IDOUP 
IDOUP 
IDOUP 
IDOUP 
IDOUP 
IDOUP 
IDOUP 
IDOUP 
IDOUP 
IDOUP 
IDOUP 
IDOUP 
ID 
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APPENDIX E 

SUMMARY OF SUPPORT 

LIBRARY FUNCTIONS 



The following pages give condensed summaries of the procedures and functions of 
the three libraries DCON87.LIB, CEL87.LIB, and EH87.LIB. You can use this 
appendix as a quick reference guide after you have assimilated the material in the 
main part of this manual. 

The 8087 Emulator is summarized in Appendix D, the table of 8087 instructions. 
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DCON87.LIB 



All input parameters to DCON87.LIB are 4-byte long pointers to 8086 memory buf- 
fers which contain the inputs and outputs. The pointer(s) are pushed onto the 8086 
stack, high byte first, before calling the procedure. The procedure returns with the 
pointer(s) popped off the 8086 stack. 

The first three procedures have one input pointer. The last four procedures have two 
input pointers. 



mqcBIN_DECLOW {Convert binary number to decimal string) 

Input: BIN_DECLOW_BLOCK_PTR - 

4-byte BIN PTR -*■ input binary number. 

1-byte BIN__TYPE: for SHORT__REAL 

1 for LONG__REAL 

2 for TEMP„REAL 
1-byte DEC_LENGTH: length of output field 

4-byte DEC PTR -* output decimal significand, decimal point at right 

2-byte DEC EXPONENT: base ten exponent of output 

1-byte DEC__SIGN: sign of output, in ASCII 

Sign and digits output for unusual inputs: 
NaN = ".." 
+INFINITY = "++" 
-INFINITY = " — " 
+0="0 " 
-0 = "-0" 

Errors: I,D,P 



mqcDEC_BIN (Convert decimal string to binary number) 

Input: DEC_BIN_BLOCK_PTR - 

4-byte BIN PTR -* output binary number. 

1-byte BIN__TYPE: for SHORT__REAL 

1 for LONG_REAL 

2 for TEMP_ REAL 
1-byte DEC_LENGTH: length of input field 
4-byte DEC_PTR — input string. 

Errors: 0,U,P 



mqcDECLOW_BIN (Convert decimal string, low-level interface, to binary 
number) 

Input: DECCLOW_BIN_BLOCK_PTR - 

4-byte BIN PTR -* output binary number. 

1-byte BIN_TYPE: for SHORT__REAL 

lforLONG_REAL 

2forTEMP_REAL 
1-byte DEC_LEN.GTH: length of input field 

4-byte DEC PTR -*■ input string, stripped down to decimal digits 

2-byte DEC EXPONENT: base ten exponent, with decimal point to 

right of input 
1-byte DEC_SIGN: sign of input, in ASCII 

Errors: 0,U,P 
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DCON87.LIB 



mqcLONG_TEMP (Convert LONG_REAL to TEMP_REAL) 

Inputs: LONG_REAL_PTR -* input number 
TEMP__REAL_PTR — output number 

Error: D 



mqcSHORT__TEMP (Convert SHORT„_REAL to TEMP_REAL) 

Inputs: SHORT__REAL__PTR -* input number 
TEMP__REAL_PTR -* output number 

Error: D 



mqcTEMP_LONG (Convert TEMP_REAL to LONG_REAL) 

Inputs: TEMP_REAL_PTR - input number 
LONG_REAL__PTR -*► output number 

Errors: I,0,U,P 



mqcTEMP_SHORT (Convert TEMP_REAL to SHORT_REAL) 

Inputs: TEMP_REAL__PTR -» input number 

SHORT__REAL__PTR -* output number 

Errors: I,Q,U,P 
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CEL87.LIB 



"x" denotes the 8087 stack top ST. 

"y" denotes the 8087 next stack element ST( I). 

"STK" denotes a 2-byte integer pushed onto the 8086 stack. 

All 8086 and 8087 stack inputs are popped on successful return. 

The errors columns give the hexadecimal code left in the 1 1-bit 8087 opcode register, 
along with the possible errors the function can produce. Unmasked errors will trap 
with the input(s) on the 8087 stack if under the "Inputs" column; with the output(s) 
on the 8087 stack if under the "Outputs" column. The first of the three hex digits 
tells how many numbers are on the 8087 stack when a trap handler is called. 









Errors, 


trap with: 


Name 


Function 


Inputs 


Outputs 


mqerACS 


x = arc cosine(x) 


175 


I 






mqerASN 


x = arcsine(x) 


174 


I 






mqerAT2 


x = arc tangent(y/x) 


277 


I 


277 


U 


mqerATN 


x = arc tangent(x) 


176 


I 






mqerCOS 


x = cosine(x) 


172 


I 






mqerCSH 


x = hyperbolic cosine(x) 


16F 


IO 






mqerDIM 


x = max(y-x,+0) 


265 


I 


265 


OU 


mqerEXP 


x = e** x 


16B 


IOU 






mqerlA2 


AX = roundaway(x) 


17E 


I 






mqerlA4 


DXAX = roundaway(x) 


168 


I 






mqerlAX 


x = roundaway(x) 


167 


I 






mqerlC2 


AX = chop(x) 


17E 


I 


17E 


P 


mqerlC4 


DXAX = chop(x) 


179 


I 


179 


P 


mqerlCX 


x = chop(x) 


166 


I 


166 


P 


mqerlE2 


AX = roundeven(x) 


180 


I 


180 


P 


mqerlE4 


DXAX = roundeven(x) 


17B 


I 


17B 


P 


mqerlEX 


x = roundeven(x) 


178 


I 


178 


P 


mqerLGD 


x = common log(x) 


16D 


IZ 






mqerLGE 


x = natural log(x) 


16C 


IZ 






mqerMAX 


x = max(x,y) 


282 


I 






mqerMIN 


x = min(x,y) 


281 


I 






mqerMOD 


x = (y mod x), has sign(y) 


269 


I 


269 


U 


mqerRMD 


x = (y mod x), close to 


27A 


I 


27A 


U 


mqerSGN 


x = (y with x'ssign) 


264 


I 






mqerSIN 


x = sine(x) 


171 


I 






mqerSNH 


x = hyperbolic sine(x) 


16E 


IO 






mqerTAN 


x = tangent(x) 


173 


IZ 






mqerTNH 


x = hyperbolic tangent(x) 


170 


I 






mqerY2X 


x = y ** x 


26A 


IZOU 






mqerYI2 


x = x** AX 


27C 


IOU 






mqerYI4 


x = x**DXAX 


27C 


IOU 






mqerYIS 


x = x**STK 


27C 


IOU 







" mqerDIM, mqerMAX, mqerMIN and mqerSGN can produce "D" errors from 
within their interiors. 
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EH87.LIB 



All EH87.LIB procedures operate on a structure ESTATE87, summarized below. 
Elements of ESTATE87: 



OPERATION WORD 
ARGUMENT BYTE 



ARG1(5)WORD 
ARG1_FULLBYTE 
ARG2(5)WORD 
ARG2_FULLBYTE 
RESULT BYTE 



RES1(5)WORD 
RES1_FULLBYTE 
RES2(5)WORD 
RES2_FULLBYTE 
FORMAT BYTE 



REGISTER BYTE 

CONTROI WORD WORD 

STATUS_WORDWORD 
TAG_WORDWORD 
ERROR_POINTERS(5) WORD 

STACK_87(40)WORD 



instruction or procedure which caused error 
two nibbles, each coded: 

0=no operand 

1=ST(0) 

2=ST(1) 

3=ST(REGISTER) 

4=see FORMAT 

5=TEMPREAL 

6=64-bit integer 

7=BCD 
bit 3 means push once 
value of bottom-nibble argument 
true if ARG1 contains a value 
value of top-nibble argument 
true if ARG2 contains a value 
formats of result, as in ARGUMENT; except 

bit 3 on means pop once, bit 7 on means 

pop twice 
value of bottom-nibble result 
true if RES1 contains a value 
value of top-nibble result 
true if RES2 contains a value 
format of type 4 ARGUMENT or RESULT: 

0=32-bit real 

1=32-bit integer 

2=64-bit real 

3=16-bit integer 
8087 stack register number for type 3 

ARGUMENT or RESULT 
8087 control word 
8087 status word 
8087 tag word 
8087 instruction pointer, opcode, operand 

pointer 
8087 stack of 8 temporary real values 



DECODE: Push 4-byte ESTATE87_PTR and 2-byte ERRORS87 before calling. 
DECODE fills ESTATE87 with information about the 8087 error that caused the 
exception. 



ENCODE: Push 4-byte ESTATE87_PTR, 2-byte ERRORS87, 2-byte 
CONTROL_WORD, 1-byte RETRY_FLAG before calling. ENCODE restores 
the 8087 to the state indicated by ESTATE87; if RETRY_FLAG is true it retries the 
error operation using CONTROL WORD. 



FILTER: Push 2-byte ERRORS87 before calling. FILTER calls DECODE, 
NORMAL, SIEVE, and ENCODE. FILTER returns AL TRUE if either NORMAL 
or SIEVE returned TRUE. 



NORMAL: Push 4-byte ESTATE87_PTR and 2-byte ERRORS87 before calling. 
NORMAL returns AL TRUE if "D" was the only error in ERRORS87; it also nor- 
malizes arguments if operation was not a load operation. 



SIEVE: Push 4-byte ESTATE87_PTR and 2-byte ERRORS87 before calling. 
SIEVE returns AL TRUE if there was a non-trapping NaN which should not have 
caused an "I" error. 
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APPENDIX F 

PUBLIC SYMBOLS IN THE 

SUPPORT LIBRARY 



Each of the three libraries DCON87.LIB, CEL87.LIB, and EH87.LIB contain 
public symbols other than the names of the procedures documented in this manual. 
They are internal names, used either within the libraries or by Intel translators. 

You should not use any of these names in your programs. 

There are no extra public names in the 8087 emulator or interface libraries, other 
than those listed in Appendix B. The names in Appendix B cannot be generated by 
Intel translators, so there is no possibility of conflict. 

Following is a list of all public names, both documented and undocumented, for 
each library. 



DCON87.LIB 

CHK_UNMSKD_0_U_ERR 

MQCBINDEC 

MQCBIN_DECLOW 

MQCDBX 

MQCDBXDB 

MQCDECBIN 

MQCDECBINLO 

MQCDECLOW_BIN 

MQCDEC_BIN 

MQCLONG_TEMP 

MQCSHORT_TEMP 



MQCSNGXDB 

MQCSNX 

MQCTEMP_LONG 

MQCTEMP_SHORT 

MQCXDB 

MQCXDBDB 

MQCXDBSNG 

MQCXSN 

POWER_OF_10 

UNMSKD__OV_OR__UN 

XCPTN_RTRN 



CEL87.LIB 

MQERACS 

MQERAIN 

MQERANT 

MQERASN 

MQERAT2 

MQERATN 

MQERCI2 

MQERCOS 

MQERCSH 

MQERDIM 

MQEREXP 

MQERIA2 

MQERIA4 

MQERIAX 

MQERIC2 

MQERIC4 

MQERICX 

MQERIE2 

MQERIE4 

MQERIEX 

MQERINT 

MQERIRT 

MQERLGD 

MQERLGE 

MQERMAX 



MQERMIN 

MQERMOD 

MQERNI2 

MQERNIN 

MQERRI2 

MQERRMD 

MQERRNT 

MQERSGN 

MQERSIN 

MQERSNH 

MQERTAN 

MQERTNH 

MQERY2X 

MQERYI2 

MQERYI4 

MQERYIS 

MQ_1 

MQ_2XM1 

MQ_63U 

MQ_63U1 

MQ_63UPI2 

MQ_AT2 

MQ_CONST 

MQ_COS 

MQ_CP2N63 



MQ_DECIDE 

MQ_EXIT 

MQ_EXM1 

MQ_I 

MQ_IRCHK 

MQ_LOG 

MQ_LOG10 

MQ_LOGDN 

MQ_MQRPI 

MQ_NAN 

MQ_NOF 

MQ_NORM 

MQ_NQ 

MQ_OF 

MQ_P0 

MQ_PI2 

MQ_PII 

MQ_Q 

MQ_RAD 

MQ__RERR 

MQ_SIN 

MQ_TXAM 

MQ_U0 

MQ_YL2X 
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Public Symbols in the Support Library 



8087 Support Library 



EH87.LIB 



DECODE 

ENCODE 

FILTER 

FQFORTRANSTATUSCHECK 

MQERACS 

MQERAIN 

MQERANT 

MQERASN 

MQERAT2 

MQERATN 

MQERCI2 

MQERCOS 

MQERCSH 

MQERDIM 

MQEREXP 

MQERINT 

MQERIRT 

MQERLGD 

MQERLGE 

MQERMAX 

MQERMIN 

MQERMOD 

MQERN12 

MQERNIN 

MQERR12 

MQERRMD 

MQERRNT 



MQERSGN 

MQERSIN 

MQERSNH 

MQERTAN 

MQERTNH 

MQERY2X 

MQERY14 

NORMAL 

SIEVE 

TQDECODE87 

TQENCODE87 

TQFETCH_AND_STORE 

TQINSTRUCTION_RETRY 

TQNANFILTER 

TQNORM87 

TQNORMALIZE 

TQPOP_THE_TOP 

TQREALMATHFILTER 

TQRESTORE_PTRS 

TQSAVE_PTRS 

TQUNPOP_THE_TOP 

TQ_312 

TQ_320 

TQ_322 

TQ_324 

TQ_326 
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INDEX 



8087. LIB, 2-3 
87NULL.LIB, 2-3, 2-4 

Accuracy of decimal conversions, 3-4 

ACS, 4-7 

Af fine mode, A-l 

Arc cosine, 4-7 

Arc sine, 4-9 

Arc tangent, 4-11,4-14 

ASN, 4-9 

AT2, 4-11 

ATN, 4-14 

Base, A-l 

Bias, A-l 

Biased exponent, A-l 

BIN_DECLOW, 3-7 

Binary Coded Decimal, A-l 

Binary point, A-l 

Binary to decimal conversion, 3-7 

Calls to PL/M-86 functions, 4-5 

CEL87.LIB,4-lff,E-4 

Characteristic, A-l 

Chop, A-l 

Chop functions, 4-26, 4-28, 4-30 

Chopping functions, list of, 4-1 

Common Elementary Function Library, 4-1 ff 

Common logarithm, 4-37 

Control word, 8087, 4-6, A-l 

COS, 4-16 

Cosine, 4-16 

CSH, 4-18 

D exception, 4-11, 5-1, 6-3 

DCON87.LIB, E-2 

DEC_BIN, 3-10 

Decimal conversion, 3-1 ff 

Decimal logarithm, 4-37 

Decimal to binary conversion, 3-10, 3-12 

Declarations, ASM-86, 3-1, 4-2 

Declarations, PL/M-86, 3-2, 4-3 

DECLOW_BIN, 3-12 

DECODE, 5-12 

Denormal, iii, A-l 

Difference, positive, 4-19 

DIM, 4-19 

Double extended format, 6-1, A-2 

Double floating point number, A-2 

E8087.LIB, 2-3 

EH87.LIB,E-5 

Emulator, 2-1 ff 

Emulator, symbols used, B-l 

ENCODE, 5-14 

Environment, A-2 

Error handler module, 5-1 ff 

Error reporting, 3-5, 4-6 

ESTATE87, 5-2, 5-12 



Exception, definition, A-2 
Exception handler, in ASM-86, 5-4 
Exception handler, in PL/M-86, 5-9 
Exceptions, emulation of, 2-1 
EXP, 4-21 
Exponent, A-2 
Exponential function, 4-21 
Exponential function, any base, 4-59 
Exponential functions, list of, 4-1 

FILTER, 5-17 
Floating point number, A-2 
Format of decimal numbers, 3-5 
Format, C-l 
FSTP-generated errors, 4-5 

Gradual underflow, A-2 

Hyperbolic cosine, 4-18 
Hyperbolic functions, list of, 4-1 
Hyperbolic sine, 4-53 
Hyperbolic tangent, 4-57 

IA2, 4-23 

IA4, 4-24 

iAPX convention, v, vi 

IAX, 4-25 

IC2, 4-26 

IC4, 4-28 

ICX, 4-30 

IE2, 4-32 

IE4, 4-34 

IEEE Standard, 3-4, 6-1 ff 

IEX, 4-35 

Implicit integer bit, A-2 

Indefinite, A-2 

Infinity format, C-l 

Infinity, A-2 

INIT87, 2-2, 2-3 

INITFP, 2-2, 2-3 

Initializing the 8087, 2-2 

Instruction address register, 8087, 4-6 

Instruction set, 8087, D-lff 

Integer power, 4-62, 4-64, 4-66 

Integer, A-3 

Interrupts, 2-2 

Invalid operation, A-3 

Inverse cosine, 4-7 

Inverse hyperbolic sine, example, 4-40 

Inverse sine, 4-9 

Inverse tangent, 4-11, 4-14 

LGD, 4-37 

LGE, 4-39 

Linkage, 2-3, 3-22, 4-69, 5-24 

Logarithm, common, 4-37 

Logarithm, natural, 4-39 

Logarithmic functions, list of, 4-1 

Long integer, A-3 
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8087 Support Library 



LONG_TEMP, 3-14 
Long real format, C-l 

Mantissa, A-3 

Manuals, related, v 

MAX, 4-41 

Maximum function, 4-41 

MIN, 4-43 

Minimum function, 4-43 

MOD, 4-45 

Modulus function, 4-45, 4-47 

NaN format, C-l 

NaN, A-3 

Natural logarithm, 4-39 

NDP, vi, A-3 

Need for Support Library, 1-1 

Nomenclature for 8086 family, v 

Non-ordered comparisons, 5-2 

Non-trapping NaN's, 5-2, 5-17, 5-20, A-3 

NORMAL, 5-20 

Normal distribution, example, 4-21 

Normal number, A-3 

Normalizing mode, 5-1, 5-17, 5-20, A-3 

Notational conventions, iv 

NPX, vi, A-3 

Opcode register, 8087, 4-6 
Overflow, A-3 

Packed decimal, A-3 

Partial emulator, 2-1 

PE8087, 2-1 

PL/M-86 function calls, 4-5 

Polar-to-rectangular conversion, example, 4-16, 4-51 

Pop, A-3 

Positive difference, 4-19 

Power function, 4-62, 4-64, 4-66 

Precision control, A-4 

Precision exception, A-4 

Projective mode, A-4 

Pseudo zero, v, A-4 

Public symbols, F-l 

Publications, related, v 

Real number, A-4 

Recovery of error information, 5-12 

Rectangular-to-polar conversion, 4-1 1 

Register usage, 4-6 

Remainder function, 4-45, 4-47 



Returning from exception handler, 5- 
RMD, 4-47 

Roundaway function, 4-23 thru 4-25 
Roundeven function, 4-32, 4-34, 4-35 
Rounding functions, list of, 4-1 

SGN, 4-49 

Short integer, A-4 

Short real format, C-l 

SHORT^TEMP, 3-16 

SIEVE, 5-22 

Sign transfer function, 4-49 

Significand, A-4 

SIN, 4-51 

Sine, 4-51 

Single extended format, A-4 

Single floating point number, A-4 

SNH,4-53 

Stack requirements, 2-2, 3-4, 4-5 

Standard, A-4 

Status word, A-4 

System software, 1-1 

Tag word, A-4 
TAN, 4-55 
Tangent, 4-55 
TEMP__LONG,3-18 
Temporary real format, C-l 
TEMP.__SHORT, 3-20 
TNH,4-57 
Transcendental, A-5 
Trapping Nan, A-5 
Trigonometric functions, list of, 4-1 
Truncation functions, list of, 4-1 
Twos complement, A-5 

Unbiased exponent, A-5 
Underflow, A-5 
Unmasked exception, A-5 
Unnormal, v, A-5 

Word integer, A-5 

Y2X, 4-59 
YI2, 4-62 
YI4, 4-64 
YIS, 4-66 

Zerodivide, A-5 
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